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Abstract 

An  integral  equation  approach  based  on  the  mixed-j>ot('ntial  forinulat ioti  i-  dovei- 
oped  for  the  analysis  of  radiation,  scattering  and  guidance  of  eh-ctromagnetic  \va\fs 
by  arbitrarily  shaped  conducting  objects  in  layered  uniaxial  media.  I  his  method 
successfully  applied  to  rigorously  analyze  coax-  and  waveguide-fed  microstrip  patch 
antennas  of  various  shapes,  as  well  as  multiconductor  t ransmissioti  littes  wjth  con¬ 
ductors  of  various  cross-section  profiles,  embedded  in  multilayered  media.  Willi  the 
approach  developed  here,  microstrip  structures  of  various,  pos.sibly  irregtilar  shajx's. 
embedded  in  multilayered  isotropic  or  uniaxial  substrates,  may  be  iitvestigated  within 
one  theoretical  framework  and  using  the  same  computer  [)rogram. 
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Preface 


This  report  comprises  five  technical  papers  written  by  the  principal  investigator  arn! 
his  collaborators.  The  unifying  tlieine  is  the  development  and  application  of  the 
mi.xed-potential  integral  equation  (MPIE)  formulation  to  the  analysis  (»f  radiation, 
scattering  and  guidance  of  electromagnetic  waves  by  conciucting  objects  t‘ml)eiided  in 
layered  material  media,  with  emphasis  on  non-planar  structures  hax  ing  irregular  or 
arbitrary  shapes.  The  papers  included  here  have  been  submitted  for  publication  in 
scientific  journals  and.  hence,  are  to  a  large  degree  .self-contained.  .As  a  result.  iIktc 
is  a  certain  amount  of  overlap  in  the  material  they  cover, 


Table  of  Contents 


1.  Formulation  of  Mixed-Potential  Integral  Elqiiations  for  Arbitrarily 

Shaped  Microstrip  Structures  With  I'niaxial  Substrates . 

2.  Mixed-Potential  Integral  Equation  (MPIE)  Formulation  for  Non-Pianar 

Microstrip  Structures  of  Arbitrary  Shape  . 

3.  RCS  Computation  of  Coax-Loaded  Microstrip  Patch  Antennas  of  Arbi¬ 
trary  Shape  . 

4.  Waveguide  Excited  Microstrip  Patch  Antenna — Theory  and  Experiment  . 

b.  Analysis  of  Multiconductor  Transmission  Lines  of  Arbitrary  Cross  Sec¬ 
tion  in  Multilayered  Uniaxial  Media . .  . 


Formulation  of  Mixed-Potential  Integral  Equations 
for  Arbitrarily  Shaped  Microstrip  Structures 
with  Uniaxial  Substrates 


Krzysztof  A.  Michalski 

Electromagnetics  <!^:  Microwave  Laboratory 
Department  of  Electrical  Engineering 
Texas  Ai:M  University 
College  Station.  Texas  77843-3128.  US.A 


Running  head'.  Formulation  of  Mixed-Potential  Integral  Equations 

Keywords:  Microstrip,  integral  equation,  layered  medium,  (Ireen's  function 

Abstract — A  systematic  development  is  presented  of  three  distinct  mixed-potential 
integral  equations  for  the  analysis  of  printed  circuits,  microstrip  resonators,  and  an¬ 
tennas  embedded  in  multi-layered,  planar,  uniaxial  media  of  infinite  lateral  ('xtent. 
with  emphasis  on  microstrip  structures  comprising  non-planar  conductors  of  irregular 
or  arbitrary  shapes. 

1  Introduction 

Most  of  the  analysis  methods  presently  available  for  the  analysis  of  microstrh>  struc¬ 
tures  are  either  limited  to  or  optimized  for  planar  geometries  with  regn',  r  shapes. 
However,  in  view  of  the  recent  advances  in  the  development  of  the  monolithic  mi¬ 
crowave  and  millimeter-wave  integrated  circuits,  it  is  increasingly  important  to  have 
at  one’s  disposal  techniques,  which  would  make  it  possible  to  acciu  ately  and  efhciently 
analyze  more  complex  microstrip  geometries,  comprising  no-. -planar  conductors.  In 
this  paper,  we  discuss  an  integral  equation  approach.  "  iiich  we  have  found  to  be 
especially  suitable  for  this  task.  More  specifically,  w'e  ievelop  three  distinct  forms  of 
the  so-called  mi.xed-potential  integral  equation  (MPIE)  [1],  [2]  for  arbitrarily  shaj^erl 
conducting  objects  embedded  in  a  planar,  multi-layered,  material  medium  of  infinite 


lateral  extent,  which  may  be  nmaxially  anisotropic.  Th<‘  (list iii}'ni.’'hiu!>  featiire  i»i  the 
.MPIE  formulations  is  that  they  employ  vector  and  scalar  potentials,  who  h  ais'  ex 
pressed,  respectively,  in  terms  of  the  current  and  charge  densities.  They  are  fuid'eraiiie 
to  other  forms  of  the  electric  field  integral  etpiation  (KFIF,  j.  because  ihev  invoi\i>  iesv 
singular  kernels  and  faster  convergent  s[)ectral  integrals  (or  series  i.  and  are  ainenabie 
to  the  well-established  method-of-moments  (M().\l)  solution  procedures.  oii”itiaii\ 
developed  for  arbitrarily  shaped  conducting  scalterers  in  fre«>  s|)ace  h$i.  :  Ij.  o  . 

.4n  MPIE  was  first  used  .n  the  present  context  by  Mosig  and  (l.udiol  dii.  who 
applied  it  to  analyze  planar  microstrip  patch  antennas  of  various  shapes  isee  also  jT’. 
[8]).  An  extension  o'"  this  approach  to  non-planar  conductors  in  layere<{  imulia  was 
later  proposed  by  .Michalski  [1].  More  recently,  Michalski  and  Zheng  have  d(>v<>loped 
three  alternative  MPIE  formulations  [9].  referred  to  as  Formulations  .\.  B,  and  C.  and 
applied  Formulation  C  to  analyze  scatterers  and  antennas  of  various  shapes  partially 
buried  in  a  material  half-space  [lOj.  as  well  as  microstrip  transmission  lines  having 
conductors  of  arbitrary  cross-section  (llj.  In  the  present  paper,  we  extend  the  MPIE 
approach  of  [9]  to  the  case  of  uniaxial  media,  which  have  important  ajtplications 
in  modern  microwave  circuits  [12],  [13],  .-Mthough  the  approach  pre.sented  here  iias 
successfully  been  applied  to  a  number  of  structures  in  layered  media  [llj.  [l  lj.  [luj. 
[16],  [17],  the  details  of  the  development  have  not  yet  appeared  in  the  open  literature. 

The  remainder  of  this  paper  is  organized  as  follows.  The  problem  statement  is 
given  in  Section  2.  In  Section  3,  we  introduce  the  Fourier  transform  apparatus  that 
considerably  simplifies  the  following  development,  and  express  the  electric  dyadic 
Green’s  function  of  the  layered  uniaxial  medium  in  terms  of  the  voltage  and  current 
Green's  functions  of  its  transmission  line  network  analog.  For  easy  reference,  the 
transmission  line  Green’s  functions  pertinent  to  a  multi-layered  uniaxial  medium  are 
given  in  the  appendix.  The  main  result  is  given  in  Section  1,  where  we  develop  and 
discuss  three  distinct  MPIE  forms  for  non-planar  conducting  structures  of  arbitrary 
shape  embedded  in  the  layered  medium,  using  the  dyadic  Green's  function  given  in 
Section  3  as  the  point  of  departure.  We  close,  in  Section  5.  with  the  summary  and 
conclusions. 
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2  Problem  Statement 


We  discuss  the  MPIE  approach  in  the  context  of  the  genera!  prol)l<‘m  illu>t rated  in 
Fig.  1.  which  shows  an  arbitrarily  shaped  cotidiictiiig  object  etn!j<*d(i*‘d  in  a  layered 
medium  and  excited  by  (prt'sumabiy  known)  <'lectric  aini  niagtH’tic  currents.  J'  and 
respectively.  1  he  material  layers  are  uniform  and  of  infinil<‘  extetit  aUiiig  the  ,r 

(2) 


Figure  1:  Arbitrarily  shaped  conducting  object  in  a  layered  medium. 

a  id  y  axes,  and  may  be  uniaxially  anisotropic,  with  the  optic  axis  parallel  to  the  c 
axis.  The  nth  layer  is  characterized  by  the  permittivity  and  permeability  dyadics 

C(n  T  ••  ytn  d"  f^zn  (  1  t 

where  is  the  unit  dyadic  transverse  to  c.  and  ctn  ifitn)  and  )  denote,  respec¬ 

tively,  the  transverse  and  longitudinal  dielectric  (magnetic.)  constants  relative  to  free 
space.  Observe  that  we  distinguish  dyadics  by  double  underlines  and  unit  vectors  in- 
carets.  The  free  .space  permeability  and  permittivity  will  be  denoted  by  /iq  and  (o. 
respectively.  For  each  layer,  we  also  introduce  its  electric  and  magnetic  anisotropy 
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ratios,  and  respectively,  given  as 

u,,  =  —  .  e,,  =  -  I J I 

^(n  Ihn 

The  top  and  bottom  layers  may  be  of  infinite  extent  along  the  ;  axis  or  as  illustrated 
in  Fig.  1 — may  be  shielded  by  ground  planes  having  specified  surface  admittances 
V5,.v+i  and  }T,i'  respectively,  where  the  latter  are  infinite  for  perfect  electrn  (oitduc 
tors. 

The  original  problem  of  fdg.  1  may  be  replaced  by  its  equivalent  ‘18.  |>.  lOtij  shown 
in  Fig.  2,  where  the  conducting  object  has  been  removed  and  its  effect  ra-placed  by  an 

{Z) 

T 

Z  — 

Z  =  2.V 

2  =  2^ 

2  =  2^ 

2  =  2/ 

Figure  2:  Problem  equivalent  to  that  in  Fig.  1. 

as  yet  unknown  electric  surface  current  with  density  Js,  residing  on  a  mathematical 
surface  S.  Because  this  current  exists  in  the  absence  of  the  conducting  object,  the 
resulting  electric  field  can  be  expressed  as  [19],  [20] 

E{r)  =  j  gf\T\T')- J{T')dr'  (3) 
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where  the  integration  is  over  the  region  occupied  by  J  the  volume  current  (]<'nsity 
associated  with  Js,  r  and  r'  denote  the  position  vectors  associated  witli  tlie  field  auri 
source  coordinates,  respectively,  and  wnere  Q^‘  is  the  electric  dyadic  (ireen's  {'unction 
of  the  layered  medium.  This  Green's  dyadic  will  be  developed  in  Section  d. 

Although  the  equivalent  current  in  Fig.  2  is  not  known,  we  may  easily  formidati' 
an  EFIE  that  J5  satisfies  by  imposing  the  impedance  boundary  condition  [21.  p.  18: 
on  the  surface  5  of  the  object.  Hence,  upon  using  (3).  we  obtain 


Zs  J  A„(r)'  Js{r)dS  -  j^An(r)-  J^gf{r\r')  •  Js{r')  dS'  d> 

=  j  An(r)-  E'{t)  dS  (1) 

where  {An}  is  a  suitable  set  of  vector  weight  functions  defined  over  and  tangential  to 
5.  E'  is  the  ‘incident’  electric  field  due  to  J'  and  M’  radiating  in  the  layered  medium 
of  Fig.  2.  and  where  Zs  is  the  surface  impedance  of  the  conducting  object,  given  as 


1+2 


8..  = 


(5) 


in  which  <7^  denotes  the  conductivity  and  Sc  the  skin  depth  of  the  conductor.  In  the 
case  of  strip-like  structures,  (5)  is  only  applicable  if  Sc  <<  tc.  where  f-  is  the  strip 
thickness.  If  this  condition  is  violated  the  expression 


may  be  more  appropriate  [22].  In  the  above  and  throughout  this  paper,  the  c-'"''  time 
variation  is  assumed. 

The  EFIE  (4)  can,  in  principle  at  least,  be  solved  for  Js  by  a  MOM  procedure, 
but  the  severe  singularity  of  its  kernel  as  r'  approaches  r  on  .8  makes  it  unsuitable 
for  this  purpose.  A  better  approach  is  based  on  a  mixed-potential  representation  of 
E,  gi  ven  as 

—  jB(r)  =  A(r) -f  V$(r)  (7) 
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in  which  the  vector  and  scalar  potentials,  .4  and  1>,  respectively.  ar<‘  ('xpressed 

A(r)  =  J Q;^{r\r')  ■  J(r')  dr'  (  S  | 

$(r)  =  J  C‘^{r\r')V'-J{r')dr' 

—  —  y  V'6"^(r|r')- J(r')  f/r'  (!)  i 

where  and  G'^  are  the  vector  and  scalar  potential  kernels,  and  where  V'  operates 
on  the  primed  (source)  coordinates.  In  view  of  the  divergence  theorem  [2:1.  p.  oO:!]. 
the  above  two  forms  of  $  are  equivalent  provided  n- J  vanishes  on  the  stirface  of 
the  object,  where  n  is  a  unit  vector  normal  to  5,  and — if  J  penetrates  a  boundary 
between  material  layers,  as  illustrated  in  Fig.  2 — provided  G^  is  a  continuous  function 
of  z'  across  the  interface.  Upon  using  (7)  in  conjunction  with  (8)  and  (9).  instead  of 
(3),  in  the  impedance  boundary  condition  on  5,  and  applying  the  divergence  theorem, 
we  arrive  at  the  MPIE 

Z5  j^An{ryJs{r)dS  +  I^Anir)-  J^g:^{r\r')-Js{r')dS'  dS 

-J^V-An{r)J^G^{r\r')V'-Jsir')dS'dS  =  J^An{ryE'(r)dS  (lO) 

provided  G*  is  a  continuous  function  of  r  across  any  interface  between  dissimilar 
material  layers  the  conducting  object  penetrates.  Observe  that  in  view  of  the  equation 
of  continuity,  V-J^  in  the  above  is  directly  proportional  to  the  surface  charge  density. 
Comparing  the  EFIE  (4)  with  the  MPIE  ( 10),  we  note  that  in  the  latter  the  differential 
operator  W'  has  in  effect  been  extracted  from  and  transferred  onto  An  and  Js. 
-As  a  consequence,  the  kernels  and  G®  appearing  in  (10)  are  less  singular  than 
G^,  and  the  Sommerfeld-type  integrals  that  arise  converge  more  rapidly. 

The  form  of  ( 10)  is  not  unique  and  leads  to  different  MPIE  formulations,  depending 
on  the  choice  of  and  G*^.  In  Section  4,  we  develop  three  distinct  MPIE  forms  that 
have  been  found  useful  in  MOM  analyses  of  the  radiation,  scattering,  and  guidance  of 
electromagnetic  fields  by  arbitrarily  shaped  conducting  structures  in  layered  media. 
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The  development  is  carried  out  in  the  spectral  domain  and  it  pnaceds  from  the 
spectral  electric  dyadic  Green's  function  of  the  layered  medium.  de\(‘ioped  in  tin- 
next  section. 


3  Preliminaries 


Following  Felsen  and  Marcuvitz  [21.  p.  188],  we  view  the  layered  meilium  in  Fiti.  2 
as  a  w-aveguide  along  the  c  axis  and  expand  the  electromagnetic  fields  and  currents 
in  terms  of  the  waveguide  modes.  Because  the  waveguide  cross  sectioti  is  of  infinite* 
extent  along  the  x  and  y  coordinates,  the  eigenvalues  corresponding  to  the  latter. 
kx  and  ky,  respectively,  form  continuous  spectra  on  the  real  axis.  ( ‘on.seqne'ntly.  the 
modal  expansions  may  be  expressed  as  the  familiar  double  Fourier  transforms,  which 
prompts  us  to  introduce  the  Fourier  transform  pair 


+  X.  +-X; 


x{/(p)} =/(*,)  = !  y /(pjeA-PA-rf., 


ill) 


X-‘{/(tu}s/(P)  = 


1 

CiTT)^ 


/  /  /(*,) 


—  x-  -  x- 


P  dkx  (Iky 


(12) 


where  kp^xkx-kyky  is  the  spectral  domain  counterpart  of  p  =  xx+yy.  If  we  now  write 
E{v)  as  E{p\z),  its  Fourier  transform  may  be  expressed  as  E{kp::).  and  similarly 
for  the  other  field  quantities.  VVe  also  note  that  in  the  spectral  domain  the  operator 
nabla  becomes  V  =  —jkp+zdjdz,  which  suggests  that  we  introduce  a  rotated  spectral 
domain  coordinate  system  based  on  kp,  as  illustrated  in  Fig.  3.  where  (cf.  [2o],  [26], 


1-271,  (281) 


V  = 


z  X  kp 

K 


K  =  ^Jk'i  +  q 


(13) 


When  the  Fourier-transformed  Maxwell’s  equations  are  projected  on  this  coordinate 
system,  it  is  found  that  the  spectral  domain  electromagnetic  field  may  efficiently  and 
elegantly  be  expressed  in  terms  of  the  voltages  and  currents  on  two  transmission  lines 
along  the  2  axis,  with  the  propagation  w'avenumbers  and  characteristic  admittances 


t 


8 


i^y) 


Figure  3:  Rotated  spectral  domain  coordinate  system. 


given  as 


and 


Z'  rjoki; 


r)oki)ft, 


( !H 


respectively,  where  ko  —  is  the  wavenumber  and  r/o  =  y/^oAo  the  intrinsic 

impedance  of  free  space.  Observe  that  at  this  stage  of  the  development,  the  per¬ 
meability  and  permittivity  (and  thus  also  the  anisotropy  ratios)  of  the  medium  may 
vary  arbitrarily  with  z  and  are  not  necessarily  restricted  to  be  piecewise-constant .  as 
is  the  case  in  Fig,  2.  The  voltage  and  current  on  the  line  with  the  pro|>aga- 
tion  wavenumber  and  characteristic  admittance  given  in  (It)  specify  the  transverse- 
magnetic  (TM)  partial  field,  while  the  voltage  V'^  and  current  on  the  line  with 
the  propagation  wavenumber  and  characteristic  admittance  given  in  ( I-’))  specify  flie 
transverse-electric  (TE)  part  of  the  field.  For  the  following  development,  it  will  lie 
convenient  to  also  introduce  transmission  line  Green's  functions,  which  represent  th'' 
voltage  or  current  excited  by  a  in  it-strength  voltage  or  current  point  source.  Hence, 
let  Vl{z\z')  and  If{z\z'),  where  the  superscript  p  stands  for  e  or  h.  denote  the  voltage 
and  current,  respectively,  at  a  point  c  due  to  a  1  .-V  current  source  located  at  r'  on 
the  corresponding  transmission  line  (cf.  [27],  [29]).  It  then  follows  that  these  Gree,i's 
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tuiu'tivms  arc  m»vcr!ic(!  by  the  coupieti  equalion" 


(Iz 

d  „  . . 

t  =  —  j  k\  \  \ !  -e  : 


where  <^1  r— c'l  is  t  he  Dirae  delta.  Similarly,  let  T"  t  1  anil  /'  <  :  vd  licriotc  t  iic  mm!  ,!uc 
and  cuna’iit.  respectively,  at  v  due  to  a  1  viil!ae<“  "^ottree  a!  It  !-■  then  iumid  tliet 
the.se  (ireen  s  Intiction'^  sattsiy  e<jiiaiioiis  du<d  to  i Ki)  i  ITt.  whn  h  are  ohtai.ned  fiom 

the  latter  bv  makinit  thc‘  snbstitntioiis  I'd'  — »  /f  I/,  . •  >  d  and  >  . ’ 

From  their  governing  etpsations.  it  can  also  be  sh<jwn  that  the  four  t  ran'.mts'.ioii  hue 
(Ireen's  tunctions  po.ssess  the  '^yimnetry  and  reciprocity  properties.  ,dl,  p  I'M' 


i;'’(cU')  =  /fide')  =  /fic'ic).  IdtcM'l 


i  1  's  1 


Once  the  transmission  line  (ireen's  fnnction.s  are  det<>rmined.  the  voitaae  and  current 
at  any  point  on  the  transrnis.sion  line,  excited  In-  an  arliitraiy.  distrilnited  voitaije 
and  current  may  be  expressed  in  terms  of  the  superposition  relations  (cf.  c{(). 
p.  63j.  [24.  p.  1!)3)) 

V^z)  =  I  V^iz\:')v^iz')d~’  ^  J  V^{z\z']r(z')dz'  i  19l 

r(z)  =  J  I^.(z\z')r^{z’jdz' d-  I  l^{z\z'}Zi:'\dz'  i20i 

In  view  of  the  uniformity  of  the  medium  of  Fig.  2  along  the  .r  and  t/  a.xes.  the 
dyadic  Green's  function  in  (3)  possesses  a  translational  symmetry  with  respi'ct  to 
the  transverse  coordinates.  Consequently,  we  may  write  Q/'dr\r'  ]  =  G^'  i  p  -  p':  ziz'l. 
where  Ofip]  may  be  given  as 


Of  ip-  c|r')  -  (0  •  c|;')  ■  g(Cl} 


Here,  is  the  spectral  domain  cotmterpart  of  expres.sed  in  the  rotated  coordi- 
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a) 


» 


I 


» 


» 


» 
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nate  system  oF  Fig.  :5,  and  tlie  dyadie  operator 


Q(^ )  =  ( ux  4-  vy )  eos  ^  -f  ( w y  —  ex )  s]  II  ^  f  i  i  i  22 1 

arui  its  transpose.  Q^is)-  ernployetl  tt)  tlirectlj  transloriii  from  tin'  r<^laied 

spectral  domain  coordinat*'  .system  to  tin*  space  domain.  Fpon  using  t  lie  •'Uperpo'! 

tion  relations  (IDl  l'iO)  in  the  Fourier-t ranst'ormed  .Ma.wvelFs  e(juat ions,  we  iind  th,!! 
'  F 

S:  may  be  expressed  in  terms  of  the  voltage  and  current  tran.smission  line  (ireenV 
functions  as 


cb')  =  uu  V;{ky.  c|c')  +  vv  K]\ky.  dc')  -  lii  ^  T' ab) 

4*.  .  d  ' 


■h, 


zu  ’■—  f^{ky.  rjc')  +  ii 
k'otz 


r  .. 


)> 

^  >iu 

<.?f; 


(2;Fi 


where  the  primed  and  unprimed  parameters  are  evaluated  at  z'  and  r,  respectively . 

a  convention  that  will  be  in  effect  throughout  this  paper,  unless  the  layer  index  is 
explicitly  stated.  Observe  that  in  the  above  we  explicitly  indicate  the  dependeine 
of  the  transmission  line  Green's  functions  on  the  transverse  spectral  wavenumbc.T  Ay. 
The  isotropic  medium  form  of  G  was  previously  given  by  Kastner  et  al.  [27]. 

In  view  of  the  fact  that  no  assumptions  have  so  far  been  made  regarding  the  d«'- 
pendence  of  the  media  parameters  on  c.  the  .spectral  dyadic  Green's  function  given 
in  (23)  is  valid,  in  particular,  for  the  layered  medium  of  Fig.  2.  where  the  perme¬ 
ability  and  permittivity  are  piecewise-constant.  The  corresponding  transmission  line 
network  is  illustrated  in  Fig.  4.  Observe  that  the  original  vector  problem  ha.‘^  thus 
been  reduced  to  a  much  simpler  equivalent  network  problem,  in  which  the  indivi<i- 
ual  material  layers  are  viewed  as  uniform  transmission  line  sections,  where  the  nth 
section  has  the  propagation  constant  ATn  and  characteristic  admittance  Y'/.  This  net¬ 
work  analog  actually  comprises  two  networks — one  associated  with  the  TM  and  the 
other  with  the  TE  partial  fields — which  have  identical  configurations,  but  differ  in  the 
characteristic  admittances  and  propagation  constants.  In  the  case  of  an  unshielded 
structure,  where  the  top  (or  bottom)  layer  is  of  infinite  extent  along  the  axis,  the 
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i 


i  1 


Ilz) 


Figure  i;  Transmission  line  network  analog  of  the  lay<‘red  medium. 


corresponding  transmission  line  section  in  the  network  analog  of  Fig.  1  also  cxiends 
to  infinity,  where  it  is  terminated  into  a  matched  load.  The  voltage  and  current 
Green's  functions  for  the  network  of  Fig.  4  are  easily  found  for  any  number  of  layers, 
as  discussed  in  the  appendix. 

We  note  that  in  view  of  (14)-(15),  the  transmission  line  Green’s  functions  appear¬ 
ing  in  (23)  depend  on  the  spectral  variables  Fj.  and  ky  exclusively  through  .As  a 
consequence,  when  the  double  spectral  integrations  in  (21)  are  performed  in  the  polar 
system  (kp,^)  illustrated  in  Fig.  3,  the  ^-integrals  can  be  evaluated  in  closed  form. 
More  specifically,  w’e  find  that 

where  we  have  introduced  the  notation 

«^n{/(^p)}  =  i^rdkpf{kp}k^^^Mkpp)  (25) 

In  the  above,  (p,  i^)  are  the  configuration  space  polar  coordinates  and  is  the  Bes.sel 
function  of  order  n.  In  this  fashion,  the  space  domain  dyadic  may  be 

expressed  in  terms  of  the  Sommerfeld-type  integrals  (25).  To  obtain  Q^{p  —  p'\  cjc') 
required  in  (3),  we  replace  ip  and  p  in  (24)-(25)  by  \)  and  p,  respectively,  where 

o  =  \/{x-x'Y  +  (y-!/)^  ( 2fi ) 


11 


i2 


Since  only  will  be  needed  in  I  lie  folKiwing  (.i<*vek>pnient .  in  the  intt‘rest  u!  breviis 
we  do  not  list  the  explicit  form  of  its  space  domain  coiinter]>att . 


4  Development  of  the  MPIE  Formulations 


Because  (3)  has  the  form  of  a  convolution  in  the  transv«‘rs<‘  plane,  its 

f  ourier  I  ransliirm 

becomes 

Eik:  z)  =  I a^'ikp:  zk'y  c')  <h' 

!  27  i 

with  given  in  (23).  Similarly,  the  spectral  domain  counterpart 

of  (7)  is  found  as 

-  Eikp-,:)  ^  A{kp::^  +  V^{kp:  :) 

(281 

in  which 

A{kp:z)  =  1  Q\kp;z\zyj{kp-.=')dz' 

(29) 

Mkp-.z)==  -  lv'G^{kp:z\zyj(kp:z')dz' 

(30) 

MB  ^ 

where  V  =  jkp  +  zdjdz'.  Comparing  (27)  and  (28).  we  arrive  at  the  relation 

-  ^'1~')  -  -I  -  vv  (7"(Jfe,:  -'ir')  (31 ) 


which  will  serve  as  the  point  of  departure  in  the  development  of  and  from 
In  what  follows,  we  continue  to  use  the  rotated  spectral  domain  coordinate  system  of 
Fig.  3,  in  w'hich  four  of  the  nine  components  of  are  zero,  as  is  evident  from  (23). 
This  and  the  fact  that  V  and  V  are  orthogonal  to  the  u  axis  greatly  simplifies  the 
development.  Hence,  we  immediately  find  from  (31)  and  (23)  that 


Ot  =  Ot  =  Gi=Gi  =  i) 


and  that  the  nonzero  elements  of  and  G"*  are  related  bv 


r'E  _  r'A  L-ir'o 


(32) 


(33) 
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For  notatioiial  simplicity,  we  omit  in  the  above  and  in  some  (■([uations  iliai  iuilow 
the  arguments  of  the  spectral  domain  kernels.  Observe  that  the  system  (ddi  (d7 i  is 
underdetermined,  because  there  are  si.x  unknowns,  viz.  the  live  components  of 

=  wwO’uu  +  +  uzG^.  +  (dS) 

and  and  only  five  eciuations.  Consequently,  we  are  free  to  impose  an  additional 
constraint  on  (38) — for  instance,  we  may  choose  one  of  the  elements  of  to  be  zero. 
Note,  however,  that  there  is  no  flexibility  in  choosing  G'X-  immediately 

follows  from  (34)  and  (23)  that 

=  (39) 

Below,  we  introduce  three  different  constraints  on  the  remaining  four  elements  of 
and  explore  the  resulting  MPIE  formulations,  w^hich  will  be  referred  to  as  Formula- 
liens  A.  B.  and  C.  following  the  terminology  introduced  in  [9]. 


4.1  Formulation  A 

To  arrive  at  Formulation  A,  we  set 


( 10) 


It  then  follows  from  (^Ki)  that 


i)(r 

1)7 


J  'Ai 


/; 


<  11 


which,  in  view  of  (  Hi),  leads  to 


//'  (A-:  )■ 


provided  the  field  point  is  within  a  homogeneous  material  layer.  I’pon  u-ing  the 
above  in  (33).  we  arrive  at 


/•■t  /  )  .  ^  -f  ,  , 


Similarly,  from  (42)  and  (35).  we  find 


A'of'  d:' 

which,  in  view  of  (16)  and  (18).  may  be  reduced  to 


6''^,(Ap;  rlc')  =  A-or/o/i'  (  1 


Finally,  from  (37)  and  (42),  we  obtain 


/,.4  fVo^p\  7  ,  no 


+ 


^(r-r')  + 


1  on- 


ko  J  t,t',  '  jkoY  ''  '  ’  dzdz' 

which,  upon  using  (16),  (18),  and  the  equation  dual  to  (17).  simplifies  to 


(  I3i 


(44) 


(45) 


(Hi) 


(17) 


We  have  thus  determined  ail  elements  o[ Q"'  in  (38),  as  well  as  the  corresponding 
To  obtain  the  space  domain  counterpart  of  we  subject  the  latter  to  the  inverse 
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Fourier  transform  in  the  manner  of  (‘Jl)  aiul  use  (21),  vvitli  the  re-^uh 


g'(p;r|r')  -  lLSolvSl^y-=\=')  +  \  ;'  (Av:cir') 


(A- 


+  -[(ii-yy)  cos  'l^  +  {xy  +  yx )  sin 2,^]  5, 


Htfz 


I  ,^(  A;,:  c;c' I  A-,';//.(,  \/(A;,:cc'i 


[iicos,^  4- yisin,^)  jA-ot/o/M  1 - TT 


iliikiiiilll 

i/MA-':)' 


A-’ 


+  zz  — 5o{/'(A>;  r|r')} 

t  - 


IS! 


Similarly,  the  space  domain  counterpart  of  Ci'^  is  readily  found  as 


6'"(p;r|c') 


Wi£!£)l 

J 


(Fi) 


Finally,  in  order  to  obtain  Qi^ip  —  p':  :\~)  and  C^ip- p':  z\:')  required  in  the  MPIE 
(10),  we  make  the  replacements  sp  i)  and  p  — o  in  (48)  and  (49).  with  i)  and  g 
given  in  (26).  Observe  that  when  the  source  and  field  points  are  within  the  same 
material  layer.  S'*  becomes  a  diagonal  dyadic  and  the  x:  and  yz  components  of 
vanish,  thus  resulting  in  a  significant  simplification  of  this  MPlF  formulation. 
Note,  however,  that  in  this  formulation  is  a  discontinuous  function  of  r  across  any 
interface  between  dissimilar  material  layers. 


4.2  Formulation  B 

To  obtain  Formulation  B,  we  set 

It  then  follows  from  (35)  that 

dG^  jrjo 


16 


which,  in  view  of  (16),  leads  to 


I  'fit-  ■ 

,/x  »,  I V.  - 


=1,')  =  ^ 


t  )-'i 


provicied  the  source  point  is  within  a  homogeneous  material  layer.  Tpon  U'-iiig  liie 
above  in  (33),  we  arrive  at 


Similarly,  from  (52)  and  (36).  we  find 


fcoCz  ^  dz 

which,  in  view  of  (16)  and  (18).  may  be  reduced  to 

kpl^(ky. 


c);')  =  koT|o^lt  ( 1  - 


Finally,  from  (52)  and  (37),  we  obtain 


I/"')  A-:')' 


A’O  / 


-z')  + 


1 


/^'(A-:')^  f)=k)z' 


which,  upon  using  (16),  (18).  and  the  equation  dual  to  (17).  simplifies  to 


GUky=\z']  =  r,^^i:.{kyz\z') 


53 1 


(51) 


(55) 


(57) 


The  elements  of  ,  as  well  as  the  corresponding  CP,  have  thus  been  determined. 
1  heir  space  domain  counterparts  are  found  in  the  same  manner  a.s  in  Formulation  .\. 
and  take  the  form 

u  r... . Wt 


gr{p-,z\z')  =  -i,So\^vt(K-z\z')  +  -^v[{k 


16 


17 


+  -{[XX- yy)  cos  2^  +  {xy  +  yx)  .sin  2^]  Sj 


/‘.Jr: 


i;''(Av;  cir'l 


-  {zxcos^  +  zy  sin  >p)jk-or)o^t  (  I  - 

+  ii7;o  ^-^Soirjky.zl:']} 

f  , 


k;  ik:'r 
rjkyzlz')] 
u"{k:’r  J 


k- 


! 


and 


cr{p:  = 


)  = 


iZiikiiJill 

v''{k'J)^  j 


(lilt) 


The  remark  following  (49)  also  applies  in  the  present  case.  Observe  that  wluni  the 
source  and  field  points  are  within  the  same  material  layer.  G  ’  becomes  a  diagonal 
dyadic  and  the  zx  and  zij  components  of  vanish,  thus  resulting  in  a  significant 
simplification  of  this  MPIE  formulation.  Note,  however,  that  G'^  in  this  formulation  is 
a  discontinuous  function  of  z'  across  any  interface  between  dissimilar  material  layers. 


4.3  Formulation  C 

To  arrive  at  P'ormulation  C,  v\e  set 


where  we  have  used  (39).  It  then  follows  from  (33)  that 


Upon  using  the  above  in  (35),  we  obtain 

noki,  J  d  . 


G 


A 

XLZ 


kot' 


v:r  -  ^  ^  ( vt  -  V" 


k,  dz'  ^  ‘ 


which,  in  view  of  (16)  and  (18),  may  be  reduced  to 

koVop't 


GtikyzW) 


k„ 


V^{kyz\z')-V:{ky=\z') 


(60) 


(61) 


(62) 
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Similarly,  upon  using  (61)  in  {'My),  wp  obtain 


riokp  je  ^ 


^'Of:  *  ^'p  <  ) 

which,  in  view  of  (16).  may  he  reduced  to 


!  ti  l  1 


k-„ 


I^{kyz\z')-  i:{k,:zW} 


i  (>■') ! 


Finally,  from  (61)  and  (d7).  we  obtain 


~  / A’  ,  ^0  c/ ,  I  _L  f  I [ 


(66) 


which,  upon  using  (16),  (18),  and  the  equation  dual  to  (17).  may  be  e.vpressed  as 


GfM:  .-I--')  =  ( ( ^ )  /;«>; + 


I 


+ 


/<(Cj  //'jcl  \  k}, 


(k-of 

1.  .  -  :  ~e.  1 

UJ, 

The  elements  of  as  well  as  the  corresponding  G"^.  have  thus  been  determitied. 
Their  space  domain  counterparts  are  found  in  the  same  manner  as  in  Formulation  .V. 
and  take  the  form 

....  ,  ^  f/yC-l-'O-WJ-,:--!--')! 

-  (zxcosi^-j-zi/sin<f)jkoT/oMtO}j^ - p - 1 


-  (xzcos</>  +  pz  sm  <^)  jkoTiofi'^S]' 


v:^(k,-zW)~v;ik,:=\z'} 


+  ZZ  Tj^fitpi't  So 


1  1 

] - - 

PtC,  ;r(t; 


(kof 

r'7  i>  ‘  -''i  > 

1  h  I 

\  / 

J 

and 


r;^(^;c|/)  =  5o' 
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The  remark  following  (49)  also  applies  in  the  present  rasi’.  Ohsf-rve  that  ( r'  in  tins 
formulation  is  a  continuous  function  of  both  r  and  z'  across  any  interface  Ijet  ween 
dissimilar  material  tavers. 


5  Summary  and  Conclusions 

We  present  a  systematic  derivation  of  three  distinct  mixed-potential  integral  erjuation 
(MPIE)  formulations,  referred  to  as  Formulations  A.  B,  and  C,  for  arbitrarily  sha[)('d 
non-planar  conducting  objects  embedded  in  a  medium  comprising  any  number  of 
planar  material  layers  of  infinite  lateral  extent,  which  may  be  uniaxially  anisotropic. 
The  development  is  carried  out  in  the  spectral  domain,  which  effectively  reduces  the 
original  vector  problem  to  a  much  simpler,  scalar  transmission  line  problem.  The 
spectral  domain  vector  and  scalar  potential  kernels  of  the  three  formulations  have 
explicitly  been  expressed  in  terms  of  the  voltage  and  current  Green's  functions  of  the 
transmission  line  network  analog  of  the  layered  medium.  It  is  found  that  the  three 
MPIE  formulations  are  of  similar  complexity  for  arbitrarily  shaped  structures,  but 
for  strictly  planar  microstrip  geometries,  Formulation  C  is  the  simplest  and  thus  the 
most  efficient  one.  This  formulation  is  also  recommended  when  the  conductor  pene¬ 
trates  an  interface  between  dissimilar  material  layers,  in  which  case  Formulations 
and  B  are  not  applicable,  unless  the  MPIE  (10)  is  augmented  by  additional  contour 
integral  terms,  as  discussed  by  Michalski  and  Zheng  [9],  [10].  However,  when  the  con¬ 
ducting  structure  is  confined  to  a  single  material  layer.  Formulations  .\  and  B  become 
attractive,  because  in  that  case  two  off-diagonal  components  of  their  vector  potential 
kernels  vanish,  which  results  in  considerable  savings  in  the  computational  effort.  The 
reader  is  referred  to  [9],  [10]  for  a  more  detailed  discussion  of  the  properties  of  the 
three  formulations. 

The  advantage  of  the  MPIE  approach  presented  here  is  that  it  involves  kernels 
which  are  less  singular  and  comprise  faster  converging  spectral  integrals  than  the  ker¬ 
nels  encountered  in  the  standard  electric  field  integral  equation  (EFIE)  formulation. 
Although  th'^  MPIE  is  more  complex  than  the  EFIE.  it  is  amenable  to  well-established 
numerical  solution  procedures  originally  developed  for  scatterers  of  arbitrary  shape 
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in  homogeneous  media  [1],  [5].  The  MPIE  approaeli  may  not  be  tiie  most  ejficient 
one  to  use  in  the  case  of  planar  structures  witli  simple,  regular  shapes,  but  it  (dfers 
the  flexibility  and  accuracy  that  make  it  well  suited  for  the  analysis  of  non[tianar 
and/or  irregularly  shaped  microstrip  transmi-ssion  lines  and  antennas  frequently  en¬ 
countered  in  modern  microwave  and  mi''imeter-wave  systems.  The  MPIK  is  also  a 
promising  candidate  for  use  in  hybrid  formulations,  which  combine  difh'rent iai  and 
integral  equation  methods  [31],  [32], 


A  Transmission  Line  Green’s  Functions 


Consider  a  transmission  line  section  corresponding  to  the  nth  material  layer  of  Fig.  2. 
characterized  by  the  parameters  (1).  Let  this  line  section  comprise  a  unit-^strength 
current  source  i  at  z\  as  illustrated  in  Fig.  .5.  It  then  follows  that  the  voltage  I/’  and 
current  If  at  any  point  s  within  this  section  obey  ihe  equations  (16)-(17).  in  which 
kP  =  kP^  and  Yp  —  V]f.  From  these  equations,  we  find  that  V]'’  may  be  expressed  in 


Figure  5:  Transmission  line  section  comprising  a  unit-strength  current  source. 
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where  =:max(c.c'),  Cq  is  an  arbitrary  refereiuf  jHjiiit  within  'in- 

ttth  section  (which  is  usually  set  to  c„  or  whichevt'r  is  more  ((ui'.emcn!  ).  am! 

r„(co)  and  r,i(cQ)  are  the  reflectiori  coelfieients  ’looking  to  tlse  left’  am!  ’lookitig  to 
the  right,"  respectively,  at  ro-  fiy  means  of  the  translation  formula 

f„{  =  )  .7!) 

where  the  upper  and  lower  signs  correspond  to  the  right  and  left  arrow--,  resfn'e 
tively,  these  reflection  coefficients  can  be  e.xpressed  in  terms  of  the  res[)ecli\e  termi¬ 
nal  reflection  coefficients  and  Fn,  which  are  related  to  the  corresponding  terminal 
admittances  Yn  and  (see  Fig.  5)  as 


Fn+  Yn  1+  F,i 


For  a  transmission  line  section  of  a  finite  length  a  particularly  convenient  form  t)f 
(70)  is 

e-jkzn\z-z'\  [i^.  r;  £-^'■^''^'■•<-<-"'•>1  [l+  "Fn 

V((c|2')  = - - 5==’=^= - ^  (7d) 

“^n  I-  r„  F„ 


which  is  obtained  by  letting  Zo  —  Zn  and  using  (71). 

The  remaining  transmission  line  Green's  functions  can  be  derived  from  (7:F). 
Hence,  /„(z|;'),  which  is  dual  to  V'((.;|c:'),  is  obtained  when  is  replaced  by  Z,., 
in  (73)  (which  causes  the  reflection  coefficients  to  change  signs),  l,{z\z')  follows  from 
(16)  and  (73),  and  Vv{z\z')  may  then  be  obtained  from  the  reciprocity  relation  given 
in  (18). 

The  voltage  and  current  on  the  nth  transmission  line  .section  that  is  source-free 
satisfy  (16)-(17)  with  the  delta  function  absent.  From  these  equations,  the  voltage 
at  any  point  within  the  line  section  may  be  e.xpressed  in  terms  of  the  voltage  \  '(To) 


across  one  of  hs  terminal  pairs.  The  result  is 


=  ^ -  14-  i;  ,71; 

1  -f  r„  f  J 

where  Zq  =  :„  or  cq  — (.see  Fig.  oj.  clepeiuling  on  vv!u>th<'r  the  source  is  hjcated  oj 
the  left  or  to  the  right,  respectively,  of  the  line  sec  tion,  and  where  the  upper  ilowcTi 
arrow  corresponds  to  c  >  Cq  (r  <  ^o).  We  have  omitted  the  siihsrripl  of  in  ^71  s, 
because  the  latter  applies  irrespective  of  the  nature  of  the  source,  and  drop[)e(i  frou! 
its  argument,  because  (74)  only  implicitly  depettds  on  the  source  location,  whic  h  i' 
outside  the  nth  line  section.  The  current  [{z)  corresponding  to  (71)  may  be  obtaineci 
by  substituting  the  latter  into  (16). 

Finally,  we  note  that  (71)-(72)  are  easily  implemented  in  a  recursive'  computc'r 
routine  to  determine  the  leftward-  and  rightward-looking  reflection  coelficients  needed 
in  (73)  and  (74).  The  computations  proceed  from  the  outward-looking  refh'clion  i,o- 
efficients  in  the  top  and  bottom  transmission  line  sections  (see  Fig.  4).  where  they 
vanish  for  unshielded  structures  or  othervvise  are  easily  determined,  provided  the  sur¬ 
face  admittances  Vs.i  and  Vs,,v-i-i  are  known.  Observe  that  the  e.xponentiai  functions 
encountered  in  (71),  (73),  and  (74)  have  nonincreasing  magnitudes,  so  there  is  no 
danger  of  overflow.  When  the  nth  transmission  line  section  extends  to  positive  ( nec- 
ative)  infinity  along  the  2  axis,  the  reflection  coefficient  r,i  f  l'„)  is  set  to  zero  in  (73! 
and  (74). 
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Abstract 

.4  straightforward  derivation  is  presented  of  three  distinct  mixed-potential  integral 
equations  (MPIEs)  that  are  well-suited  for  the  analysis  of  arbitrarily  shaped,  planar 
or  non-planar  microstrip  structures  in  multi-layered,  planar,  uniaxially  anisotropic 
media  of  infinite  lateral  extent. 

1  Introduction 

The  integral  equation-method  of  moments  (lE-MOM)  approach,  implemented  either 
in  the  spectrum  or  space  domain,  continues  to  be  the  method  of  choice  for  the  analysis 
of  open  microstrip  structures  (see  [1]  and  [2]  for  recent  examples  of  the  application  of 
these  techniques).  The  spectral  and  space  domain  IE  formulations,  both  being  ba.sed 
on  the  rigorous  Green’s  function  of  the  background  medium,  are  formally  identical, 
and  only  differ  in  the  order  in  which  the  integrals  that  arise  are  performed.  When 
applied  to  planar  geometries,  the  spectra!  domain  approach  (SD.4)  results  in  very 
simple  expressions  for  the  MOM  matrix  elements,  which,  however,  comprise  slowly 
convergent  double  spectral  integrals.  Consequently,  to  make  this  technique  prar  tical. 
a  considerable  effort  is  required  to  accelerate  these  integrals,  and  the  initial  simplicity 
of  the  formulation  is  lost  in  the  proce.ss  [3],  [4],  [5].  .Although  the  SDA  can  in  princi¬ 
ple  be  extended  to  non-planar  conductors,  the  space  domain  integrals  that  must  be 
analytically  evaluated  involve  in  that  case  the  layered  medium  Green’s  function,  and 
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the  cieveiopment  is  tedious  and  results  in  several  s[)ertrai  inte^rais  lur  each  luarnx 
element,  thus  further  eoruplicating  the  acceleration  procedure. 

The  space  domain  IE  is  usually  first  transformed  into  a  mixed-potential  form,  in 
which  the  electric  field  is  expressed  in  terms  of  the  vector  and  scalar  |)ot('nt  ials,  1..  tor.- 
it  is  subjected  to  the  MOM  [6].  Tiie  resulting  IE.  which  is  often  referre.l  to  in  ilu-  liter 
ature  as  the  mixed-potential  integral  equation  (MPIE)  [7|.  [Sj.  has  a  mor.'  coiiifiii.  ate.i 
form  than  its  spectral  domain  counterpart,  hut  has  many  advantages,  making  it  \er’. 
attractive,  particularly  for  non-planar  and  arbitrarily  shaped  microstrip  structures. 
First,  the  MPIE  potential  kernels  are  less  singular  wtu'ii  the  source  and  .)bsfr\at  ion 
points  coincide  oa  the  surface  of  the  object  than  the  dyadic  kernel  of  the  siandar.l 
electric  field  integral  equation  (EFIE)  (in  the  SD.-\  implementation  of  tin'  EFIE.  tlie 
severe  kernel  singularity  of  the  latter  is  manifested  by  a  divergent  behavior  of  the  dou¬ 
ble  spectral  integrals).  Second,  the  .MPIE  is  amenable  to  well-established  numerical 
solution  procedures  [6],  [9],  [10],  [11]  and — in  contrast  to  the  SD.-\--it  only  requires 
single  spectral  integrals  of  the  Sommerfeld  type  [12.  p.  210|.  which  are  independent 
of  the  testing  and  basis  functions  employed  in  the  MOM  procedure.  'I'hird.  these 
Sommerfeld  integrals  are  amenable  to  various  approximation  methods,  such  as  the 
interpolation  and  table  look-up  scheme  [13]  (see  also  [14].  [15]).  or  the  comph'x-image 
representation  [16],  [17|.  which  drastically  reduce  the  computational  effort. 

In  a  layered  medium,  different  .MPIE  formulations  are  possible,  as  a  result  of  the 
non-uniqueness  of  the  potentials.  There  are  at  least  three  possible  useful  choices  of  the 
magnetic  vector  potential  (or  electric  Hertz  vector)  associated  with  a  horizontal  time- 
harmonic  dipole  [18],  which  also  lead  to  different  scalar  potentials.  Of  the  three  vector 
potential  choices,  the  classical  Sommerfeld's  form  [12.  p.  258]  (termed  ■fraditionar 
in  [19])  has  been  and  remains  the  most  popular  and  useful.  .Another  form,  termed 
‘alternative,’  has  been  shown  to  have  some  advantages  over  the  traditional  one  in 
the  MPIE  formulation  [19].  As  was  first  pointed  out  by  Mosig  and  (lardiol  [13], 
the  scalar  potential  kernel  of  the  MPIE  may  be  interpreted  as  the  potential  of  a 
single  point  charge  associated  with  a  time-harmonic  dipole.  In  a  layered  medium, 
the  potentials  associated  with  the  horizontal  and  vertical  dipoles  are  different  [19], 
which  poses  a  difficulty  in  the  case  of  arbitrarily  shaped,  non-planar  objects,  because 


only  one  scalar  potential  appears  in  the  stamiard  form  of  tiie  MPIK.  A  lueUiod  to 
overcome  this  difficulty  was  first  proposed  by  Michalski  [71.  flie  two  fortiis  of  the 
vector  potential  (traditional  and  alternative)  and  two  scalar  potentials  fhori/ontal 
and  vertical)  lead  to  three  distinct  Mf^IE  formulations,  which  were  first  dev('!oped 
and  explored  in  detail  by  Michalski  and  Zheng  [20].  [21]. 

The  MPIE  approach  was  first  applied  to  planar  microstrip  structures  with  isot  topic 
substrates  by  Mosig  and  Cardiol  [13].  Their  formulation  was  later  adopt('d  lyv  others 
[22],  [23],  [24],  [25],  [26],  [27]  and  was  extended  to  multi-layered  media  [1  Ij.  [15]  and 
to  electrically  uniaxial  media  [28].  One  of  the  .MPIE  forms  developed  in  [20],  referred 
to  as  Formulation  C  (which  corresponds  to  the  ‘traditional  formulation’  of  [10] h  lias 
recently  been  employed  to  analyze  coax-fed  microstrip  patch  antennas  of  arbitrary 
shape  [29],  [30],  [31].  as  well  as  microstrip  transmission  lines  of  arbitrary  cross-section 
[32],  [33].  This  formulation,  which  is  based  on  the  .Sommerfeld’s  form  of  the  vector 
potential  and  the  scalar  potential  associated  with  a  horizontal  dipole,  reduces  for 
planar  microstrip  structures  with  isotropic  substrates  to  that  originally  rleveloped  by 
.Mosig  and  Cardiol  [13]  (see  also  [bo],  [14]).  A  modification  of  Formulation  C  was 
adopted  by  .Montgomery  and  Wilton  [.34]  (see  also  [35])  for  the  analysis  of  arbitrary 
conducting  periodic  structures  embedded  in  layered  media.  For  mirrostrip  geometries 
that  only  comprise  planar  and  vertical  components,  such  as  coax  probe-fed  microstrip 
patch  antennas,  hybrid  MPIE  formulations  have  been  developed,  which  utilize  two 
scalar  potential  kernels  (associated  with  the  horizontal  and  vertical  dipoles)  and.  as  a 
result,  require  a  point  charge  at  the  probe-to-patch  junction  [36],  [37],  [38].  [39],  [lOj. 
A  still  different  MPIE  formulation  has  recently  been  introduced  by  Mandenbosch  and 
Van  de  Capelle  [41], 

The  purpose  of  this  paper  is  to  present  a  new.  straightforward  development  of 
three  distinct  MPIE  formulations  (corresponding  to  Formtdations  .A,  B.  and  (‘  of  [20]  ). 
that  are  suitable  for  arbitrarily  shaped,  non-planar.  microstrip  structures  emb('dd(>d 
multi-layered,  planar,  uniaxial  media  of  infinite  lateral  extent. 


2  Preliminaries 


Consider  a  medium  comprising  .V  planar,  iiniaxialiy  aiiisot  ropi<-  Iay(‘r>-  of  infiniP'  lat 
eral  t'Xtent.  as  illustrated  in  Fig.  1.  The  nth  layer  is  characterized  bv  th(‘  perrnittivitv 

(^) 

A 


Figure  1;  .Medium  comprising  A'  planar,  uniaxially  anisotropic  layers, 
and  permeability  dyadics.  £„  and  /Lt„.  of  the  form 

~  ^  ^  (-zn  ^  ~  Xf  f^tn  //;n  (  1  i 

where  C(n  ifi-tn)  and  (fizn)  denote,  respectively,  the  transverse  and  longitudinal 
dielectric  (magnetic)  constants  relative  to  free  space,  and  X<  is  the  unit  dyadic  trans¬ 
verse  to  z.  Observe  that  we  distinguish  dyadics  by  double  underlines  and  unit  vectors 
by  carets.  The  free  space  permeability  and  permittivity  will  be  denoted  by  /(o  and  fo- 
respectively.  For  each  layer,  we  also  introduce  its  electric  and  magnetic  anisotropy 
ratios,  and  respectively,  given  as 


The  top  and  bottom  layers  may  be  of  infinite  extent  alont^  the  :  axi^  dp  a>  lihi^tra 
ted  in  Fig.  1  —may  be  shielded  by  ground  planes  liaving  spr-citied  suriact'  admittan¬ 
ces  Ks-.  v+i  and  Ys.i,  respectively,  where  the  latter  are  infinite  for  perfect  e!e(  U!< 
conductors  (PfXf). 

Let  an  arbitrarily  shaped  PEC  object,  whose  snrfac*'  may  be  open  or  closed,  be 
embedded  in  the  layered  rnedinmof  Fig.  1.  Furthermore.  let  E’  denot<'  the  impre'^-'i-d 
(primary)  time-harmonic  (the  time  convention  is  adopted  herei  tdectric  H<'ld  that 
exists  in  the  layered  medium  in  the  absence  of  the  object.  With  the  olpjecr  present, 
this  field  excites  on  5  the  surface  current  Js.  which  is  the  source  of  the  scattered 
(secondary)  electric  field.  The  latter  may  be  expressed  as 

E{r)  ^ 

where  Qf'  is  the  dyadic  Green's  function  of  the  layered  medium,  r  and  p  denote, 
respectively,  the  position  vector  and  its  projection  on  the  xy  plane,  and  where  primes 
are  used  to  distinguish  source  coordinates.  For  later  convenience,  we  have  introduced 
in  (3)  the  shorthand  notation  (;)  for  an  integral  of  a  product  of  two  functions  sep¬ 
arated  by  the  comma.  The  dot  over  the  comma  signifies  a  ‘dot  product'  of  vector 
arguments  and  the  prime  over  )  indicates  that  the  integration  is  over  the  primed 
(source)  coordinates,  while  the  subscript  5  designates  the  domain  of  integration. 

The  dyadic  Green's  function  in  (3)  may  be  expressed  as  the  inverse  X’ector  Fourier 
Transform  (cf.  [42])  of  its  spectral  domain  counterpart.  of  the  form 

.-I--')  =  -'I--')  •  g(e} 

-|-rX5  -k-oo 

-oo  -oo 

Here.  kj.  and  ky  are  the  spectral  domain  counterparts  of  x  and  //.  respectively,  and 
kp  =  xki+yky  is  the  radial  wave  vector.  Based  on  kp.  we  introduce  the  rotated  spectral 
domain  coordinate  system  illustrated  in  Fig.  2,  where  the  unit  vectors  u  and  v  are 
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(ky) 


given  as  (cf.  [43],  [44],  [45],  [46]) 

*  A  2  X  ASp  I  — 

“=77’  ^  =  ~T~'  =  + 

In  (4),  the  dyadic  operator 


Q(a  =  iux  +  vy)cos^  4- (uy  —  rx)  sin  ^  +  ii  (6) 

and  its  transpose,  Q  (^),  are  employed  to  directly  transform  Q_  from  the  rotated 
spectrum  space  to  the  configuration  space.  Upon  projecting  Fourier-transformed 
Maxwell’s  equations  on  the  rotated  coordinates  of  Fig.  2,  we  find  that  these  equations 
split  into  two  independent  sets  comprising,  respectively,  fields  that  are  transverse- 
magnetic  (TM)  and  transverse-electric  (TE)  to  c.  Furthermore,  it  is  found  that 
these  partial  fields  may  efficiently  be  expressed  in  terms  of  the  voltages  and  currents 
on  two  transmission  lines  along  the  c  axis,  with  the  propagation  wavenumbers  and 
characteristic  admittances  given  as 
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and 
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respectively,  where  k(^  =  u-’^/ioCo  is  the  wavenumber  and  ;?o  =  y'/'oAo  the  intrinsic 
impedance  of  free  space.  Observe  that  at  this  stage  of  the  development,  the  perme¬ 
ability  and  permittivity  (and  thus  also  the  anisotropy  ratios)  of  the  medium  may  vary 
arbitrarily  with  c  and  are  not  necessarily  restricted  to  be  piecewise-constant ,  as  is 
the  case  in  Fig.  1.  In  the  above  we  have  adopted  ♦ '  e  convention  that  the  quantities 
associated  with  the  TM  and  TE  transmission  line,  .re  distinguished  by  superscripts 
e  and  /i,  respectively.  For  the  following  development,  it  will  be  convenient  to  also 
introduce  transmission  line  Green's  functions,  which  represent  the  voltage  or  current 
excited  by  a  unit-strength  voltage  or  current  point  source.  Hence,  let  lA(cjc')  and 
/f(2|2'),  where  the  superscript  p  stands  for  e  or  h,  denote  the  voltage  and  current, 
respectively,  at  a  point  2  due  to  a  1  A  current  source  located  at  z'  on  the  correspond¬ 
ing  transmission  line  (cf.  [47],  [46]).  These  Green’s  functions  are  governed  by  the 
coupled  equations 


^VnAz')  =  -}k’Z>  (9) 


where  6{z—z')  is  the  Dirac  delta.  Similarly,  let  Vqf(2|2')  and  I^.{z\z')  denote  the  voltage 
and  current,  respectively,  at  2  due  to  a  1  V  voltage  source  at  z'.  It  then  follows  that 
these  Green's  functions  satisfy  equations  dual  to  (9)- (10).  which  are  obtained  from 
the  latter  by  making  the  substitutions  —>■  /,G  If  — *•  V[f,  and 

From  their  governing  equations,  it  can  also  be  shown  that  the  four  transmission  line 
Green’s  functions  possess  the  symmetry  and  reciprocity  properties  [48,  p.  194] 


V7(2|2')  =  Vf’iz'lz),  If(z\z')  =  /n~'i^-),  V^{z\z')  =  -/f(2'i2) 
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In  the  rotated  coordinate  system  of  Fig.  2,  the  spectral  domain  dyadic  in  (4) 
may  be  expressed  in  terms  of  the  transmission  line  Green's  functions  as  (see  also 


# 


34 


Kastner  et  al.  [47]) 


c|c')  =  till  V7(^;  z\:')  +  iv \]'^{ky.  rjc') 

-  .  ^/O^p  ,p.,  I  /V  .  .  ’I()kp  ,  -p  , 

-  zu  - - I-{ky  rj'  )  -  uz  j-j  \  „  (A-,,. 


k'ijt. 


i-  2Z 


'<iy^  KiK-A--’) 

.<■0/ 


A'of  A 

'/u 

jk\j(k 


H: 


(12) 


where  8  is  the  Dirac  delta  function,  and  where  the  primed  and  iinprinied  media  pa¬ 
rameters  are  evaluated  at  z'  and  respectively — a  convention  that  will  also  he  in 
effect  throughout  the  remainder  of  this  paper,  unless  the  layer  index  is  explicitly 
stated.  Observe  that  in  the  above  we  explicitly  indicate  the  dependence  of  the  trans¬ 
mission  line  Green’s  functions  on  the  transverse  spectral  wavenumber  kp.  Since  no 
assumptions  have  so  far  been  made  regarding  the  dependence  of  the  media  parame¬ 
ters  on  the  spectral  dyadic  Green's  function  given  in  (12)  is  valid,  in  particular,  for 
the  layered  medium  of  Fig.  1,  where  the  permeability  and  permittivity  are  piecewise- 
constant.  The  voltage  and  current  Green’s  functions  of  its  transmission  line  network 
analog  are  easily  found  for  any  number  of  layers,  as  discussed  in  the  appendix. 

VVe  note  that  the  transmission  line  Green's  functions  appearing  in  (12)  depemd  on 
the  spectral  variables  ky  and  ky  exclusively  through  kp.  As  a  consequence,  when  the 
double  spectral  integration  in  (4)  is  performed  in  the  polar  system  (Ap.^f)  shown  in 
Fig.  2.  the  if-integrals  can  be  evaluated  in  closed  form.  .More  specifically,  we  find  that 


where  we  have  introduced  the  notation 

SnUikp)}  =  rdkpfikp)k';+\h(kpg) 
Ztt  jo 

In  the  above, 

i)  =  arctan  + 


8 


35 

and  Jn  is  the  Bessel  function  of  order  ri.  [n  this  fashion,  the  space  <iuinain  dvadic  Q!' 
may  be  expressed  in  terms  of  the  Somrnerfeld-type  integrals  ( 1  1).  Its  (e\pli<  it  form  is 
easily  found  with  the  help  of  (4)  and  (I:}),  hut  will  not  be  required  in  the  iVdhnving 
development,  and  is  not  listed  here  in  the  interest  of  brevity. 

With  determined,  the  electric  field  scattered  by  tin'  conduct  ins  oliject  iiuiv 
be  found  from  (3).  provided  the  surface  current  J,-  is  given.  Wlnni  the  latter  is  not 
known  a  priori  as  is  usually  the  case,  we  may  use  (3)  to  formulate  an  IE  for  J< 
enforcing  the  condition  that  the  total  tangential  electric  field  must  vanish  at  .\s  a 
result,  we  obtain  the  EFIR 

-(Ak:{gf:Js%X  =  {A,:E%  (16) 

where  {A,}  is  a  suitably  chosen  complete  set  of  linearly  independent  vector  weight 
(testing)  functions  defined  over  and  tangential  to  .S'.  This  EFIE.  however,  is  not 
suitable  for  a  direct  application  of  the  MOM  [6]  because  the  dyadic  kernel  exhibits 
a  severe  singularity  when  the  source  and  observation  points  coincide  on  and  the 
Sommerfeld-type  integrals  that  arise  are  slowly  convergent. 

When  the  spectral  representation  of  given  in  (4)  is  substituted  into  ( 16)  and 
the  order  of  the  spectral  and  space  integrals  reversed,  the  latter  may  be  evaluated 
analytically,  provided  the  testing  functions  and  the  basis  functions  used  to 

represent  Js  are  Fourier-transformable  in  closed  form.  This  procedure  leads  to  the 
previously  mentioned  SDA,  which  has  been  popular  in  the  analysis  of  planar  rni- 
crostrip  structures.  Unfortunately,  the  double  spectral  integrals  associated  with  this 
approach  suffer  from  an  extremely  slow'  convergence  and — except  for  simple,  regular 
geometries,  for  w'hich  the  MOM  matrix  has  a  small  size — the  method  is  inefficient. 

The  difficulties  encountered  when  the  MOM  is  applied  to  the  EFIE  may  be  alle¬ 
viated  if  the  secondary  electric  field  (3)  is  expressed  as 


-E{t)^  Js%  +  V(G*^  V'-  Js% 


in  W'hich  V'  operates  on  source  coordinates,  and  where  and  are  referred 
to,  respectively,  as  the  vector  and  scalar  potential  kernels.  Because  both  potentials 
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^6 


appear  in  (17).  the  latter  is  termeil  the  niixetl-potent ial  representation  ih  the  i-iec  i-ir 
field  [7j.  [49],  I  poii  using  (17).  instead  of  (3).  in  the  procedure  tfiat  led  ttj  tiie  Kf-'Il-. 
(16),  and  applying  a  Gauss'  theorem  [50.  p.  503|.  we  arrive  at  !h<>  .MPll'. 


-  (v-a,.  (US  V'.  ^  {,u:  es. 


1  Vi 


provided  G'^  is  a  continuous  function  of  c  across  any  interfact'  itetween  divsimifa.r 
material  layers  the  object  pt'iietrates.  Gornparing  the  EFIF.  (16)  with  th(>  MIMF  i  1  sj. 
we  note  that  in  the  latter  the  differential  operator  W'  has  in  efft'cl  been  exiracteti 
from  and  transferred  onto  yl*;  and  [.35].  As  a  consequence,  the  kernels  G  '  a 


and 


appearing  in  the  .MPIE  (18)  are  less  singular  (or  ‘less  discontinuous'  [50.  p.  1!)5A 
than  and  the  resulting  Sommerfeld-type  integrals  converge  more  ra[)idly.  Clearly, 
for  this  equation  to  be  applicable,  the  testing  functions  and  the  liasis  functions  used 
to  represent  Js  must  be  at  least  once  differentiable  — a  requirement  tnat  can  <'asily  be 
accommodated  [10].  Finally,  we  note  that  (18)  may  readily  be  extended  to  irnpt'rhx  t 
conductors  via  the  surface  impedance  concept  [8],  [38]. 

In  the  next  section,  we  employ  the  spectral  dyadic  given  in  (12)  as  the  point 
of  departure  in  the  development  of  the  vector  and  scalar  potential  kernels.  G  '  and. 

respectively,  associated  with  three  distinct  MPIE  formulations,  for  rnulti-laycrerl, 
planar,  uniaxial  media  of  infinite  lateral  extent. 


3  MPIE  Development 

The  development  is  carried  out  in  the  spectral  domain  and  it  proceerls  from  the 
spectral  electric  dyadic  Green's  function  given  in  (12).  First,  we  add  to  and 
subtract  from  (12)  the  term  tiii  V]''(Ay:  cjc')  and  recognize  that  uu  +  vv  Next, 

we  observe  that  Vt=:  -jk^u,  where  V,  is  the  spectral  domain  counterpart  of  -  the 
trarnsverse  (to  z)  part  of  the  operator  nabla.  and  that  v/  =  jk^u.  .Xs  a  results  of  thest' 
steps,  we  obtain 

-  g^(fcp;  -  Lt  -  v,v; 


V/^{k,;z\z')-V:{k,:z\z') 
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Next,  \vt’  use  tl',','  [)r()j)«'rt it‘s  of  the  !r.U!S!ii!S“-ii>ii  iiiu-  (Irci-it'^  iniH  iuju-'  in  (< 
with  the  !(lent!t\' 


!'(< 


n  »!i 


vv'  =  v,v/  -  ^  z  —  v'  - 

(f:'  (h  ifzih’ 


e_'U^ 


to  express  in  tin'  form 


=  GUk, -  VV  (r  i  kz 


i  s  i 


which  in  tfie  space  domain  becomes 

-  gf-ip-p':  ;j;')  =  Q'ip-p'- -  W'C-ip-p'-.  i-n) 

I'pon  using  the  abo\'e  in  {'■])  and  apitealing  to  the  diversience  th('orem  !■)().  p. 
we  arrive  at  the  mixed-potential  representation  o['  E  given  in  (  ITi  provided  ( !'■'  is  a 
continuous  function  of  z'  across  any  interface  Ix'tween  dissimilar  material  layers  the 
object  penetrates,  which  then  leads  to  the  desirefl  MPIK  I  18). 

The  decomposition  of  according  to  (21)  is  not  unique  and  leads  to  dilferent 
MPIE  forms,  depending  on  the  choice  oi  S’*  '"ctid  Below,  we  develoj)  three 

distinct  .\IPIE  formulations  that  have  been  found  useful  in  MOM  analyses  of  the 
radiation,  scattering,  and  guidance  of  electromagnetic  fields  by  arbitrarily  "shaped 
conducting  structures  in  layered  media,  in  the  interest  of  brevity,  simie  details  of 
the  derivations — which  repeatedly  invoke  the  equations  fb)  (10)  and  their  duals,  as 
well  as  the  reciprocity  and  symmetry  properties  (11)  of  the  transmission  line  (Ireen  s 
fu  net  ions — are  orn  i  t  ted. 
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Formulation  A 

By  invokiug  the  properties  oi  the  transniissioti  line  (ireeii'*'  fuiietions.  we  iiiav  expr<*s^ 
(19)  as 


=  L  y/'{k^=\=')  -  v.tt 


\';‘{L:z\z')  -  v:\k- z^z'i' 


+  c|c')  +  i£  /^(A- :  ri j 

(  .0  ->  .  .  d  d  \  \;^(Av;r|y) 


3— V  +  ££ - '  ‘  : 

dz  ^^  dzdz') 


provided  the  field  point  is  within  a  homogeneous  material  layer.  The  form  of  the 
above  suggests  that  vv  '  may  choose 

as  the  spectral  domain  scalar  potential  kernel  (this  choice  corresponds  to  Formula¬ 
tion  A  of  [20j).  Next,  we  make  use  of  the  identity  (20)  and  combine  terms  in  (29  i  to 
arrive  at  (21).  in  which 


a>,<,  vak.:]:'} 

‘  A-2  /  ••£  A'-^ 


The  space  domain  counterparts  of  (25)  and  (24)  are  found  as 

.-|r')  =  jX,  5„|  VyAv;  --|=')  +  ^  Kd.:  -'I--')} 

+  h(xx-i/y)cos2iJ+(ii/+yi)sin'2t?]52|  — I  _  -~_i  '  — - 

2  I  ^  ( k- )  V 

_  .a/'.  \  „  f  V"(  1 '' )  1 


-  (  x£  COS  d  +  y i  sin  d )  jkorfo/,  ^  ^  j  *^1 1 
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+  ’io~  -I-')} 


G'^(p-p'::\=')  =  So 


Vv(A^:~|c') 


respectively,  where  we  have  used  the  notation  introduced  in  ( i:5)  ( lo).  ()bserv<>  iliat 
when  the  source  and  held  points  are  within  the  same  material  layer,  the  .v:  and  >/: 
components  of  S"'  vanish,  thus  resulting  in  a  signihcant  simplification  of  this  MIMK 
formulation.  Note,  however,  that  CG  in  this  formulation  is  a  discontinuous  fimciion 
of  r  across  any  interface  between  dissimilar  material  layers. 

Formulation  B 

Upon  using  the  properties  of  the  transmission  line  Green's  functions,  we  may  also 
express  (19)  as 


-  ^-1-'')  =  L  z\z')  - 


-  iV,'  +  ii  ,5^  nay. 

.  .y  d\V,'(l.yz\-:) 


~  o  ..  a  a  \ 


provided  the  source  point  is  within  a  homogeneous  material  layer.  The  form  of  the 
above  suggests  that  we  may  choose 

as  the  spectral  domain  scalar  potential  kernel  (this  choice  corresponds  to  I'orrnula- 
tion  B  of  [20]).  Next,  we  make  use  of  the  identity  (20)  and  combine  terms  in  (2S)  to 
obtain  (21).  in  which 


a\KrAz')  =  L\\HK-z\z') 
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-  zVf  jk'oTjQ^t  1 


[ktf  kj 

^iW\l:ifc,-=\z') 


‘r;  \  ‘i  ''V'  -!-  '  ,  •  j/0  ,,  ,  ,  ,  ,, 

-  1 - —  +  ZZ  //,,  —  /,  I  A',, ;  c  1  c  ) 

IhtJ  u^'{k-r  — 

In  the  space  domain,  (30)  and  (29)  become 

a^p-p': + 


1 


+  ~[{xx  —  yy)  cos  2d  +  (ac^d-y®)  sin2dj S2< 


1^:') 


,2 


\  ^  1  mk}y::W) 


-  {zx  cos  d  +  zysmt})j koUo fit  \  I - ^  h^i  1  ,,2 

+  ii,3^5„{/:(*i;cl.-')) 

i:' 


t. 


and 


G'^{p-p'-.z\z')^So< 


vtiK^zW) 


(30) 


v^K^zW)  v;{k,-\z') 


(31 


(32 


t  J 

respectively.  Observe  that  when  the  source  and  field  points  are  within  the  same 
material  layer,  the  zx  and  zy  components  of  vanish,  thus  resulting  in  a  significant 
simplification  of  this  MPIE  formulation.  Note,  however,  that  in  this  formulation  is 
a  discontinuous  function  of  z'  across  anv  interface  between  dissimilar  material  lavers. 


Formulation  C 

The  form  of  (19)  suggests  that  we  may  also  choose 

v!'{K,z\z')-vy{k,-,. 


I~M 


G\K-z\z') 


¥ 


(33) 


as  the  spectral  scalar  potential  kernel  (this  choice,  which  was  previously  used  in  [28, 
for  a  planar  microstrip  geometry,  corresponds  to  Formulation  C  of  [20]).  Next,  we 
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make  use  of  the  identity  (20)  and  combine  terms  in  (10)  to  obtain  121  ).  in  u  liich 


-  iV(  jkorjoi.1, 

-1-  V,ijl’o//o/4 

[v-;"(A^ 

+  zzT]^fitn',  |( 

(rj^. 

A-r 


V-;f(A^:r|c')- V-(A>:c!-') 


1  I 

—  + 


i:{k,-.z\:']\  CM) 


The  space  domain  counterparts  of  (34)  and  (33)  are  found  as 

a'ip-p'-.  =  L  S„[V‘‘(ky. 


(ix  cos  d  +  iy  sin  d)  jkorjofitSi 


lSkp:=W)-i:(k,:z\:') 


...  ,  ,  .  .  .  V^ky,  .-|i')  -  V-(A>:  z\=') 

-  {xzcosx}  +  yzsiTn})jkQT)ofifS\< -  - 


1  1 

- 1 — — 

/be,-  /i<e'. 


— 

k. 


Kikp::W)\ 


and 


G'^ip~p';z\z'}^So^ 


V/^{k,:z\z')-Vnkp:z\z') 


135) 


(36) 


respectively.  Observe  that  in  this  formulation  is  a  continuous  function  of  both  r 
and  z'  across  anv  interface  between  dissimilar  material  lavers. 


Discussion 

Although  the  three  MPIE  formulations  are  of  similar  complexity  for  arbitrarily  shaped 
conductors,  Formulation  C  is  the  simplest  (and  thus  the  most  efficient)  one  for  strictly 
planar  microstrip  geometries.  This  formulation  is  also  recommended  when  the  con¬ 
ductor  penetrates  an  interface  between  dissimilar  material  layers,  in  which  ca.se  For¬ 
mulations  A  and  B  are  not  applicable,  unless  the  MPIE  (18)  is  augmented  by  addi- 
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^  tional  contour  integral  terms,  as  discussed  by  Michaiski  and  Zheng  iiOl.  lb»\ve\rT. 

wlien  the  microstrip  structure  is  confined  to  a  single  mati'rial  layer  (as  is  often  the 
case).  Formulations  A  and  B  become  attractive,  because  in  that  case  two  off-diageuia! 
components  of  their  vector  potential  kernels  vanish,  which  may  result  in  signilicant 
H  savings  in  the  computational  effort.  The  reader  is  referred  to  [201.  121]  for  a  mor<' 

detailed  discussion  of  the  properties  of  the  three  formulations. 


4  Conclusion 

.A  new.  straightforward  development  has  been  presented  of  three  distinct  mixed- 
potential  integral  equation  (.VIPIE)  formulations  for  non-planar.  arbitrarily  shaped 
microstrip  structures  embedded  in  planar,  laterally  open,  multi-layered,  uniaxial  me¬ 
dia.  These  MPIEs  may  be  used  in  conjunction  with  the  method  of  moments  to  ana¬ 
lyze  three-dimensional  microstrip  discontinuities,  such  as  transitions  between  printed 
traces  on  different  levels,  vias.  air-bridges,  bond-wires,  etc.,  as  well  as  microstrip  patch 
antennas  with  shorting  pins  and  coaxial  probe  feeds. 

A  Transmission  Line  Green  s  Functions 

The  transmission  line  analog  of  the  layered  medium  invoked  in  Section  2  comprises  two 
transmission  line  networks — having  identical  configurations,  but  different  propagation 
wavenumbers  and  characteristic  admittances — associated  with  the  TM  and  TE  partial 
fields.  In  each  of  these  two  networks,  a  homogeneous  material  layer  of  the  medium 
of  Fig.  1  is  represented  by  a  uniform  transmission  line  section.  In  the  rase  of  an 
unshielded  structure,  where  the  top  (or  bottom)  layer  is  of  infinite  extent  along  the 
2  axis,  the  corresponding  transmission  line  section  also  extends  to  infinity,  where  it 
is  terminated  into  a  matched  load. 

Consider  a  transmission  line  section  corresponding  to  the  nth  material  layer,  char¬ 
acterized  by  the  parameters  (1).  Let  this  line  section  comprise  a  unit-strength  current 
source  i  at  2',  as  illustrated  in  Fig.  3.  It  then  follows  that  the  voltage  V[''  and  current 
/f  at  any  point  2  within  this  section  are  governed  by  (y)-(lO).  in  which  AT  —  AT,, 
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Figure  3:  Transmission  line  section  comprising  a  unil-slrerigt h  current 


sourc('. 


and  y'P=  y'P,  where  the  latter  are  obtained  by  substituting  the  matenal  parameters 
of  the  layer  into  (7)-(S).  From  (9)-(10).  we  find  that  V,'’  may  be  expressed  in  the 
traveling- wave  form  (cf.  [48.  p.  213|) 


_ _ _  J: _ _ J_L _  } 

—  ^  n 


1 n( -o) In( ~o) 


(37) 


where  =  min(:r,c'),  r>  =  max(c.g')'  ~o  is  an  arbitrary  reference  point  within  tht' 
nth  section  (which  is  usually  set  to  or  whichever  is  more  convenient),  and 

r„(go)  and  r„(go)  are  the  reflection  coefficients  ’looking  to  the  left’  and  ’looking  to  the 
right,'  respectively,  at  ^o-  In  the  above  and  in  what  follows,  we  omit  the  superscript 
p  for  notational  simplicity.  By  means  of  the  translation  formula 


(38) 


w’here  the  upper  and  lower  signs  correspond  to  the  right  and  left  arrows,  respec 
tively,  these  reflection  coefficients  can  be  expressed  in  terms  of  the  respective  termi 
nal  reflection  coefficients  r„  and  F^i.  which  are  related  to  the  corresponding  terminal 
admittances  ^  and  (see  Fig.  3)  as 


Fn  +  y  n  I  +  I  ri 
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For  a  transmission  line  section  of  a  finite  lengtft  a  particularly  convenient  form  uf 
(37)  is 


=  - tt: - ^ - — 

iV.  |1_  ]'  r 


(  lOi 


which  is  obtained  lyv  letting  ro=r,i  and  using  (38). 

The  remaining  CJreen's  functions  can  easily  be  derived  from  (40).  Htuice.  I,  iTi:'). 
•  which  is  dual  to  l',(r|.;'),  is  obtained  when  }„  is  replaced  by  Zr,  in  (40)  (which  rau.ses 

the  reflection  coefficients  to  change  signs).  J,(c\z')  follows  from  (40)  and  (O).  and 
l'i,(c|c')  may  then  be  obtained  from  the  reciprocity  relation  given  in  (11). 

The  voltage  and  current  on  the  nth  transmission  line  section  that  is  source-free 
®  satisfy  (9)-(10)  with  the  delta  function  absent.  From  these  equations,  the  voltage  at 

any  point  c  within  the  line  section  may  be  expressed  in  terms  of  the  voltage  V  (co) 
across  one  of  its  terminal  pairs.  The  result  is 


V(‘)  =  Vizo) 


-'I 


1+  r,j  e-j2*zntin 


1+  t; 


(41) 


i 


w’here  ^o  =  -n  or  co  =  ;n-(-i  (see  Fig.  3).  depending  on  whether  the  source  is  located  to 
the  left  or  to  the  right,  respectively,  of  the  line  section,  and  where  the  upper  (lower) 
arrow  corresponds  to  s  >  Cq  (;  <  Cq)-  We  have  omitted  the  subscript  of  V  in  (41). 
because  the  latter  applies  irrespective  of  the  nature  of  the  source,  and  dropped  z'  from 
its  argument,  because  (41)  only  implicitly  depends  on  the  source  location,  which  is 
outside  the  nth  line  section.  The  current  /(c)  corresponding  to  (41)  may  be  obtained 
by  substituting  the  latter  into  (9). 

Finally,  we  note  that  (38)-(39)  are  easily  implemented  in  a  recursive  computer 
routine  to  determine  the  leftward-  and  rightward-looking  reflection  coefRcients  needed 
in  (40)  and  (41).  The  computations  proceed  from  the  outw’ard-looking  reflection  co¬ 
efficients  in  the  top  and  bottom  transmission  line  sections,  where  they  vanish  for 
unshielded  structures  or  otherwise  are  easily  determined,  providf'd  the  surface  admit¬ 
tances  >5,1  and  >s,jv+i  are  known.  Observe  that  the  exponential  functions  encountered 
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45 


in  (3S).  (40),  and  (41)  have  nonincreasing  magnitudes,  so  there  is  no  danger  of  over- 
tiow.  When  the  rth  transmission  line  section  extends  to  positive  (negative)  infinitv 

►  ^1 II <■ 

along  the  ;  axis,  the  reflection  coefficient  (  F,,)  is  set  to  zero  in  (40)  and  (41). 
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Abstract 

A  space-domain  mixed-potential  integral  ecpiation  approach  i<  applied  in  conjuiu  tion 
with  the  method  of  moments  to  compute  the  radar  cross-section  iH(’S)  of  coax- 
loaded  microstrip  patch  antennas  having  arbitrary  or  irregular  shapf's.  i  he  effects 
of  the  substrate  -  which  may  be  electrically  thick  and  may  consist  of  any  number  of 
planar,  possibly  uniaxially  anisotropic  dielectric  layers,  backed  by  a  ground  [)!ane  are 
rigorously  incorporated  in  the  analysis  by  means  of  the  vector  and  scalar  i>otential 
Green's  functions.  The  latter  are  expressed  in  terms  of  the  voltages  and  currents 
on  transmission  line  analogs  of  the  layered  medium,  associated  with  TM  and  I  K 
partial  fields.  The  current  distribution  on  the  microstrip  patch  is  approximated  usimt 
vector  ba^is  functions  defined  over  triangular  elements  and  the  coax  probi'  current 
is  expanded  in  terms  jf  piecewise-linear  subdomain  basis  functions.  .A  simple  prolu'- 
to-patch  attachment  mode,  compatible  with  the  triangular  element  model  of  the 
microstrip  patch,  is  used  to  enforce  current  continuity  at  the  junction,  and  the  coax- 
aperture  is  modeled  by  a  magnetic  current  frill.  The  far  zone  fields  are  found  by  iIh* 
stationary  phase  method,  and  are expres.sed  in  terms  of  the  Fourier-transformed  basis 
functions  and  the  transmission  line  voltages  and  currents  evaluated  at  the  stationary 
phase  point  value  of  the  transverse  wavenumber.  Computed  Ff CS  results  are  presented 
for  several  loaded  and  unloaded  microstrip  patch  antennas  of  various  shaf)es  an<l  are 
shown  to  be  in  agreement  with  published  measured  data  and  with  comf>uted  results 
obtained  by  specialized  techniques,  which — unlike  the  method  pre.serUed  lu're  are 
not  easily  extendable  to  arbitrary  shapes. 


1  Introduction 
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The  radar  cross-section  (RX'S)  of  rectangular  niicrostrip  {)atc!i  antennas  was  tirsl  --tiHl 
ied  by  Newman  and  Forrai  [1],  who  used  a  spectral  domain  integral  eijuatiun  [iictltud 
and  validated  it  by  measurements.  Pozar  [2]  used  a  similar  approarli,  t-xtimded  to  the 
case  of  a  uniaxial  substrate  and  combined  with  an  idealized  feed  nnxhT  to  compute 
the  RCS  of  both  loaded  and  unloaded  rectangular  patch  antentias,  .Jackson  kJ’  in 
vestigated  the  superstrate  effects  on  the  RCS  of  rectangular  microstrip  patches,  also 
using  a  spectral  domain  integral  equation  method,  .\berle  et  al.  [ll  cotnputeci  tin* 
RCS  of  loaded  and  unloaded  rectangular  and  circular  patch  antennas  usinsi  a  six'ctral 
domain  integral  equation  technique  incorporating  a  rigorous  feed  model,  and  corrtjbo- 
rated  their  analysis  by  measurements.  More  recently.  King  and  Bow  [.u|  used  a  similar 
approach  to  compute  the  RCS  of  finite  arrays  of  rectangular  rnicrosirip  patch  anten¬ 
nas.  The  RCS  of  cavity  backed,  loaded  and  unloaded  microstrip  patch  antennas  and 
antenna  arrays  was  studied  by  Jin  and  Volakis  [6]  using  a  hybrid  approach  combining 
finite-element  differential  and  integral  formulations. 

In  this  paper,  we  present  an  RCS  analysis  of  coax-loaded  microstrip  patch  antennas 
of  arbitrary  or  irregular  shape,  residing  in  a  grounded  dielectric  substrate  of  infinite 
lateral  extent.  The  analysis  is  based  on  the  mixed-potential  integral  equation  formu¬ 
lation  in  conjunction  with  the  method  of  moments  (MOM)  [7|  utilizing  a  triangular 
element  model  of  the  patch.  In  this  approach,  which  was  pioneered  by  Pichon  et  al.  l-'sj 
and  adopted  by  others  [9],  [10].  [11].  [12].  [13].  [14].  [15].  [16].  [17],  the  effects  of  the 
substrate  and  superstrate  (if  present)  are  rigorously  taken  into  account  by  means  of 
the  vector  and  scalar  potential  Green's  functions.  The  latter  are  <'X|)resse<i  m  terms 
of  the  voltages  and  currents  on  transmission  line  analogs  of  the  favored  medium,  as 
sociated  with  TM  and  TE  partial  fields.  The  current  liisi rib\itioii  en  iir  nm  lawtnf) 
patch,  which  may  have  an  arbitrary  or  irregular  shape,  is  appro\im..,e<l  m  terms  ot 
vector  basis  functions  defined  over  triangular  subdomains  lb'':,  I  he  '.ub'-trati*  ami 
the  cover  layer  (if  present)  may  be  electrically  thick  and  uniaxiallv  am^'otropic.  1  he 
current  on  the  coax  probe  is  expanded  in  terms  of  piecewise  linear  ''ulidomain  l)asls 
functions  and  the  coax  aperture  is  modeled  by  a  magnetic  current  frill  :19i.  [12:. 
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simple  probe-tO"patch  attaelu.u'iil  mode  il  1  i.  rompaiible  with  tlic  t ria.neislar  ol<>!n<  iit 
model  of  the  microstrip  patch,  is  impleinetited  to  otdorcc  cmreiil  cotitiiiuit}  at  i  he 
junction.  The  far  zone  fields  are  found  by  the  stationary  pha.s»*  method,  and  are  ex 
pressed  in  terms  of  the  Fourier-transformed  basis  fmicticuis  and  tin*  transmission  lute 
voltages  and  currents  evaluated  at  the  stationary  {jhase  point  value  uf  the  transverse 
wavenumber. 

The  remainder  of  this  paper  is  organized  as  follows.  In  St'ction  2.  a  (■<;Uj)h>d  s,u 
of  mi.xed-potentia!  integral  e<iuations  for  the  currents  induced  on  a  ( oax-loa<led  mi 
crostrip  patch  antenna  is  formulated.  Also,  the  coa.x  load  model  is  developed  and  the 
formulas  for  the  far  fields  and  F{(.'S  are  stated.  The  formulation  consistently  emjdov^ 
a  transmission  line  analog  of  the  layered  medium  and  is  thus  directly  applicable  to 
antennas  with  any  number  of  layers.  The  relevant  details  of  lh<“  transmission  line 
analysis  are  enclosed  for  easy  reference  in  .\ppendi.x  .\.  I  he  numerical  procedures 
are  described  in  considerable  detail  in  Section  3.  including  the  coax  probe-to-patch 
junction  treatment.  The  approach  relies  heavily  on  concepts  <ievelo[)ed  by  the  finite 
elements  community  to  efficiently  organize  the  computations  and  thi'  a-ssemb!y  of  the 
resulting  matrix  equation.  The  Fourier  transforms  of  the  shape  functions  encoun¬ 
tered  in  the  far  fields  computation  are  given  in  closed  form,  with  the  details  of  the 
development  relegated  to  .Appendix  B.  In  Section  4.  sample  computed  IK’S  results 
are  presented  for  several  loaded  and  unloaded  rectangular  and  circular  microstri[) 
patch  antennas  and  are  shown  to  be  in  agreement  with  published  measured  data  and 
with  computed  results  obtained  by  specialized  techniques,  which  -  tinlike  the  method 
presented  here — are  not  easily  extendable  to  patches  of  arbitrary  shape.  .\ew  R(’S 
results  are  also  included  for  a  coax-loaded  pentagonal  microstrip  patch  antenna  de 
signed  to  radiate  circularly  polarized  field  from  a  sdigle  feed.  Finally,  the  summary 
and  conclusion  are  given  in  Section  5. 

2  Formulation 

A.  Problem,  statement  and  assumptions 

The  geometry  of  the  problem  under  consideration  is  illustrated  in  Fig,  1.  where  the 
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Figure  1:  Geometry  of  a  coax-fed  microslrip  patch  antenna. 


dielectric  layers,  which  are  assumed  to  be  nonmagnetic,  are  characterized  by  their 
transverse  and  longitudinal  dielectric  constants,  tm  and  c,,,.  respectively,  relative  to 
free  space,  and  anisotropy  ratios  u„  =  t^nhtn-  If  dielectric  losses  are  present.  Cm  and 
t^n  are  multiplied  by  (1-  j  tand„).  where  tan<5n  is  the  loss  tangent  of  the  nth  layer. 
(Here  and  throughout  this  paper,  the  e-'"'  time  convention  is  implied.)  The  layered 
medium,  which  may  consist  of  an  arbitrary  number  of  planar  layers,  and  the  ground 
plane  are  assumed  to  be  of  infinite  lateral  extent.  The  coaxial  probe,  which  has  the 
inner  and  outer  radii  a  and  6.  respectively,  is  centered  at  (x.-,. (/,-)•  The  patch,  the 
ground  plane,  and  the  probe  are  assumed  to  be  perfectly  conducting. 

By  invoking  the  equivalence  principle  [20],  the  original  problem  of  Fig.  I  may  be 
replaced  by  its  equivalent,  shown  in  Fig.  2.  In  the  latter,  the  conducting  patch  and 
the  coax  probe  are  replaced  by  equivalent  electric  currents.  J,  and  /.  respectively. 
Also,  an  equivalent  magnetic  current  frill,  M,,  is  placed  over  the  coax  aperture,  which 
is  shorted.  To  simplify  the  analysis,  we  assume  that  the  distribution  of  the  aperture 
electric  field  is  that  of  the  TEM  coax  mode.  We  further  assume  that  the  coax  probe 
current  is  z-directed  and  is  azimuthally  invariant  on  its  circumference.  I'he  stn'ngth 


plane  wave 


(top  view) 


-ground  plane 
Figure  2:  Problem  equivalent  to  that  in  Fig.  i. 


of  the  magnetic  current  will  be  related  to  the  load  impedance.  Zi..  and  to  the  value  of 
/  at  the  base  of  the  coax  probe.  As  a  result,  the  patch  current,  J^.  which  resides  on 
the  surface  5  and  represents  the  vector  sum  of  the  currents  that  exist  on  the  bottom 
and  top  sides  of  the  microstrip  patch,  and  the  total  c-directed  coax  probe  current.  /. 
which  is  uniformly  distributed  on  the  perimeter  of  the  cylindrical  surface  C\  remain 
the  basic  unknowns  of  the  problem.  Once  these  currents  are  found,  other  quantities 
of  interest,  such  as  the  far  fields  and  the  RCS  of  the  antenna,  are  readily  determined. 


B.  Integral  equations 

Since  Green’s  functions  for  the  layered  medium  of  Fig.  2  are  available,  the  electric 
field  produced  by  the  coax-fed  patch  may  be  expressed  in  terms  of  integrals  over  the 
as  yet  unknown  currents  J,  and  /,  weighted  by  the  appropriate  kernel  functions.  The 
fields  excited  by  the  plane  wave  and  Af,  can  also  be  easily  determined.  .A  coupled  set 
of  integral  equations  for  J,  and  I  may  then  be  obtained  from  the  condition  that  the 
total  tangential  electric  field  must  vanish  on  the  patch  and  the  probe.  On  the  coax 
probe,  since  I  is  assumed  to  be  azimuthally  invariant,  we  only  enforce  the  vanishing 
of  the  *  component  of  the  electric  field,  averaged  over  the  probe  circumference.  In 
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what  follows,  we  use  the  following  notation:  r  is  the  position  vector  of  an  arbitrary 
point  with  respect  to  the  global  coordinate  origin,  p  is  the  projection  of  r  on  the  ./// 
plane.  is  the  radial  position  vector  of  the  coax  probe  axis,  primes  denott'  sourc«“ 
coordinates,  unit  vectors  are  distinguished  by  carets,  and  the  subscript  .s  denotes 
components  tangential  to  S. 

Upon  using  the  procedure  described  above,  we  arrive  at  a  coupled  set  of  intc'tiral 
equations  for  Js  and  I.  which  may  be  expressed  as  (the  details  of  the  derivation  are 
omitted  here  due  to  lack  of  space) 

^_a4(r|r')  J,(p')  +  V,/^Go(r|r')  V;.  J,(p')  d.s'' 

+  Qc I^I<,Ar\z)I(z')dz'  +  V,lj<^{r\z')  dz' 

-  Qc  Egir)  ^  Elir).  re>  ill 

-  eiiz)  =  e\[z).  z^C  (2) 

where  g.  =  p  —  Pc^  q'c  =  Pc~P'-  where  the  script  symbols  are  u.sed  tor  quantities 
averaged  over  the  coax  probe  circumference.  The  kernel  functions  in  ( l)-(2)  may  be 
expressed  as 

a4(r|r')  =  ^  ri-Hk,:z\zVo{lc.\p-p’\)k,dk,  (dj 

iTT  Jo 

G4r\r')  =  ~  rP[k,:z\z')Jo(k,\p-p'\)  ^  (4) 

1T\  Jo  fdp 

KpArW)  =  ^  rQ(K-,z\z%k{Ka)Jy{k,\p-pA)dk,  (b) 

ZTT J  Jo 

K^rW)  =  ~  rPikp,  z\z’)Jo{k,a)Jo{k,\p-pA)  (b) 

ITT  JQ 

Q.o{z\r')  =  ^  rR(k,;z\z')Mkpa)MkAp.-p\)di),  (Ti 

2wj  Jo 
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Oo{=\f')  =  —  f  P{ky.  :\:').Juik-^ti)Jo{k,\p.- p'])  ^  Oi 

k., 

K.,AA~')  -  " >■( k>-  k><‘ )  I '» ! 

Jtt  Jo 

A.'o(-|-')  -  ^  r  P(ky.=\z’)J^{kyi}  ^  nOi 

JttA)  k, 

in  which  denotes  the  Bessel  function  of  order  n.  kp  is  the  spectral  wavernirnber 
variable  corresponding  to  p,  and  where  t]Q  -  \Jpq/(o  and  k:^i  =  ^'y^pyo  are  the  intrinsic 
impedance  and  wavenumber  of  free  space.  The  Joikyi)  factor  appearing  in  fht'  ai)ove 
equations  is  the  result  of  averaging  of  the  respective  integral  kernel  over  the  t  oa.x 
probe  circumference  and  the  application  of  Graf's  addition  theorem  [21.  p.  .Ibd].  In 
(4)-(10).  we  have  introduced  the  auxiliary  functions 


S(kyzW) 


P(ky=\=')  =  V;'ikyz\z')-V:ikyz\=') 
Q{kp:z\z')  = 

R(kyz\z')^li^(kp-.z\z')-i:{kyz\z') 


M  nikyzlz')^ 


c,  +  e; 

(ko\'' 

\kj  _ 

m:=\=’) 


(12) 

(Id) 

(ID 


where  e,  and  e',  denote  the  longitudinal  dielectric  constants  of  the  observation  and 
source  layers,  respectively.  This  convention  is  in  effect  throughout  this  paper,  unless 
the  layer  index  is  explicitly  shown. 

It  is  important  to  note  that  the  integral  equations  (l)-(2)  are  in  the  rni.xed-potential 
form  [22],  [23],  [24],  [25],  which  is  amenable  to  the  existing  numerical  solution  proce¬ 
dures  developed  for  scatterers  of  arbitrary  shape  residing  in  a  free  space  [18].  [25]. 


C.  Transmission  line  network  analog  of  the  layered  medium 

As  an  aid  in  deriving  (3)-(10).  we  have  employed  a  transmission  line  network  analog 
of  the  layered  medium,  in  which  each  layer  is  represented  by  a  transmission  line  sec¬ 
tion.  as  illustrated  in  Fig.  3.  This  analog  compri.ses  two  networks,  which  arise  from 
the  decomposition  of  the  electromagnetic  field  into  partial  fields  that  are  transverse- 
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Note  that  the  integral  equations  (1)  (2)  are  vali<i  for  a  medium  with  an  arbitrary 
number  of  layers,  provided  the  appropriate  transmission  line  (Ireen's  functions  are 
used  in  {■^)-  { 10).  .\lso,  it  can  be  shown  that  the  kernel  functions  appearing  in  (  1 )  (2) 
exhibit  only  mild  singularities  when  the  source  and  observation  points  <oiiicide  on 
5  and  C.  which  makes  these  ecpialions  particularly  amenable  to  numerical  solution 
procedures. 

D.  Coax  load  model 

We  assume  that  the  aperture  field  in  the  problem  of  Fig.  1  is  that  of  a  TEM  coax 
mode  with  a  voltage 

=  2'l  4  (17) 

where  Z/,  is  the  load  impedance  presented  to  the  antenna  by  the  coaxial  transmission 
line  and  4  is  the  load  current.  In  view  of  this  assumption,  with  respect  to  a  polar 
coordinate  system  (o,  d)  centered  at  the  magnetic  surface  current  M.,  in  Fig.  2 

has  only  a  d  component,  given  as  [20,  p.  112] 


A4(o)=  -4-. 

0 


Zl 

\n.{b/a) 


in  which  a<  g<b.  The  electric  field  excited  by  this  magnetic  current  frill  embedded 
at  z  =  z'  in  the  layered  medium,  which  appears  in  (l}-(2).  can  be  expres.sed  as 


E^{r)  =  kT{p:z).  Stiz)  =  kU{z)  (19) 

where,  for  notational  convenience,  we  have  introduced  the  auxiliary  functions 

r(p;z)  -  K  rV:{K^z\z')\MKa)-Mk,b)]j,[k,\p-p^)dk  (20) 

JQ 

U{z)  =  4^  (21) 

The  second  Joikpo)  factor  in  (21)  results  from  the  averaging  of  El  over  the  coax 
probe  circumference  In  the  case  of  a  short-circuit  load  (Zi  =  0)  the  magnetic  current 
frill  is  absent.  In  the  other  extreme,  when  the  antenna  terminals  arc  open-circuited 
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(Z^,  =  cc ),  the  frill  is  absent  as  well,  but  we  impose  the  condition  4  =  0  in  t  lie  solution 
procedure.  F'or  any  other  load,  4  is  not  known  a  priori  and  must  be  computed 
together  with  and  /.  Note  that  (20)“(21)  are  in  a  general  form  applicable  to  a 
medium  with  an  arbitrary  number  of  layers,  provided  the  appro|U!ate  transmission 
line  Green's  functions  I’/  and  /'  ^^e  employed  (see  Appendix  A). 

E.  Incident  field 

The  structure  of  Fig.  1  is  excited  by  a  plane  wave  field 

JE;'"4r)  =  ^7*o[psinfi,cos(,.— ,.,)+(j-<nrosS,l 

incident  from  the  direction  on  the  grounded  layered  medium  from  the  upper 

half-space,  which  is  assumed  to  have  free-space  parameters.  Observe  that  in  (22)  the 
phase  reference  point  is  chosen  on  the  z  axis  at  z  =  d,  where  d  specifies  the  location 
of  the  uppermost  interface.  The  plane  wave  (22)  is  the  source  of  the  'incident'  field 
appearing  in  (l)-(2),  which  may  be  expressed  as 

Elir)  =  [p^'^k;-z\d)  +  ^,VHk;-,z\d)j  e^KP  (23) 

E'fiz)  =  -  ^sin0,4(A:;;c|d)  ('24) 

where  the  direction  of  the  unit  vectors  p,  and  and  =  p^kfi  with  k^  =  kosinO,. 
are  specified  by  the  direction  of  arrival  of  the  plane  wave.  In  the  above.  and 
denote,  respectively,  the  voltage  and  current  on  the  corresponding  transmission  line 
network  analog  of  the  layered  medium  (see  Fig.  3),  excited  in  the  uppermost  section 
by  the  incident  voltage  waves 

VT  ( ^ )  =  cos  e,  E'^^  ( 2.^1 ) 


_  j^inc  ^jko(z—ri)crtsi 


propagating  in  the  — r  direction.  In  (23)--('24),  we  explicitly  indicate  in  the  arguments 
of  and  P  their  dependence  on  k'^  and  on  the  phase  reference  point  {z  =  d)  of  the 


10 


61 


plane  wave  field.  (.Although  similar  notation  is  used.  an<l  /''  appearing  in  (2di 
(24)  should  not  be  confused  with  the  transmission  line  (ireen's  functions  introdur(>d 
earlier.)  The  Jo(A’^‘a)  factor  in  (24)  is  the  result  of  the  averaging  of  E[  over  the  coax 
probe  circumference.  Observe  that  (23)-(24)  are  in  a  gener  form  applicable  t<j  a 
medium  with  any  number  of  layers,  provided  the  appropriate  transmission  lin<’  xolt- 
ages  and  current  are  employed  (see  Appendix  .A). 


F.  Far  fields  and  RCS 

The  'scattered’  far-zone  field  radiated  by  the  patch  and  probe  currents  can  be  deter¬ 
mined  by  the  stationary-phase  method  [28].  It  is  then  found  that  the  patch  contri¬ 
bution  to  the  far  field  in  the  direction  is 


£e  ~  e-Ao(r-<fcos<t,)  .  f 

'2xjr  ^  Js 


■  n 


^’o 

Strjr 


,-jko{r-dcos9o) 


cos  e.  If  (A:;;  d\z')  Jfip')  e^^P'  dS'  (28) 


where  the  orientation  of  the  unit  vectors  po  and  (po-  and  k°  =  Pokp,  with  —  kosinO, 
are  specified  by  {6o^<~po)-  It  is  assumed  in  the  above  that  the  layered  medium  and  the 
radiating  structure  are  confined  to  the  region  z  <  d.  The  contribution  of  the  coax 
probe  to  the  far  field  is  found  by  a  similar  procedure,  with  the  result 


^9 


ko 

27rjr 


g-J^oir-dccsdo)  gJ<-0Pcsin5o  cos(v 


sin^„cos^.^  I  I{z’)n.{k^:d\z')dz' 
e'  Jc 


(29) 


where  {pc^c)  are  the  polar  coordinates  of  the  probe.  .Analogous  expressions  for 
the  far  field  of  the  magnetic  frill  current  can  also  be  given,  but  are  omitted  here 
because  the  contribution  of  the  coax  aperture — which  is  assumed  to  be  electrically 
small — to  the  total  radiated  field  has  been  found  to  be  negligible.  We  note  that  (27) 
(29)  are  applicable  to  multilayered  media,  provided  the  appropriate  transmission  line 
Green's  functions  are  employed  (see  .Appendix  A).  However,  these  simple  expressions 
lose  validity  near  the  horizon  [Oo  ~  90°),  in  which  case  the  stationary  phase  point 
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approaches  a  branch  point  singularity  in  the  A'p-plane.  and  guided  wave  phenomena 
dominate  the  far  field. 

If  the  plane  wave  field  incident  in  the  upper  half-.space  is  polarized  in  the  u  direriion 
and  the  c  component  of  the  scattered  field  is  considered,  where  n  and  c  stand  for  0 
or  the  RCS  of  the  antenna  is  given  as  [29,  p.  18] 


CTuv  = 


|£'.,.c!2 


(80) 


It  should  be  noted  that  the  scattered  field  in  the  above  does  not  include  the  plane 
wave  (geometrical  optics)  field  reflected  by  the  layered  medium. 


3  Numerical  Method 

A.  Weak  form  of  the  integral  equations 

The  MOM  is  applied  to  the  weak  forms  of  (1)  and  (2),  which  are  obtained  by  ’testing’ 
them  with  suitably  selected  weight  functions  {/lit}  and  {A’*},  respectively,  where  the 
former  are  defined  over  5  and  the  latter  over  C.  As  a  result,  upon  using  (lauss' 
theorem  [30,  p.  -503]  and  integrating  by  parts,  we  obtain  (cf.  [31]) 


where  we  have  introduced  the  notation  ( : )  for  an  integral  of  a  product  of  two  functions 
separated  by  the  comma.  The  dot  over  the  comma  signifies  a  ‘dot  product'  of  vector 
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arguments  and  the  prime  over  )  indicates  that  the  integration  is  over  tlu-  primed 
(source)  coordinates,  while  the  subscript  5  or  C  designates  the  domain  of  integration 
(i.e.,  the  surface  of  the  patch  or  the  length  of  the  coax  probe-  see  Fig.  2).  For 
simplicity,  the  arguments  of  the  integrand  functions  are  omitted  in  tdl  )  (d2). 

B.  Patch  current  expansion 

The  microstrip  patch  is  modeled  by  triangular  elements,  as  indicated  in  Fig.  2.  flie 
probe-to-patch  junction  can  be  located  anywhere  on  the  patch,  including  edges  and 
corners,  but  must  coincide  with  a  node  of  the  triangular  element  mesh.  The  nodes  of 
each  triangular  element  are  assigned  indices  /,  j,  and  h  in  a  counterclockwise  direction, 
as  illustrated  in  Fig.  4.  VVe  adopt  here  a  local  indexing  scheme,  in  which  these  indices 
assume  the  values  1,  2,  or  3,  in  a  cyclic  manner.  The  sides  of  a  triangle  vvith  an 

k 


Figure  4:  Local  coordinates  associated  with  a  triangular  element. 

area  A  are  formed  by  three  edge  vectors.  I,,  £j,  and  £>,.  where  £,  is  oriented  from 
node  j  to  node  k.  The  position  of  the  ith  node  with  respect  to  the  global  coordinate 
origin  is  specified  by  the  vector  r,,  whose  projection  on  the  xy  plane  is  p-.  Since  the 
microstrip  patch  lies  in  a  2  =  constant  plane  with  z  known,  the  location  of  an  arbitrary 
point  r  within  an  element  may  also  be  uniquely  specified  by  its  radial  position  vector 
p  =  P,  +  Qi,  where  g,  is  the  local  position  vector  originating  at  the  Th  node  of  the 
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element.  As  indicated  in  Fig.  4.  the  three  local  position  vectors  further  divid*‘  l  he 
element  into  three  triangles,  where  the  area  of  the  triangle  opposite  node  /  is  denoted 
by  .4,.  To  facilitate  the  integrations  over  irregularly  shaped  triangh's  encountered  in 
(31) -(32).  we  introduce  for  each  triangle  a  ’naturar  coordinate  system  (/.,.  /. ..  7.^  ). 
where  Li  is  known  as  the  ‘area  coordinate'  or  “shape  function"  associated  with  node 
i  of  the  element  [32.  p.  110].  and  is  defined  as 

''  1=1 

In  terms  of  the  area  coordinates,  the  local  position  vector  may  he  expicssed  as 

g,  =  tkLj  -ijLk  (34) 


To  represent  the  patch  current  on  each  triangular  element,  we  introduce  a  s  ector  basis 
function  ylj,  given  as  (cf.  [18]) 


Noting  that  [33] 


V,L, 


(36) 


where  n,  is  a  unit  vector  normal  to  edge  i  in  the  plane  and  pointing  out  of  the 
triangular  element  (see  Fig.  4),  we  find  the  divergence  of  A,  as 


V,  •  A, 


1 

A 


(31 


VVe  also  find  that  n,-  A,  is  constant  on  edge  t.  which  makes  it  easy  to  enforce  the  con¬ 
tinuity  of  the  normal  component  of  J,  between  the  elements  that  share  this  edge.  The 
patch  current  density  and  its  divergence  on  each  element  may  now  be  approximated 
as 

08) 

1=1  t=l 

where  /j  is  the  current  leaving  the  element  through  edge  i. 


14 


ft 


Proht  furi’t  nt  I  fpuii.<ion 

i  he  coaxial  [trobe  is  inodt'led  by  a  nuiidxT  ut  i!t.>!i  <n  crla[)|>!n‘;,  inie  ''i‘euicn!  «  Icnn-ii!  •• 
I  he  end  poinls  ol  a  .st'gmeiU  with  length  h  are  .e-signed  lu<  al  indicc*.  ;  and  /.  whii  ii 
assume  the  values  1  or  2  in  a  evciic  manner.  1  he  global  c<ajrdin<i! ol  tiu  'c  [i  Ui!- 
are'  speeiiied  by  and  As  indieale'd  in  b  ig.  tin'  di-'lance  td  a  point  ;■  v.  nin!;  an 


h, 


fw 


n. 


■^1 


I'igur*'  a;  Local  cuorilinate's  associal('<l  with  ;t  line  s('i>nieiii  elernent. 


element  from  node  j  is  denoted  by  //..  For  «‘a<  h  ('lenn'iit,  w<-  introduce  a  [eair  oi  litmar 
shape  functions  {.V.,-\j).  defined  as  jcf.  idl.  p.  119]) 
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The  location  of  a  point  c  on  an  elenn'iit  may  now  be  specitieei  as 


I  ddi 


C,hS: 


z- n, 


(  lOi 


where  n,  is  a  unit  vector  pointing  out  of  the  element  at  node  i.  Observe  tliat  the' 
value  of  Cl  is  either  +1  or  —1,  and  that  c,—  —  Frotn  (  !0i.  it  follows  that 
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The  axial  probe  current  and  its  derivative  on  an  eletnent  may  now  lie  appr(^xima!<'<i 
by  the  expansions 

i=l  1=1 

where  /,  is  the  current  leaving  the  element  through  node  i . 


D.  Attachment  mode 

The  current  expansion  (38)  is  not  suitable  for  a  patch  element  having  !h<'  probe  to¬ 
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patch  junction  point  as  one  of  its  nodes,  i'o  repre.sent  t  fie  current  dist  ril^ut  ion  mi  ■-iicli 
’junction  elements."  a  special  ’attachment  mode"  is  re<juired.  which  is  su|)erposed  on 
the  non-junction  part  of  the  current,  [(‘presented  by  (.iKi.  I,<-t  there  In'  .\"|  junctimi 
elements  attached  to  a  junction  node,  where  the  total  current  Ij  enters  th,e  microsirip 
patch  (see  Fig.  (i.  where  only  one  patch  junction  element  is  shown  for  simplicitv  i. 
Also,  let  the  junction  node  of  each  junction  element  be  assigned  a  local  index  /.  and 


k 


Figure  6:  Geometry  of  the  probe-to-patch  junction. 

let  the  angle  between  the  element  edges  that  meet  at  this  node  f)e  designated  a,. 
Furthermore,  let  the  sum  of  the  junction  verte.x  angles  be  denoted  by  a.  iGlearly. 
a  =  27r.  unless  the  junction  node  is  on  the  edge  of  the  microstrip  patch. )  To  reprc'sent 
the  junction  part  of  the  current  on  an  element  whose  node  i  is  the  junction  node.  w(‘ 
use  the  basis  function  [1’2] 

III  =  “ 

The  junction  current  density  and  its  divergence  on  this  element  may  then  be  expressed 
as 

Js  ~  Ij  III,  V,  •  J,  =  —  /;  “  f  11  i 

0.1 

where  the  second  expression  follows  from  (43)  and  (37).  The  union  of  the  currents 
represented  by  (44),  w'hich  only  exist  on  the  patch  elements  adjacent  to  the  junc¬ 
tion  node,  form  the  patch  part  of  the  attachment  mode.  As  mentioned  above,  this 
attachment  mode  exists  in  addition  to  the  non-junction  part  of  the  patch  current, 
represented  by  (38).  The  coax  probe  part  of  the  attachment  mode,  which  ('xteiids 
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over  thf’  line  st'gnient  adjacent  to  t  lie  juiu  tion  node,  ■'iinplv  mvcit  a'^ 


/  =  Ij.W , 


,/ 
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where  it  is  assunu'd  that  the  junction  node  ha.s  a  local  inde.x  i.  Ue  iioie  that  tiii' 
attachment  current  is  easilv  i!icor{)orated  in  the  ex[)an^iun  i  I'Ji.  hy  'ettinLt  /  -  //  tm 
the  junction  line  segment .  \\V  als(j  not<' that  therurrent  gict-ri  hy  i  i  i  i  doi-s  not  ciiif! 
the  patch  through  the  junction  node;  rather,  it  enters  througli  t  In*  elcnicrii  edo<*N  that 
meet  at  that  node.  Consequently,  unlike  suine  more  rigorous  atiachmeut  tnoih^  d')  . 
[36],  (41)  does  not  correctly  model  the  diverging  oirrent  hehavii>r  near  the  jutit  tiuii, 
.Nevertheless,  the  current  continuity  at  the  junction  is  satisfic'd  in  a  glohal  sense,  i.e,. 
the  total  current  entering  the  .\’j  junction  patch  elements  is  ecpial  to  the  cmrc'tit  /; 
leaving  the  coa.x  probe,  and  the  net  charge  associated  with  i  Mi  i  loi  is  /.ouj.  This 
attachment  mode  is  easily  implemented  in  the  .\1().\I  procedure-  because  it  (omprisc's 
regular  basis  functions  already  used  to  e.xpand  the  non-junction  |)ar1  of  the  t  urn'tit. 


E.  Glohal  MOM  matrix  a.i~ie.mbli)  and  solution  proctdiirt 

The  testing  functions  in  the  integral  equations  (31)  (32)  are  drawn  from  the-  '•ame 
.sets  as  the  basis  functions  used  to  repre,sent  the  current  on  the  [)atch  and  on  the  coa.x 
probe.  Hence,  {Ak }  consists  of  A\’^\  with  ?  =  1.  2,  and  3.  and  n  =  1 .....  NT.  where'  .N's 
is  the  total  number  of  triangular  patch  elements.  Similarly,  (.NT)  consists  of  c;’‘'.V,* 

with  ^  =  1,  2,  and  n  =  1 . Xc.  where  .\"c  is  the  total  number <jf  line  segment  (dement'^ 

on  the  coax  probe.  In  addition,  (31)  is  tested  with  the  patch  part  of  the  attachment 

mode,  comprising  the  union  of  the  bcisis  functions  /j]”*.  whe-re  n=  1.2 . Nj.  I  he 

resulting  equation  is  then  combined  with  the  equation  that  results  from  the-  tc'sting  of 
(32)  with  the  coax  part  of  the  attachment  mode.  where  thc'  nth  prolie  c-lc-ment 

is  assumed  to  be  adjacent  to  the  junction.  .As  before,  we  assume  here  that  the  local 
index  of  the  junction  node  is  i.  In  the  above,  it  was  nc'cc'ssary  to  introduce'  the 
superscript  (n)  to  distinguish  the  local  basis  functions  as.soriatc'd  with  elemc'nt  In 
what  follows,  where  there  is  no  danger  of  confusion,  this  element  .supc’rscri[)t  will  be 
omitted  for  notational  simplicity. 

When  the  expansions  (38).  (42),  and  (44)  (4o)  are  substitutc'd  into  the’  integral 
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equations  (.31)  (Ji),  the  coefficients  /*  *  are  constrained  by  the  houndar\'  conditions, 
which  require  continuity  of  the  normal  components  of  J,  at  ross  tlu*  edt^c's  shared 
by  adjacent  elements,  or  their  vanishing  at  tlie  boundary  edges  td  .s’,  as  w<dl  as 
continuity  of  the  coa.x  probe  current  between  adjacent  .segments.  On  the  npp<'rmosl 
probe  segment,  if  the  junction  node  has  a  local  index  i.  is  set  Uj  /j.  which 
explicitly  enforces  the  continuity  of  the  attachment  mode  current  at  the  [nobe-Uj- 
patch  junction.  Similarly,  on  the  line  segment  at  the  base  of  the  coax  [trobe.  if 
the  load  node  has  a  local  index  i,  /*"’  is  set  to  4.  In  general.  Ij  and  //,  are  not 
known  a  priori  and  must  be  computed  together  with  the  other  current  expansion 
coefficients.  In  the  unloaded  case  (when  Z^,  =  oci.  ’  is  set  to  zero.  If  the  number  of 
the  ncm-boundary  triangle  element  edges  is  denoted  by  Np  and  the  total  number  of 
the  unknown  current  coefficients  by  N,  then  N  =  +  l  (or  one  less  that  number, 

if  Zi  =  oc).  .As  a  result  of  this  procedure,  the  coupled  integral  equations  (dl)  -(d’2) 
are  converted  into  an  algebraic  system 


[Zmn]Un]  =  [Vm]  (4b) 

where  [Zmn]  is  the  .N-by-N  global  MO.\I  matrix  (also  referred  to  as  the  global  impe¬ 
dance  matrix),  [/„]  is  the  .\-by-l  global  vector  of  the  sought  after  current  expansion 
coefficients,  and  [Vm]  is  the  N-by-l  global  voltage  excitation  vector. 

To  assemble  the  global  system  (46),  we  consider  one  source  elernent-^test  element 
pair  at  a  time.  Let  the  global  indices  of  these  elements  be  and  n).  respectively. 
Then,  for  each  such  pair  we  assemble  a  local  system 

=  lb”’l  M7) 

where  [Z,,?"']  is  the  local  impedance  matrix  and  [Vi*""’]  is  the  local  voltage  excitatiun 
vector.  Here,  the  primed  local  source  element  indices  j'.  cind  (in  case  of  a  triangular 
element)  k'  follow  the  same  cyclic  convention  as  the  unprimed  ones.  Observe  that 
there  exists  a  unique  mapping  betw'een  the  local  coefficients  {/)'  *}  of  each  element 
and  the  global  current  coefficients  {/„},  where  the  reference  directions  of  the  latter 
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are  specified  by  the  order  in  which  the  element  nodes  appear  in  the  input  geometry 
data.  This  mapping  determines  to  which  entries  of  the  global  system  (16)  should  the 
elements  of  (47)  be  added,  and  with  what  signs. 

Consider  first  a  test  element,  5^,  on  the  microstrip  patch.  Then,  if  the  source 
element,  5^.  is  also  on  the  patch,  the  entries  of  the  residting  3-by-3  local  im|>edance 
matrix  are  found  as 


If,  on  the  other  hand,  the  source  element.  C„.  is  on  the  coax  probe,  there  results  a 
.3-by-2  local  impedance  matrix  with  the  entries  given  as 


In  addition,  when  Cn  is  the  line  segment  at  the  base  of  the  probe  and  L  is  the  local 
index  of  the  load  node,  this  matrix  is  augmented  by  the  load  contributions 

(50) 

which  only  affect  the  column  of  (49)  corresponding  to  the  load  node.  Finally,  the 
elements  of  the  3-by-l  local  voltage  excitation  vector  resulting  from  testing  over  Sm 
are  found  as 

=  (.51) 

Next,  consider  a  test  element,  Cm,  on  the  coax  probe.  Then,  if  the  .source  element, 
5„,  is  on  the  patch,  there  results  a  local  2-by-3  impedance  matrix  with  the  entries 
given  as 

y(jnn) 

^Xl' 

If,  on  the  other  hand,  the  source  element,  C„.  is  also  on  the  coax  probe,  the  entrif's 
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of  the  resulting  the  2-by-2  local  impedance  matrix  are  found  as 
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^11' 
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In  addition,  when  C„  is  the  line  segment  at  the  base  of  the  probe  and  /.  is  the  tucal 
index  of  the  load  node,  this  matrix  is  augmented  by  the  load  contributions 
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(mn) 
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which  only  affect  the  column  of  (53)  corresponding  to  the  load  node.  Finally,  tlu' 
elements  of  the  2-by-l  local  voltage  excitation  vector  resulting  from  testing  over  f 
are  given  as 


y(m) 


(551 


The  local  impedance  matrices  given  above  must  be  further  modified  to  include  the 
contribution  of  the  patch  part  of  the  attachment  mode.  Hence,  let  J  be  the  local 
index  of  the  junction  node  on  the  coax  probe.  Then,  for  every  patch  junction  element 
5n,  w'hose  junction  node  is  assumed  to  have  a  local  index  the  entries  of  (19)  are 
augmented  by 


and  the  entries  of  (53)  are  augmented  by 


.Note  that  (56)  and  (57)  only  affect  the  column  of  the  respective  impedancf'  matrix 
corresponding  to  the  junction  node.  Furthermore,  for  every  patch  junction  (dement 
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whose  junction  node  is  assunu'd  to  havt^  a  local  index  /.  tJic  inqx'dant  c  niatrix 
given  by  (52)  is  augnieiitecl  by 


^{rnn)  _ 


^A,:((;,.A, 


n.  l  \  .  5 /v., 


In  addition,  the  element  ZjJ"^  of  (55)  is  augmented  by  the  attachment  mode  ■-elf 


Z'r'-  ra.-a. 


^A.-./cU.-Ay) 

2q  \  2a  /s'„ 


2o  \  2n  /v„  /^. 

2a  \  2a  />,. 


O.d  \  •  0,1 /s. 


in  which  i'  is  the  local  junction  node  inde.x  on  element  Similarly,  for  every  patch 
junction  element  5m.  whose  junction  node  has  a  local  index  i.  the  imjjedance  matrix 
given  by  (53)  is  further  augmented  by 


and  bv 


Amn)  / 


in  which  L  is  the  local  index  of  the  load  node  on  the  coax  probe.  Also,  tin-  voltage 
excitation  vector  given  in  (55)  is  augmented  bv 


Note  that  (r)8)--{62)  affect  only  the  row  of  t  lie  res|>ect  i\'e  local  systou!  (  17  i  idi  ropond 
ing  to  the  junction  node  on  the  coax  probe. 

It  can  be  shown,  by  using  the  reciprocity  properties  of  the  t ransinissi(jn  line  ( iiccn'-' 
functions  (see  .Appendix  .-\),  that  *  apart  from  the  load  t*‘rnis  the  local  impedance 
matrices  defined  above  posses  symmetry  properties,  which  may  be  exploited  to  nearlv 
halve  the  computational  effort  involved  in  fillitig  the  glol)al  im|>edance  tnatrix  in  i  |t>i. 
The  integrals  over  the  soitrce  coordinates.  whi<'h  appear  in  the  impedance  inatij(es 
above,  involve  kernels  that  are  singular  when  the  test  and  source  eleitients  coincide. 
These  singularities  are  extracted  and  integrated  analytically  [381.  lea\  ing  well-ix-haved 
integrals  over  triangular  and  line  segment  elements,  which  are  numerically  evaluated 
by  Gaussian  quadratures  (see.  e  g..  [32.  p.  1 13]  regarding  Gaussian  rules  for  triangular 
domains).  On  the  other  hand,  ttie  testing  (exterior)  ititegrals  have  regular  and  slowly 
varying  integrands,  and  may  thus  be  approximated  using  one-point  (juadrature  rules 
[IS],  For  example,  the  integral  in  (51)  is  approximated  as 


ib3i 


where  and  denote,  respectively,  the  local  (with  respect  to  node  /  -^ee  Fig.  1 1 
and  global  position  vectors  of  the  centroid  of  element  S,rf  This  procedure  r('sults  in 
significant  savings  in  the  computational  effort,  even  though  it  sacrifices  the  symmetry 
properties  of  the  impedance  matrices. 

Once  the  complex-valued  matrix  equation  (16)  is  assembled  and  solved,  which  is 
accomplished  by  standard  procedures  (LT  factorization  with  partial  pi\oting,  followed 
by  a  forward  and  back  substitution  [37.  p.  160]).  the  current  density  within  each 
rnicrostrip  patch  element  may  be  obtained  from  (38)  (and  (14),  in  the  case  of  a 
junction  element).  Similarly,  the  current  within  each  coax  probe  element  may  be 
found  from  (42).  The  resonant  frequencies  and  Q-factors  of  a  microstrip  antenna  tna\ 
be  found  as  zeros  in  the  complex  frequency  plane  of  the  determinant  of  the  global 
impedance  matrix  (this  is  accomplished  by  means  of  the  Miiller's  search  procedure  [37. 
p.  T20|).  The  corresponding  modal  currents  may  then  be  obtained  from  (46)  with  the 
excitation  voltage  vector  set  to  zero. 


•>■> 


F.  Computation  of  spectral  integrals 

For  t!ie  solution  procedure  described  above  to  be  practical,  the  spectral  ii!t<-2rals  that 
occur  in  (3)-(10)  and  in  (20)-('2l)  must  be  ellicieutly  evaluated.  I'o  acconii)!ish  that, 
these  integrals  are  accelerated  by  asymptotic  integrand  subtraction  ami  the  method 
of  averages  [39],  [40|.  [41].  In  addition  to  these  techniques,  an  interindatiou  ami  tabh’ 
look-up  scheme  is  implemented  to  further  reduce  the  computation  time  [39;.  :  lOi.  IF. 
The  integration  path  is  properly  deformed  to  avoid  the  integrand  singularities,  whicli 
occur  on  or  near  the  real  a.'<is  in  the  Ay  plane  [14]. 

G.  Far  field  and  RCS  computation 

Once  the  coefficients  /,  are  found  for  each  element,  the  current  expansions  (38).  (42). 
and  (44)  are  substituted  into  (27)-(29)  to  determine  the  far  fields.  The  integrals 
encountered  in  (27)-(28)  are  then  recognized  as  Fourier  transforms  of  the  \ecior  basis 
functions  (35).  evaluated  at  In  view  of  (35)  and  (34).  the  Fourier  transform 

of  A,  associated  with  a  triangular  element  .8',,  may  be  expressed  as 

where  L,  denotes  the  Fourier  transformed  shape  function  I  pon  using  the  proce¬ 
dure  of  .Appendix  B.  we  may  express  /,,  as 

iy=l 

+  ^  JiifijkJ-l)  -  jjo{ajk'J2)  er^^^  Pj 

-  +  JJoiakkJ'l)]  eJ^'  P-  m) 

2Ay  I 

where  p,,  is  a  position  vector  of  the  midpoint  of  edge  i  of  the  element,  a,  -  u 
6,  =  v-i,.  and  j„  denotes  the  spherical  Bessel  [unction  of  order  n.  We  note  that 
although  (65)  can  be  put  in  a  simpler  form,  this  is  counterproductive  iiecause  it 
introduces  removable  singularities,  which  occur  when  happens  to  be  orthogonal  to 
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any  of  the  three  edges  of  tlie  element,  i.e.,  when  «,=0  for  /  —  I.  2.  or  d.  In  contrast, 
the  expression  for  Li  given  in  (65)  is  clearly  nonsingular,  unless  k),  =0.  Even  in  the 
latter  case,  it  can  be  shown  that  (65)  approaches  a  finite  limit  .1/5  as  Ay  — 0.  W’e  also 
note  that  when  (65)  is  evaluated  for  =  as  is  required  in  (27)  (2S).  then  Ay,  =  A' 
u  =  po.  and  v  =  (po- 

The  integrals  that  arise  when  the  probe  current  expansion  (12)  is  siibstitutecl  into 
(29)  depend  on  the  form  of  the  transmission  line  Creen's  function  /(.(A;,;  c/j;').  LcU  us 
suppose  that  the  coax  probe  extends  over  layer  /  of  thickness  f/;.  between  the  ;-axis 
coordinates  c;  and  c/+i.  where  <d.  Consequently.  /^^( Ay;  f/|r')  can  be  expressed  as 
a  constant  factor  times  /^(Ay:  -f+ilr').  Using  the  traveling- wave  form  of  the  latter  (see 
.Appendix  A),  we  then  find  that  the  integrals  over  a  line  segment  element  C,,  within 
layer  /.  required  to  evaluate  (29).  are  given  as 


c, 
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cuoihkli/2)+jMhky2) 


cjo{hk:i/2)  -  jji(hklJ2] 


(66) 


where  /;  is  the  length  and  the  midpoint  coordinate  of  tlie  element,  c,  is  defintHl  in 
(40).  A'!,  and  are  given  in  (16).  and  where  the  terminal  reflection  coefficients  F; 
and  r,  are  found  from  (71).  When  insed  in  (29),  the  above  is  evaluated  for  Ay  —  AC 
Once  the  microstrip  patch  and  coax  probe  contributions  to  the  far  zone  fields  are 
determined,  the  RCS  is  readily  obtained  from  (50). 


4  Sample  Results 

The  techniques  developed  in  the  previous  sections  have  been  implemented  in  a  FOR¬ 
TRAN  computer  code.  In  this  section,  we  present  sample  computed  monostatic  (i.e.. 
0„~0,,  -Po  =  •pi)  results  for  loaded  and  unloaded,  coax-fed  rectangular,  circu¬ 

lar  and  pentagonal  microstrip  patch  antennas  (see  f'ig.  7).  and  where  [)ossibl(' 
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Figure  7:  Geometries  of  probe-fed  (a)  rectangular,  (6)  circular,  and  (c)  pentagonal 
microstrip  patch  antennas. 

compare  them  with  available  published  data.  In  all  cases  included  here  the  patch 
resides  on  a  single- layer  substrate  of  varying  thickness  h.  without  a  cover  layer  (see 
Fig.  1).  The  incident  field  is  ^-polarized,  with  the  electric  field  amplitude  =  1  \’/m. 
All  RCS  results  are  plotted  vs.  the  frequency  /,  and  are  referred  to  1  m^  and  given 
in  dBsm  (‘decibels  above  a  square  meter’)  [29,  p.  160], 

The  results  shown  in  Fig.  8  are  for  an  unloaded  (Zf,  =  oo)  rectangular  patch  antenna 
(see  Fig.  7a)  on  an  isotropic  substrate  characterized  by  €r(l— jtan6).  with  and 
tan^  specified  in  the  figure  caption.  Newman  and  Forrei’s  [1]  and  Pozar's  [2]  results 
for  this  antenna  are  included  for  comparison. 

In  Fig.  9  we  show  results  for  a  coax-fed  rectangular  patch  antenna  (see  Fig.  7a, 
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where  {Xp.t/p)  are  the  feed  point  coordinates)  on  a  lossless  niiiaxial  siibstrat<\  I  hrce 
sets  of  data  are  presented,  corresponding  to  Z/,  =  oc  (no  load).  Z/,  =  bO  H  (inaicheil 
load),  and  Zl  =  0  (short-circuit  load).  For  the  first  two  cases,  our  r(>sults  are  compared 
with  those  of  Pozar  [2].  We  note  that,  as  e.xpecled.  the  first  and.  to  a  lesser  d<>gree, 
the  second  RCS  peak  are  suppressed  by  a  matched  load,  and  that  the  short  in"  [mu 
causes  the  RCS  peaks  to  shift  up  in  frequency. 


Figure  8;  RCS  of  an  unloaded  rectangular  microstrip  patch  antenna.  The  parameters 
are:  £  =  36.6  mm,  l'F  =  26mm,  £  =  1.58mm,  er  —  2.l7.  tan  (5  =  0.001.  .,^,=  45'’.  and 
61,  =  60“. 

The  results  shown  in  Figs.  10  and  11  are  for  coax-fed  circular  microstrip  patch 
antennas  (see  Fig.  7b,  where  Xp  specifies  the  feed  point  location)  on  isotropic  sub¬ 
strates,  with  the  loads  Zi  =  oc,  50,  and  OfT  In  Fig.  10  we  also  show  for  comparison 
Aberle’s  [42]  computed  results,  and  in  Fig.  11  we  include  the  results  computed  and 
measured  by  Aberle  et  al.  [4]  for  Zl  =  oo  and  5012.  In  reference  to  Fig.  11.  it  is  of 
interest  to  note  that  our  results  are  closer  to  the  computed,  rather  than  the  measured 
data  of  Aberle  et  al.  [4]. 

Finally,  in  Figs.  12  and  13  we  show  results  for  a  coax-fed  pentagonal  microstrip 
patch  antenna  (see  Fig.  7c,  where  P  denotes  the  feed  point)  on  an  isotropic  substrate. 
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f  (GHz) 

Figure  9:  RCS  of  a  rectangular  microstrip  patch  antenna  for  Zi  =  oc.  50.  and  OQ. 
The  parameters  are:  L  —  10mm,  W  =  15mm,  Xp  =  —L/A,  ijp  =  —W/A.  h  —  1.27 mm, 
Cf  =  13,  tz  —  10.2,  a  =  0.432mm,  b  =  1.397mm,  .p,  =  45°,  and  0^  =  0°. 

With  the  dimensions  given  in  the  caption  of  Fig.  12  (which  correspond  to  a  ca.se 
considered  in  [43,  p.  245])  and  no  load  attached,  the  first  three  excitation-free  modes  of 
this  antenna  were  found  at  the  resonant  frequencies  fri  =  1.123GHz.  frz=  1.194  GHz. 
and  /r3  =  2.184  GHz,  respectively.  The  modal  currents  of  the  first  two  modes,  which 
resonate  at  very  close  frequencies,  are  orthogonal  to  one  another.  With  a  properly 
selected  feed  l.jation,  these  two  modes  are  excited  with  equal  amplitudes,  resulting 
in  a  circularly  polarized  (CP)  radiation  field  at  a  frequency  that  lies  between  fri 
and  fr2  [43],  The  coax  feed  is  used  to  minimize  the  degradation  of  ellipticity  by 
unwanted  radiation  from  the  feed  network.  In  the  case  of  Fig.  12o.  where  p,  —  0°.  the 
polarization  of  the  incident  field  is  such  that  it  excites  the  first  mode  of  the  antenna, 
but  is  orthogonal  to  the  second  mode.  As  a  result,  the  co-polarized  RCS  component. 
agg.  dominates  the  cross-polarized  component,  crg.p,  except  near  the  CP  frequency,  at 
which  the  first  and  second  modes  are  coupled  through  the  loaded  coax  feed,  and  the 
value  of  (Tg^  approaches  that  of  rrgg.  The  second  peak  in  the  plot  of  rrop  corresponds 
to  the  third  mode  of  the  antenna.  In  the  case  of  Fig.  126,  where  p,  =  45°.  the  incident 
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F'igure  10:  RCS  of  a  circular  microstrip  patch  antenna,  (a)  Z/,  cc.  (b)  Z, 
and  Oil.  The  parameters  are:  R  —  Xp  —  9.2 mm.  h  —  1.58mm.  (r 

tan  6  —  0.0009.  a  —  0.4.32  mm,  b=  1..397  mm,  9,  =  0°.  and  0^  =  60°. 
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Figure  li:  RCS  ot  a  circular  rnicrustrip  patch  antenna,  ini  Z/  --  x.  \ht  Z/,  =  hi) 
and  Of).  The  parameters  are:  =  7.1  mm,  .r,,  —  2..”)  mm.  h  --  0.7^7  iTun.  ,  -  2.2. 

tan  0.0009.  a  =  0.  l22mm.  /)=1.297min.  .^.'.=0°.  and  th  =  r.2°. 
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tieUi  couplfs  approxiniateiy  ecpiaHy  Kj  tlic  lirsi  iwo  nuxic-.  hi  Im.  i;5  v.t-  -huw  '.cc  i.,t 
plots  of  the  real  ami  imaginary  i-iarts  of  the  j.iateh  (■urreiii  liensify  fur  Ui<‘  r,,'.*-  ol 
{■  ig.  I2u.  at  a  frecjuency  /  =  1 . 1!)  1  ( !llz,  which  cvirre-.pDnds  to  i  Ik-  lii  |/fak  of  i  Ik-  f7„ 
curve.  We  note  that  near  tlie  coax  |)rolK'  tin-  imaginary  jiar!  of  li-e  eurren!  lien^i!’,  e 
<lirectecl  towards  the  feed  point,  as  ex[)ecied. 

I'lu'  results  presented  in  I'  igs,  12  and  1-1  wen' ohi ain<‘d  usiim  a  ioti  elouient  niodei 
ol  the  pentagonal  patch  (which  can  he  seen  in  l-iu.  Ihi  and  a  i  hree  i-lenien!  niod<l 
of  the  coax  probe.  The  resulting  si/'e  of  the  <>iobal  impedance  maTix  m  >  ho  v,a.- 
22Q-hy  220.  The  computation  tirm-  was  umh-r  ill  minutes  fu-r  fref|u<*nc\  jioin’.  on  ai; 
elS6  PC  running  at  2o.\IHz. 

5  Conclusion 

We  have  developed  an  elegant  and  efficient  integral  equation  approach  for  the  RCS 
computation  of  coax-loaded  microstrip  patch  antenna.s  of  arlutrary  shape  with  sub¬ 
strates  that  may  be  electrically  thick.  Th<'  method  has  been  validated  for  rei  tammlar 
and  circular  patch  antenna.s,  for  whicii  result.s  are  available  in  the  literature.  Nh-w  K’(  'S 
results  have  also  been  presented  for  a  pentagonal  microstrip  patch  antenna  (iesigm'd 
to  radiate  a  circularly  polarized  field  from  a  single  coax  feed.  With  the  approacii  [ire- 
sented  here,  microstrip  patch  antennas  of  various,  [lossibly  irregular  shapes,  emln-d- 
ded  in  a  multilayered  uniaxial  substrate,  may  be  analyzed  within  a  singh'  theorefical 
framework,  using  the  same  computer  code. 

A  limitation  inherent  in  this  approach  is  the  assumption  of  a  laterally  iidinite  sub¬ 
strate  and  ground  plane.  .Although  the  patch  current  distribution  is  rather  insi-nsitive 
to  the  finite  substrate  and  ground  plane  effects  (assuming  tliat  0,  <  ')l)°i.  tin'  same 
cannot  be  said  of  the  scattered  field  (and  RCSl.  especially  for  0.  approac  hirsg  ftO'h 
where  the  diffracted  field  from  the  edges  of  the  substrate  and  ground  plane  may  be 
the  dominant  contribution  to  the  far  field  of  the  antenna.  .Also,  the  liark  lolies  of  ihe 
radiation  pattern  obviously  cannot  be  predicted  based  on  the  infmiti' gromnl  [)!aiK'  as 
sumption.  One  way  to  remedy  this  i.s  by  using  hybrid  techniques,  which  comlune  ilie 
infinite-substrate  integral  equation  method  with  the  geonu't rical  theory  of  diifraction 

:n 


((ITD)  [  t-l]  or  tho  Feist  Fourier  I'ransforru  (FFl)  inetliod  'ioi. 


A  Transmission  Line  Green’s  Functions 

Tlie  voltage  Ij  aiui  current  /,  at  a  |>oint  c  in  the  nth  transmission  line  '.ec  iiuu  (  (iii 
taining  a  unit-strength  current  source  i  at  z' .  as  illustrated  in  i'ig.  11.  the 

equations  [27,  pp.  7-47] 

[,{z\z')  f»7i 

az 

— =  -jk:.„)^,\]iz\z’]+d{Z~z'}  it)''! 

dz 

where  is  the  propagation  constant.  (Z,,)  is  the  characteristic  admittance  (inqie- 
dance).  and  S  denotes  the  delta  function.  For  notational  simplicity,  we  omit  here  the 
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Figure  1-4:  Transmission  line  section  comprising  a  unit-strength  current  sourn'. 

superscript  p  (which  stands  for  t  or  h).  and  we  do  not  e.xplicitlv  indicate  the  <h>pen- 
dence  of  the  voltage  and  current  on  the  transverse  wavenumber  k,,.  whic  h  determine 
k;n  through  (15)-(16).  From  (67)-(68).  Vifeje')  within  the  nth  transmission  lint 
section  containing  a  source  may  be  expressed  in  the  traveling- wave  form 
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wtiere  ~<  =in!n(c,c  ),  ;>  =iuax(c.c  ),  is  an  arbitrary  rcfcrfiict'  poip.!  uiihin  iri<- 
nth  section  (which  is  usually  set  to  c.  or  whichever  is  more  <  oiivt-nien!  i.  an>i 

I„(^o)  <tn(l  lri(~o)  are  the  rc'tleclion  coeHicients  ’lookins  t(j  ihe  left'  and  luDkim!  to 
the  right.'  respectively,  at  Cy.  By  means  of  the  translation  formula 


,,(;)  =  I'd 'o)t 


where  the  upper  and  lower  signs  (orrespond  to  the  right  and  left  arruW'-.  respe. 
tiveiy.  these  reflection  coefficients  can  he  expressed  in  terms  the  respt'cii\e  leirm 
rial  reflection  coefficients  F,,  and  F,,.  which  are  relati'd  to  the  corresponding  terminal 
admittances  and  (see  Fig.  FI)  as 


K  +  V  n 


For  a  transmission  line  section  of  a  finite  length  d^.  a  particularly  i  (invenient  form  of 
(69)  is 


v;(rb')  = 


\  jl+  F,  (  I  |_p  - 

2V,  1_  rr 


which  is  obtained  by  letting  =  and  using  (70). 

The  voltage  V’l  and  current  excited  by  a  unit-strength  voltage  source  r .  'satisfy 
equations  dual  to  (67) -(68).  which  areobtainerl  from  the  latter  by  making  the  subsii 
lutions;  V,  — >  .  /,  — *•  V).,  Zn  — »  V„.  and  >k  — +  Z„.  Furtlu’rmorm  it  can  be  shown  that 
the  following  symmetry  and  reciprocity  relations  hohi  (cf.  ;16.  p.  7F.  (27.  p.  lOlF 


V:izW)  =  V\{z'\z).  /,dc|Z)  =  /,(c'ic).  Fd-jZl-  -/pr' 


.•\s  a  result.  li.{z\z')  ran  be  obtained  from  (72)  by  replacing  in  the  lattei  >,i  by 
inverse  (which  causes  the  reflection  coefficients  to  change  signs).  /,!  :k')  follows  from 
(67)  and  (72),  and  V',.(~|::')  may  then  be  obtained  from  the  last  reiation  in  f7d!. 

The  voltage  and  current  on  the  nth  transmission  line  section  that  is  sonri  ('  free 


satisfy  tlu'  iKMiiogeneous  tortn  ot  i(>7i  (tiSj.  l-roni  thcM'  i'(|i!,ti  !'n<'  a;  a.n\ 

point  c  within  llio  line  sertioti  may  !)»•  oxpn“s.st*<i  in  tfrm‘>  of  ilir  \olt,ii>o  \  >  .  a<  no- 

one  of  its  terminal  pairs.  The  resnli  is 

\-(ci  =  I  io.)'  —  fi^  t;  f  -  1  ,7!, 

I  r-  17  (  ^ 

where  ~o  =  c„,  or  co  =  ;>i  +  i  0*‘(‘  Fig.  1  1).  liepeixlirig  on  wiietiiei  i  lio  -oiir<  c  )-  loratei!  to 
the  left  or  to  the  right,  res[)eetively.  of  the  line  section,  ami  where  t  ho  upper  t  lowers 
arrow  corresponds  to  c  >  ~o  (c  <  rij).  We  have' omitted  the  siiF'-cript  of  I  sn  iflc 
because  the  latter  applies  irrespect iv»‘ of  the  naiureo:  the  source,  and  dro[>j!eii  fiom 
its  argument,  because  (71)  only  implicitly  d<*[)<-nds  on  the  source  location,  which  !> 
outside  the  «th  line  section.  The  current  /(cl  corresponding  to  t  7!  i  may  be  obtained 
by  substituting  the  latter  into  (()7). 

Finally,  we  note  that  (70)  (71)  are  easily  implemeutcHl  in  a  recursive  r(/rnpuier 
routine  to  determine  the  leftward  and  rightward  reliection  coelfi*  ient-'  needc-d  in  f72i 
and  (74).  The  computations  proceed  from  the  outward-looking  rcdlc'ctioti  (  ^>e|hcie^t^ 
in  the  top  and  Itottom  transmission  line  sections  (s<>e  Fig.  4).  ObscTve  that  the  ex[)o- 
nential  functions  encountered  in  (70).  (72).  and  i74)  lia\e  noniticreasing  magnitudes, 
so  there  is  no  danger  of  overflow.  Wlien  the  nth  transmission  line  mtIiou  extemls  to 
positive  (negative)  infinity  along  the  -  axis,  the  reflection  coefficient  1'.;  ;  F,  >  is  mU  to 
zero  in  (72)  and  (74). 

B  Fourier  Transform  of  a  Triangle 
Shape  Function 

The  Fourier  transform  of  a  shape  function  L,  associatc'd  with  a  fnangular  element 
is  given  as 

Z,  =  /  P  (IS  i  7bi 

where  kp  —  ukp.  To  evaluate  the  integral  in  (75).  it  is  helpful  to  first  convert  it  to 
a  line  integral  around  the  boundary  contour.  ()S,,.  of  .S,  ('f.  (17;.  jlsj.  ildji.  ibis  is 
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most  oasilv  accomolishod  hv  uotiii"  tliat  mOi 


p  =  T,  •  (  u 


,  P  , 


jk,  y 

and  l)y  making  use  of  tlie  <!ivergen<e  ilieorem.  As  a  n'snlt.  ulien  k  ill,  \xf  (ihiaHi 


I  ,  I. 
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where  v  =  z  x  u  and  n  denotes  a  unit  vector  tiormal  t(j  ()S,,  at  i  in  iIh-  plane  <ind 
pointing  out  of  >’,1.  The  integral  in  (77)  is  easily  evaluated  in  tlie  !(ical  coordinates, 
if  one  notes  that  when  t  is  on  edge  t.  n  =  r^^.  0  <  <1.  /.,  .=  I).  /. ,  =  1  -  .  arid 
dk  —  f^dLk  (see  Fig.  1).  .\lso.  in  that  rase  p  =  p,^  Q,-  where  is  a  x'ecitu'  fruin  node 
j  to  the  point  Z  on  edge  i  of  the  element.  The  resulting  closed  form  e.xpression  for  /., 
is  given  in  (60). 
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ABSTRACT 

An  arbitrarily  shaped  microstrip  patch  antenna  excited  through  an  arbitrarily  shaped 
aperture  in  the  mouth  of  a  rectangular  waveguide  is  investigated  theoretically  and 
experimentally.  The  metallic  patch  resides  on  a  dielectric  substrate  grounded  by  the 
waveguide  flange,  and  may  be  covered  by  a  dielectric  superstrate.  The  substrate  (and 
superstrate,  if  present)  consists  of  one  or  more  planar,  homogeneous  layers,  which 
may  exhibit  uniaxial  anisotropy.  The  analysis  is  based  on  the  space  domain  integral 
equation  approach.  More  specifically,  the  Green's  functions  for  the  layered  medium 
and  the  waveguide  are  used  to  formulate  a  coupled  set  of  integral  equations  for  the 
patch  current  and  the  aperture  electric  field.  The  layered  medium  Green's  function  is 
expre.ssed  in  terms  of  Sommerfeld-type  integrals  and  the  waveguide  Green’s  function 
in  terms  of  Floquet  series,  which  are  accelerated  to  reduce  the  computational  effort. 
The  coupled  integral  equations  are  solved  by  the  method  of  moments  using  vector 
basis  functions  defined  over  triangular  subdomains.  The  dominant  mode  reflection 
coefficient  in  the  waveguide  and  the  far  field  radiation  patterns  are  then  found  from 
the  computed  aperture  field  and  patch  current  distributions.  The  radar  cross  section 
(RCS)  of  a  plane-wave  excited  structure  is  obtained  in  a  like  manner.  Sample  numer¬ 
ical  results  are  presented  and  are  found  to  be  in  good  agreement  with  measurements 
and  with  published  data. 
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1  Introduction 

Microstrip  patch  antennas,  which  belong  to  a  large  class  of  printed  circuit  aiuennas. 
are  widely  used  in  the  microwave  fre<iuency  range,  both  as  single  elements  ami  tmore 
often)  in  array  configurations.  Their  advantages  are  well  known;  low  cost,  (  onfomiity. 
ease  of  fabrication  and  integration,  reproducibility,  ruggedness,  light  weight,  and  low 
profile  [1],  [2].  In  recent  years,  the  utilization  of  millirneterwave  systems,  with  smaller, 
lighter  components  and  antennas,  has  provided  a  wider  bandwidth,  and  conseciuently 
higher  data  rate  communication  and  better  resolution  than  microwave  systems  [3], 
[4].  However,  the  feed  structures  that  operate  w'ell  at  microwave  frequencies  are 
not  always  viable  in  the  millirneterwave  range.  For  example,  microstrip  line  losses 
become  significant  and  coaxial  feed  components  are  not  available  above  about  50  G  Hz. 
which  renders  the  direct  feed  techniques  impractical  in  this  frequency  range.  On 
the  other  hand,  since  the  waveguide  bulkiness  becomes  less  of  a  factor,  while  its 
losses  remain  smaller  than  those  of  a  microstrip  line  [5],  the  indirect  waveguide  feed 
becomes  an  attractive  option  for  millirneterwave  antennas  and  antenna  arrays  [6]. 
Moreover,  in  some  applications  the  aluminum  waveguide  may  also  serve  as  a  heat 
sink  and  support  for  active  devices  that  may  be  integrated  with  the  antenna.  Another 
antenna  configuration  that  possesses  the  advantages  listed  above  is  a  waveguide-fed 
slot  antenna.  Its  gain,  however,  is  significantly  lower  than  that  of  a  waveguide-excited 
microstrip  patch. 

In  this  paper,  we  present  a  rigorous  integral  equation  analysis  of  a  waveguide- fed 
microstrip  patch  antenna,  as  illustrated  in  Fig.  1.  The  arbitrarily  shaped  microstrip 
patch  is  coupled  to  the  rectangular  waveguide  through  an  aperture,  which  may  also 
be  of  arbitrary  shape.  The  dielectric  medium  above  the  ground  plane  may  consist  of 
one  or  more  planar,  homogeneous  layers,  which  may  exhibit  uniaxial  anisotropy.  The 
nth  layer  is  characterized  by  the  permittivity  and  permeability  dyadics 

^  ~  I^^tn  T  M  ~  T  ZZfl^n  (1) 

where  is  the  unit  dyadic  transverse  to  z,  and  tm  ifJ'tn)  and  t^n  (M^n)  denote,  respec¬ 
tively,  the  transverse  and  longitudinal  dielectric  (magnetic)  constants  relative  to  free 
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Figure  1:  Geometry  of  a  waveguide-fed  microstrip  patch  antenna. 


space.  Observe  that  we  distinguish  dyadics  by  double  underlines  and  unit  vectors  by 
carets.  The  free  space  permeability  ajid  permittivity  will  be  denoted  by  Ho  and 
respectively.  For  each  layer,  we  also  introduce  its  electric  and  magnetic  anisotropy 
ratios,  and  respectively,  given  as 


^zn 

^in 


f^zn 

f^tn 


(2) 


The  structure  is  excited  either  by  the  dominant  (TEio)  waveguide  mode  or  by  a  plane 
wave  incident  in  the  upper-half  space.  The  primary  quantities  to  be  computed  are 
the  aperture  field  and  patch  current  distributions.  From  these,  other  quantities  of 
interest  will  be  found,  including  the  dominant  mode  reflection  coefficient,  the  far  field 
radiation  patterns,  and  the  radar  cross  section  (RCS).  Although  attention  is  limited 
to  a  single  antenna  element,  it  is  expected  that  the  results  of  this  study  will  also  be 
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useful  in  the  analysis,  by  an  approximate  teehnique,  of  finite  waveguide  fed  auiet da 
arrays. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Sec.  2  we  formulate  <  <aip!(>d 
integral  equations  for  the  patch  current  and  the  ap<>rture  electric  fi»‘h!.  where  the 
latter  is  represented  by  an  equivalent  magnetic  current.  1  he  kernels  of  th<-se  integral 
equations  arc  expressed  in  terms  of  the  voltage  and  current  (Ireen's  functions  of 
a  transmission  line  network  analog  of  the  layered  medium,  which  is  discussed  in 
Appendix  A.  In  Sec.  3  we  give  the  numerical  procedure  for  the  solution  of  the  integral 
equations,  and  in  Sec.  I  we  describe  the  experimental  verification  of  the  theory.  VVe 
present  sample  computed  and  measured  results  in  Sec.  5,  and  give  conclusions  in 
Sec.  6. 

2  Formulation 

2.1  Integral  equations 

To  facilitate  the  analysis  of  the  structure  shown  in  Fig.  1,  we  invoke  the  equivalence 
principle  [7]  to  in  effect  decouple  the  original  problem  into  two  simpler  ones,  referred 
to  as  the  ‘interior’  (inside  the  shorted  waveguide)  and  ‘exterior’  (above  the  ground 
plane)  problems,  as  illustrated  in  Fig.  2.  This  decoupling  is  achieved  by  first  shorting 
the  aperture  Sa  through  which  the  two  regions  interact,  and  then  placing  over  it  an 
equivalent  magnetic  current  Ms,  which  represents  the  tangential  electric  field  in  the 
aperture.  The  negative  of  Ms  is  placed  on  the  opposite  side  of  the  shorted  aperture, 
thus  explicitly  enforcing  the  continuity  of  the  tangential  electric  field  across  5^.  In 
the  equivalent  problem,  the  effect  of  the  metallic  patch  Sp  is  replaced  by  an  equivalent 
electric  current  J5. 

The  conditions  of  vanishing  tangential  electric  field  on  the  patch  (which  is  assumed 
to  be  perfectly  conducting)  and  continuity  of  the  tangential  magnetic  field  across  the 
aperture  may  now  be  stated  as 

z  X  [jE;;(r)  +  £;v(r)]  =  O,  r^Sp  (3) 
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Figure  2:  Problem  equivalent  to  that  in  Fig.  1. 

i  X  [H;(r)  +  H\{r)]  =  i  x  [Hl{r)  +  ,  r  6  (-1) 

where  the  subscripts  +  and  —  refer  to,  respectively,  the  exterior  and  the  interior 
regions,  r  is  the  position  vector,  {E\H')  are  the  the  ‘short-circuit'  electric  and 
magnetic  fields  of  known  sources,  computed  in  the  absence  of  the  patch  and  with  the 
aperture  shorted  [8],  and  {E^,  H“)  are  the  ‘scattered’  fields  produced  by  Js  and  Ms- 
Rearranging  (3)  and  (4),  we  obtain 

-  z  X  E\,{r)  =  z  X  E\(r) ,  r  ^  Sp  (.o) 

-  i  X  (H;(r)  -  HHt)\  =  i  X  [HV(r)  -  .  r  €  .%  (6) 

The  scattered  fields  in  (5)-(6)  can  be  expressed  in  terms  of  the  dyadic  and  scalar 
kernels  of  the  exterior  and  interior  problem  as  [9] 

=  -  f  g:'{r\r')-Jsip')dS'-vf  G^ir\r')V'-Js{p')dS' 

JSp  JSp 

+  /  Q^^{r\r')-Ms{p')dS' 

J  Sei 

=  I  G‘^-'{r\r')-Js{p')dS'  -  I  G^Hr\r')- Ms{p')dS' 

JSp  JSa 


(7) 
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V  /  V'-A/s(p'bi.s'' 


I  ^ ) 


HUr)  =  l_  G‘''-{r\r')-M<{p')dS' 

+  V/  a'^-(r\r')V'-M,{p')dS' 


CM 


where  primes  indicate  source  coordinates  and  p  is  the  projection  of  r  on  t  lie  r //  phme. 
Upon  substituting  (7)  (9)  into  (5)  (6),  we  obtain  a  coupled  set  of  integral  e<juations 
for  the  unknown  equivalent  currents  J.s  and  Ms,  given  as 


+  V,  (G"  ,  -  (Sf ' :  =  E;.,(r),  r  e  ( 10) 

(a"^ ;  -  (a"*  +  a''-  ^  m.); 

-V,  +  G^-  ,  V;-Ms)'  =  -  H’Jr),  r  e  .%  (11) 

Here,  we  have  introduced  the  notation  (;  )  for  an  integral  of  a  product  of  functions 
separated  by  the  comma.  The  dot  over  the  comma  signifies  a  'dot  product'  of  vector 
arguments,  while  the  prime  over  )  indicates  that  the  integration  is  over  the  primed 
coordinates.  The  subscript  t  distinguishes  components  transverse  to  z.  and  the  sub¬ 
scripts  Sp  and  Sa  indicate  the  integration  domains  (see  Hg.  ‘2).  The  dyadic  kernels 
in  (lO)-(U)  are  given  as  [9] 


Sf(rl'')  =  LG^Jrlr') 


(12) 


Gf"(rlr') 


=  (ii-5r!()G“(r|r') 
+i!(G™(r|r')  +  »iC“(r|r') 


(13) 


0  (rjr')  =  (-XX  +  yy)G;;'^(r'jr) 

-iifG^f(r'lr)  -  yxGf^'''(r'|r) 

Gfnr|r')  =  /Gf^(r|r') 


(11) 

(15) 
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when* 


Sf  "(»■  i  »•')  =  xi'('»’[~(r  j  r')  f  yy(!‘^:^  ir  j  r') 
r;^,(r|r')  =  .S)(v;'‘(Av:cjc'!} 
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&■"'(■•  I'-') 


sin2C  J 

\V:{k,-\z')^V:^{k,::\:') 
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[  ^ 
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2  ‘M 

[  ^ 
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t  1  -t 
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^  mn  ^  •'  mn  •'  mn  mn 


and  the  scalar  kernels  are  given  ais 

C*(r|r')  =  5o 


K 
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ItH 


( 


( i'M 


(21) 

(22) 


(23) 


(24) 

(25} 


G"-(r|r')  = 


jr(/)  +  jr(%" )  .  4. 

mn  '  mn  *  ^  mn  4  mn 


In  the  above,  we  have  introduced  the  Sommerfeld  integral  notation 
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^  =  yTr  -  +  (i/  ;/)•■  C  =  iin  taii 

where  Jn  is  the  Bessel  function  of  order  n.  and  the  Moquet  series  notation 

E  E  ^-"11 

=  —  >0  =r  —  X.’ 

with 

^fimn  “  ^'x,n®  +  =  \J  ^r,r,  "h  ~  (77  ^  ' 

Here,  =  d^dy.  with  dr  —  2a  and  dy  =  26,  where  a  and  6  are  the  inner  waveguidt' 
dimensions  along  the  x  and  y  directions,  respectively.  The  symbols  h'f  and  /''  m  M?) 
(26),  where  p  stands  i  a  e  {E  mode)  or  6  (//  mode)  and  a  stands  for  i  (current  sour<  e) 
or  V  (voltage  source),  denote,  respectively,  the  voltage  and  current  transmission  line 
Green’s  functions  discussed  in  Appendix  A.  The  four  series  in  (22)  (231  and  (26)  arise 
from  the  multiple  images  of  a  transverse  magnetic  dipole  in  the  walls  of  a  rectangular 
waveguide. 

Once  the  current  distributions  Js  and  Ms  are  determined,  the  secondary  (}uan- 
tities  of  interest,  such  as  the  waveguide  dominant  mode  reflection  coefTicient  at  the 
aperture,  the  far  zone  field  radiation  pattern  and  the  RCS,  can  be  found  with  litth' 
extra  effort. 

2.2  Short-circuit  incident  fields 

In  the  driven  antenna  analysis,  the  incident  electric  field  in  the  air-filled  waveguide  is 
taken  to  be  that  of  the  dominant  TEio  mode  (we  assume  that  a  >  6).  and  is  given  a.s 

=  14610(^)6-^*^-°'  (31) 

with  _ 

e,o(p)  =  wy^cos(^^)  ,  -  (t)  02) 

where  ko  =  ^Jpo(-o  is  the  free  space  wavenumber.  This  incident  wave  excites  the 
short-circuit  magnetic  field,  whose  transverse  component  is 

H\  ^r]  =  [l-  fwiz)]  (33) 


4'N/(-4)} 
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with 
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where  is  the  free  spae«'  intrinsic  impedance,  in  1  i„(-i  i'- 

dominant  mode  reflection  coefficient  in  the  short-circuited  waveguide.  wlii<  h  a.sMiines 
the  value  —  1  at  -  0. 

In  the  RCS  analysis,  the  structure  is  illuminat»’d  by  a  uniform  plane  wave  iicdd 


E+‘^{r)  =  [d,E'g  +  <i>,  K'J)  f  cost V' cos i.i  .j- , 

impinging  from  the  direction  on  the  grounded  layered  medium  from  the  upper 

half-space,  which  is  assumed  to  have  free  space  parameters.  Observe  that  in  (d')) 
the  phcLse  reference  point  is  chosen  at  (0,0, d),  where  d  is  the  c  — coordinate  of  the 
uppermost  interface.  This  plane  wave  is  the  source  of  the  short-circuit  electric  field, 
whose  transverse  part  is  given  as 

£;.+(<•)=  (:)6) 

where  fc*  =  p^k^,  with  k^  =  koS'mOi.  In  the  above,  V^{ky,  z  \d)  denotes  the  voltage 
on  the  respective  transmission  line  analog  of  the  layered  medium,  excited  in  the 
uppermost  section  by  the  voltage  waves 


Vl{k;-,  z\d)  =  cos 0,  El  (37) 

( /fc; ;  2 1  d)  =  ( 38 ) 

propagating  in  the  —z  direction.  The  transverse  short-circuit  magnetic  field  associ¬ 
ated  with  (36)  may  be  expressed  cis 

HU(r)  =  -  [p,IHK.  ^  I'')]  I®) 

The  voltages  and  currents  in  (36)  and  (39)  can  easily  be  found  in  any  transmission 
line  section,  as  discussed  in  Appendix  A. 
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2.3  Secondary  quantities 


From  the  equivalent  currents  and  1V/5.  other  quantiti<‘s  of  interest  nia\  easily  l>e 
obtained.  Hence,  for  the  driven  antenna  we  compute  tlie  dominant  niod.e  .tperiure 
reflection  coefficient  F  and  the  corresponding  normalized  aperture  admit tanee  V  . 
given  as 

r  -  ~  1  .  V  =  — 4:  (If)i 


For  the  plane-wave  excited  structure  we  also  compute  the  dominant  mode  aperT3ire 
transmission  coefficient,  given  as 

{Ms',  ^io)s 


T  -  — 


(in 


’<  1-2 


ymsp  +  TE 

For  either  excitation,  the  far  zone  fields  in  the  direction  {0o,'Po)  ^^re  obtained  a.s 


[  («) 

where  and  are  the  radial  and  azimuthal  unit  vectors  evaluated  at  {Oa-  po),  and 
k°  =  Pgkp,  with  k°  =  Atosin^o.  It  is  assumed  in  (42)'-(43)  that  the  aperture  is  in 
the  z  =  0  plane  and  the  patch  resides  on  a  substrate  with  thickness  h.  Observe  that 
(42)-(43)  include  contributions  both  from  the  patch  and  the  aperture. 

For  the  plane-wave  excited  antenna  wc  also  compute  its  monostalic  RCS.  Hence, 
if  the  incident  field  in  the  upper  half-space  is  polarized  in  the  u  direction  and  the  c 
component  of  the  scattered  field  is  considered,  the  RCS  of  the  antenna  is  given  as  [9] 

It  should  be  noted  that  the  scattered  field  in  the  above  does  not  include  the  plane 
wave  (geometrical  optics)  field  reflected  by  the  layered  medium. 
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3  Numerical  Procedure 

3.1  Weak  form  of  the  integral  equations 

The  method  of  moments  (MOM)  is  applied  to  the  weak  forms  of  { 10)  and  I  1 1  j.  wlii<  h 
are  obtained  by  testing  them  with  suitably  selected  weight  functions  deiined 

over  Sp  and  As  a  result,  upon  using  the  Gauss  theorem  [10.  p.  50;$|.  we  obtain 

(  A„  ;  (Gi  .  -  (v. .  A„  .  (G-  .  v; . 

-  (a„  ;  (gf  '> ;  ^ 


An. ;  (a'"  ;  -  (a„  ;  (Gfji,  +  Sf  -  ; 

+  (v,- Am,  +  O'"-  .  v;  ■  Ms)[  ) 

\  /Sa 

=  (a„  ;  -  H\\ 


Here,  the  incident  fields  E\.^_  and  are  absent  when  the  waveguide  aperture 

reflection  coefficient  is  calculated,  and  H\  _  is  zero  in  RCS  computations. 


3.2  Patch  and  aperture  current  expansions 

As  indicated  in  Fig.  1,  we  model  the  microstrip  patch  Sp  and  the  aperture  .Ti  by 
triangular  elements.  The  nodes  of  each  triangular  element  are  assigned  indices  i.  j. 
and  fc  in  a  counterclockwise  direction,  as  illustrated  in  Fig.  3.  The  sides  of  a  triangle 
are  formed  by  three  edge  vectors,  li,  Ij,  and  ik,  where  i,  is  oriented  from  node  j  to 
node  k.  The  position  of  the  ith  node  with  respect  to  the  global  coordinate  origin  is 
specified  by  the  vector  ?•;,  whose  projection  on  the  xy  plane  is  p,.  Since  the  rnicrostrip 
patch  and  the  aperture  each  lie  in  a  ^  =  constant  plane  with  c  known,  the  location 
of  an  arbitrary  point  within  an  element  may  uniquely  be  specified  by  its  radial  vector 
P  —  P,  +  6ti  where  q-  is  the  local  position  vector  originating  at  the  ith  node  of  the 
element.  As  indicated  in  Fig.  3,  the  three  local  position  vectors  further  divide  the 
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Figure  3:  Local  coordinates  associated  with  a  triangular  element 

element  into  three  triangles,  where  the  area  of  the  triangle  opposite  node  i  is  denoted 
by  4,-.  To  facilitate  the  integration  over  triangular  domains,  encountered  in  (15)  (46). 
we  introduce  for  each  triangle  a  ‘natural’  coordinate  system  (L,,  Z.;.  Lk)-  where  L, 
is  known  as  the  ‘area  coordinate’  or  ‘shape  function'  associated  with  node  i  of  the 
element  [11,  p.  llOj,  and  is  defined  as 

1=1 

where  A  is  the  triangle  area.  In  terms  of  the  area  coordinates,  the  local  position 
vector  g,  may  be  expressed  as 


Si  —  ^kLj  tjLk 


(48) 


To  represent  the  patch  and  aperture  currents  on  each  triangular  element,  we  introduce 
a  vector  basis  function  Ai,  given  as  (cf.  [12]) 


A,= 


2A 


(49) 


Noting  that  the  gradient  of  the  shape  function  is 


VtL, 


-n, 


2.4 


(50) 
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where  h,  is  a  unit  vector  [)erpen(licnlar  to  edge  ;  in  llu'  phue'  and  |)uin.iuig  out  ui  ’Sir 
triangle  element  {.set>  Fig.  3),  we  find  tin*  divergence  of  /I,  as 

V, -4.  =  ^  idh 

\V^  also  find  that  n,  •  yl,  is  constant  on  edge  t.  which  makes  it  easy  !o  enicjnc  'he 
continuity  of  the  normal  component  of  current  density  between  the  elentents  that 
share  this  edge.  The  patch  current  density  and  its  divergence  on  ea<  h  eh ‘meat  may 
now  be  approximated  as 

*  =  =  (Vi, 

,=1  .=  1 

where  /,  is  the  electric  current  leaving  the  element  through  edge  i.  Similarly,  the 
aperture  current  density  and  its  divergence  on  each  element  are  approximated  as 

Ms  =  Y1  VrMs  = 

1=1  t=i 

where  K,  is  the  magnetic  current  leaving  the  element  through  edge  i. 

3.3  Global  system  assembly  and  solution  procedure 

The  testing  functions  in  the  integral  equations  (45)-(46)  are  drawn  from  the  .same  sets 
as  the  basis  functions  used  to  represent  the  patch  and  aperture  currents.  Hence.  {/!„, } 
consists  of  with  i  =  I,  2,  and  3,  and  n  ~  I,  •  •  • ,  N,,  where  N,  is  the  total  number 
of  triangular  elements  of  Sp  and  S^-  In  the  above,  it  was  necessary  to  introduce  the 
super.script  n  to  distinguish  the  local  basis  functions  associated  wi»h  element  n.  In 
what  follows,  where  there  is  no  danger  of  confusion,  this  element  superscript  will  be 
omitted  for  notational  simplicity. 

When  the  expansions  (52) -(53)  are  substituted  into  the  integral  equations  (45) 
(46),  the  coefficients  and  are  constrained  by  the  boundary  conditions,  which 
require  the  continuitj  of  the  normal  components  of  J5  and  Ms  across  the  edges 
shared  by  adjacent  elements,  or  their  vanishing  at  the  boundary  edges  of  S'p  and  .h,,. 
respectively.  If  the  numbers  of  the  non-boundary  triangle  element  edges  on  the  patch 


i  U') 


aiul  the  aperture  ar«‘  (l('notetl  bv  .V,,  atul  A re>,|)<-etiv('iy.  atid  the  nuiu*  •  r  >.!  the 
unknown  current  expansion  coetfieieius  by  .V.  tiien  .V  =  A.,  r  A  ,,  As  a  le'-uh  of  ilu' 
procedure,  the  couphxi  integral  iMjualions  (  lul  t  lb»  .ire  <unveit(‘d  :n!<>  an  alaidiran 
system 


1 '  -  •  fin  j 

r  .  3 

[/•d 

’  f'-d  ' 

If)  1  f  V  ‘  ! 

[  ]  !.  '  'n  .'i  j 

.  i^'-d . 

.  [H,.j 

Here,  the  A -by- A’  system  matrix  consists  ot  tour  submatrices,  wiicre  is  the 

Ap-by-.V.,  global  impedance  matrix,  is  the  Ap-by  A,  ghdial  (oufiling  matrix. 

[Ann]  -  —  [Ain]^-  where  the  superscript  /  indicatc's  the  matrix  transpose,  and  I  ^ 
is  the  A'a-by- A'a  global  admittance  matrix.  I  he  A -by-1  vector  of  the  u'’kno.vn  current 
expansion  coetFicienls  in  (.u  l)  comprises  the  Ap-by-1  global  vector  j/,i;  of  (latch  current 
coefficients  and  the  A’,,-by-i  global  vector  [An]  of  aperture  current  coeilicienls.  I  he 
.V-by-1  global  excitation  vector  consists  of  the  Ap-by-1  global  patch  voltage  <‘xcitation 
vector  and  the  Ad-by-1  gi  bal  aperture  current  excitation  vector  [II., ij- 

To  assemble  the  global  system  (54),  w'e  consider  one  source'  element-test  element 
pair  at  a  time.  Each  such  pair  will  in  general  contribute  to  nine  elements  of  the 
global  system  matrix.  It  is  convenient  to  view  these  contributions  as  the  entries  ol 
a  local  3-by-3  system  matrix  corresponding  to  the  element  pair.  Furthermore,  each 
element  will  in  general  contribute  to  three  entries  of  the  global  excitation  vector, 
and  these  contributions  may  be  assembled  into  a  local  3-by-l  excitation  vector.  lo 
b"  more  specific,  let  the  global  indices  of  the  source  and  test  elements  be  n  and  n?, 
respectively.  Also,  suppose  that  both  elements  are  on  the  microstnp  patch.  Then, 
the  associated  local  system  will  take  the  form 


where  is  the  local  impedance  matrix.  [/,*”*)  is  the  local  vector  of  current  coeifi- 

cients,  and  [Vj"'*]  is  the  local  voltage  excitation  vector.  Here,  the  primed  local  source 
element  indices  i'  and  k'  follow  the  same  cyclic  convention  as  the  unprirned  ones. 
Observe  that  there  exists  a  unique  mapping  between  the  local  coefficients  {/,*'’’}  of 
each  element  and  the  global  current  coefficients  {/n}.  where  the  reierence  directions 
of  the  latter  are  specified  by  the  order  in  which  the  element  nudes  appear  in  the  input 
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geometry  data.  I  his  mappin"  <!eterniiues  lo  wliich  entries  of  tlie  glof»al  sysKnti  ! 's  i  i 
should  the  elements  of  (-55)  l)e  a<ided,  aiul  with  what  signs,  l.oi  .d  s',  stems  sitnilat  to 
(55)  arise  when  both  the  source  and  test  elements  are  in  the  aperture.  au<!  when  one 
of  them  is  on  the  patch  and  the  other  in  the  aperture.  Below,  we  dt's.  rihe  ihes*'  hjcal 
systems  in  more  detail. 

Consider  first  a  test  element,  5'“.  on  the  patch,  'riien,  if  the  sourci'  element, 
is  also  on  the  patch,  the  entries  of  the  resulting  local  impedance  matri.x  are  found  as 


If,  on  the  other  hand,  the  source  element,  .S’",  is  in  the  aperture,  there  results  a  local 
coupling  matrix,  with  the  entries  given  as 


rp(mn) 

I 


Oi 


Each  patch  element  also  contributes  a  local  voltage  excitation  vector,  whose  entries 
are 

158) 

Next,  consider  a  test  element,  5”,  in  the  aperture.  Then,  if  the  source  element. 
5p,  is  on  the  patch,  there  arises  a  local  coupling  matrix,  which  may  be  shown  (from 
reciprocity  considerations  [7,  pp.  116-120])  to  be  the  negative  transpose  of  the  cou¬ 
pling  matrix  given  by  (57).  If,  on  the  other  ha.nd,  the  source  element,  5,",  is  also  in 
the  aperture,  there  arises  a  local  admittance  matrix,  whose  entries  are  given  as 


m 

y.(mn)  _ 

(A,.(G5,A„);,)^_^-(i,(G'*, 

1)') 

ASJ  /gm 

(60) 

^(mn)  __ 

'ill 

(61) 

Finally,  each  aperture  element  also  contributes  a  local  current  excitation  vector,  with 
the  entries 

np'  =  (a.  ;  HP  -  (62) 
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The  integrals  over  source  coordinates  in  (.%)  and  (60)  (61)  invt)lve  kernels  that 
are  singular  when  the  test  and  source  elements  coincide.  These  singularities  are  ex 
tracted  and  integrated  analytically  [13|.  leaving  well-hehave<l  integrals  over  triangular 
elements,  which  are  numerically  evaluated  by  a  Cau.ssian  (juadrature  [11.  p.  ll.lj.  On 
the  other  hand,  the  testing  (exterior)  integrals  in  (56)  and  (60)  (61).  as  well  as  the 
integrals  in  the  excitation  terms  (58)  and  (62),  have  regular  and  slowly-varying  inte¬ 
grands.  and  may  thus  be  approximated  using  a  one-point  quadrature  rule  [12].  Tor 
example,  the  integral  in  (56)  is  approximated  as 

.  (m) 

(a, ;  (Gi  .  ^  ■  (Gi(r<"l  1 C) ,  AAp%, 


where  and  denote,  respectively,  the  local  (with  respect  to  node  z— -see  Fig.  3) 
and  global  position  vectors  of  the  centroid  of  element  5™.  This  one-point  approxima¬ 
tion  results  in  significant  savings  in  th<  computational  effort,  even  though  it  sacrifices 
the  symmetry  properties  of  the  impedance  and  admittance  matrices.  The  integrals 
appearing  in  (57)  have  regular,  but  rapidly  varying  kernels  (which  represent  electric 
field,  rather  than  potentials),  especially  when  the  substrate  between  the  microstrip 
patch  and  the  aperture  is  electrically  thin.  For  this  reason,  both  the  interior  and 
exterior  integrals  in  (57)  are  evaluated  by  Gaussian  quadratures  without  further  ap¬ 
proximations. 

Once  the  complex- valued  matrix  equation  (54)  is  assembled  and  solved,  which  is 
accomplished  by  standard  procedures  (LU  factorization  with  partial  pivoting,  followed 
by  a  forward  and  back  substitution  [11,  p.  120]),  the  current  density  within  each 
triangular  element  may  be  obtained  from  (52)  or  (53). 


3.4  Evaluation  of  spectral  integrals  and  Floquet  series 

For  the  solution  procedure  described  above  to  be  practical,  the  Sommerfeld-type 
spectral  integrals  that  occur  in  (17)-(21)  and  (24)-(25),  as  well  as  the  Floquet  series 
that  appear  in  (22)-(23)  and  (26),  must  be  efficiently  evaluated.  To  accomplish  that, 
the  Sommerfeld  integrals  are  accelerated  by  asymptotic  integrand  subtraction  and 
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the  method  ot  averages  [11].  In  adtiitioii  to  these  terhni<]ues,  an  iiiterpolaticai  and 
table  look-up  scheme  is  implemented  to  lurther  reduce  the  computation  time'  [lol 
The  integration  path  is  properly  deformed  to  avoid  the  integrand  singularities,  whit  h 
occur  on  or  near  the  real  axis  in  the  Ap  plane  [16], 

The  Floquet  series  are  accelerated  as  well,  by  a  combination  of  the  Kunmier 
and  Poisson  transformations  [17].  As  a  result  of  this  procedure,  the  original  slowly 
convergent  Floquet  series  is  converted  into  the  sum  of  an  accelerated  spectral  serit's 
and  an  exponentially  convergent  spatial  series.  Observe  that  in  (61)  these  series  are 
integrated  against  the  basis  functions  over  the  source  and  test  triangular  elements. 
The  source  element  integrals  are  introduced  inside  the  spectral  sums  and  evaluated 
analytically  [17,  pp.  107--1 14],  thus  further  accelerating  the  convergence  of  these  series. 
The  same  integrals  over  the  spatial  series  are  evaluated  by  a  Gaussian  quadrature,  as 
discussed  in  Sec.  3.3. 

3,5  Far  field  and  RCS  computation 

Once  the  coefficients  T  and  Ki  are  found  for  each  element,  the  current  expansions 
(52)-(53)  are  substituted  into  (42)-{43)  to  determine  the  far  zone  fields.  The  integrals 
encountered  in  (42)-(43)  are  then  recognized  as  Fourier  transforms  of  the  vector  basis 
functions  (49),  evaluated  at  kp  =  k°.  In  view  of  (49)  and  (48),  the  Fourier  transform 
of  Ai  associated  with  a  triangular  element  may  be  expressed  as 

A,  = 

=  (64) 

where  the  T,  denotes  the  Fourier  transformed  shape  function  Li.  Upon  using  the 
procedure  of  Appendix  B,  we  may  express  Li  as 

+  ^  -  iio(nA/2)l 
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where  p^,  is  a  position  vector  of  the  midpoint  of  edge  i  of  the  element  .  <i,  =  it  t,. 
6,  =  y  •  and  denotes  the  spherical  f3essel  function  of  order  it.  V\'e  note  that 
(65)  is  evaluated  for  kp  =  fc^,  as  is  required  in  (42)  (Id),  then  k,,  -  k^.  it  -  p  .  and 
V  =  It  can  be  shown  that  (65)  approaches  a  finite  limit  .4/4  as  X;,  —  0. 

Once  the  microstrip  patch  and  the  waveguide  aperture  contributions  to  tlie  fat- 
zone  fields  are  determined,  the  RCS  is  readily  obtained  from  (44). 

4  Experimental  Setup  and  Measurement  Proce¬ 
dures 

The  experiment  was  done  in  X-band,  rather  than  in  the  millimeterwave  range,  to 
reduce  the  effect  of  fabrication  tolerances  on  the  results.  The  components  of  the 
measured  structure  are  shown  in  Fig.  4.  Circular  and  rectangular  microstrip  patch 
antennas  were  investigated,  each  excited  through  a  waveguide- backed  concentric  rect¬ 
angular  slot,  as  illustrated  in  Figs.  5a  and  5b,  respectively.  The  substrate  and  super- 
strate  materials  used  in  the  experiment  (and  also  in  the  numerical  examples  presented 
in  Sec.  5)  are  non-magnetic  and  isotropic.  Therefore,  the  nth  layer  may  conveniently 
be  characterized  by  a  complex  number  c„(l  -  jtan(5„),  where  c„  is  the  real  relative 
dielectric  and  tan^n  is  the  loss  tangent.  The  substrate  used  in  antennas  of  Fig.  5  has 
Cl  =  2.2,  tan^i  =  0.001,  and  thickness  h  =  3.15  mm.  The  cover  layer,  if  present, 
is  made  of  the  same  material  with  thickness  t  —  1.57  mm.  The  microstrip  antennas 
were  mounted  on  the  waveguide  flange  using  plastic  screws  (see  Fig.  4). 

The  dominant  mode  reflection  coefficient  F,  referred  to  the  aperture  plane,  was 
measured  using  the  HP-8510B  network  analyzer.  Prior  to  the  measurement,  the 
thru-reflect-line  (TRL)  two-port  calibration  method  [18]  was  used  to  eliminate  the 
systematic  errors  due  to  the  coax-to-waveguide  adapter  and  to  establish  the  mea¬ 
surement  reference  plane  (MRP)  at  the  aperture  location.  In  order  to  implement  the 
TRL  calibration  procedure,  three  sets  of  measurements,  referred  to  as  thru,  reflect. 
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Figure  4:  Components  of  the  measured  structure. 
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Figure  5;  (a)  Circular  and  (b)  rectangular  microstrip  patch  antennas  used  in  the  mea¬ 
surements.  Each  antenna  is  fed  through  a  waveguide-backed  concentric  rectangular 
slot  in  the  ground  plane. 
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Figure  6:  (a)  Thru,  (b)  reflect,  and  (c)  line  (delay)  standards  used  in  the  HP-85 lOB 
TRL  calibration  procedure. 

and  line  (or  delay)  were  taken.  The  thru,  reflect,  and  line  calibration  standards  were 
fabricated  from  four  WG-90  waveguide  segments,  as  illustrated  in  Fig.  6.  The  thru 
was  arbitrarily  selected  to  be  84.2  mm  long  (the  length  of  the  thru  should  be  at  least 
two  guide  wavelengths,  to  reduce  the  interference  between  the  coax  adapters),  which 
resulted  in  the  reflect  length  of  42.1  mm.  A  waveguide  feed  of  the  same  length  was 
used  in  the  actual  measurements  (see  Fig.  4),  to  ensure  that  the  MRP  coincides  with 
tne  aperture  plane.  Because  no  two  coax  adapters  are  identical,  two  reflect  standards 
were  fabricated  and  used  in  the  HP-8510B  calibration,  as  indicated  in  Fig.  6b.  The 
line  (delay)  was  made  9.9  mm  longer  than  the  thru.  This  9.9  mm  length  difference 
between  the  thru  and  the  line,  which  at  the  center  frequency  =  10 GHz  is  approx¬ 
imately  a  quarter  of  the  guide  wavelength  Xg  =  39.7  mm,  results  in  a  time  delay  of 
24.94  ps,  and  this  value  was  keyed  into  the  HP-8510B  network  analyzer  during  the 
TRL  calibration  procedure. 

The  procedure  for  the  far  field  measurement  is  rather  standard,  and  is  not  de¬ 
scribed  here  to  conserve  space. 
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5  Sample  Computed  and  Measured  Results 

In  this  section  we  present  sample  computed  results  ol)taine<l  usinir  a  coinpiila-r  pro- 
gram  implementing  the  procedures  dt?veloped  in  Sec.  d.  For  sevcnai  antenna  Knifigu 
rations,  these  results  are  compared  with  the  corresponding  nu'asured  data  ubtained 
by  the  experimental  procedures  described  in  Sec.  1.  Although  no  publislied  results 
were  available  for  comparison  for  waveguide- backer!  microstrip  patch  antennas,  (wer  v 
effort  was  made  to  validate  at  least  the  most  important  modules  of  the  developed 
computer  code  against  independently  obtained  data.  First,  the  part  of  the  program 
that  deals  with  an  arbitrarily  shaped  microstrip  patch  on  a  grounded  substrate  (in 
the  absence  of  the  aperture)  was  extracted  and  used  to  compute  the  scattering  rhar- 
acteristics  for  rectangular  patch  antennas  on  both  isotropic  and  uniaxial  substrates. 
The  computed  RCS  results  were  found  to  closely  agree  with  tho.se  published  by  New- 
mann  and  Forrai  [19]  and  Pozar  [20].  Second,  the  part  of  the  code  that  deals  with 
the  aperture  was  isolated  and  used  to  compute  the  RCS  of  a  narrow  slot  in  a  ground 
plane.  This  problem  was  then  related  via  the  B  binet’s  principle  [21,  p.  500]  to  that 
of  a  thin  wire  scatterer  with  equivalent  radius  [22],  and  the  latter  was  analyze  u.sing 
a  commercial  code  PCAAD  [23].  Again,  close  agreement  between  the  corresponding 
RCS  results  was  observed.  Third,  a  stripped  down  version  of  our  program,  which  did 
not  include  the  microstrip  patch  part,  was  used  to  analyzed  a  rectangular  waveguide 
radiating  through  a  centered  rectangular  slot  into  a  half-space.  The  computed  equiv¬ 
alent  magnetic  current  in  the  aperture  was  found  to  closely  agree  (both  in  magnitude 
and  phase)  with  the  corresponding  result  obtained  by  Harrington  and  Mautz  [24. 
Fig.  15a].  Finally,  a  rectangular  waveguide  radiating  through  a  centered  rectangular 
aperture  covered  by  a  dielectric  layer  was  analyzed.  For  this  problem,  in  Fig.  7  we 
compare  our  aperture  admittance  results  with  the  computed  and  measured  data  ob¬ 
tained  by  Bodnar  and  Paris  [25].  Although  the  agreement  between  the  three  sets  of 
results  is  judged  to  be  good,  we  note  that  our  data  are  closer  to  the  measured  results 
than  the  data  computed  by  Bodnar  and  Paris. 

VVe  next  pre.sent  measured  and  computed  results  for  four  waveguide-backed  mi¬ 
crostrip  patch  antennas  with  isotropic  substrates,  excited  through  a  centered  rect- 
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Figure  7:  Aperture  conductance  G  and  susceptance  B  normalized  to  the  dominant 
mode  wave  admittance  of  a  flanged  rectangular  waveguide  radiating  through  a  cen¬ 
tered  rectangular  aperture  covered  by  a  slab  with  dielectric  constant  ci  =  2.25 
and  thickness  h  =  3.201mm.  The  waveguide  has  dimensions  a  =  22.86  mm  and 
b=  10.16  mm,  and  the  aperture  size  is  0.7a-by-0.86. 
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Figure  8;  Plots  of  (a)  VSWR  and  (b)  phase  of  F  for  a  waveguide-excited  circular 
patch  antenna  without  a  superstrate  (Antenna  1). 

angular  slot  by  the  dominant  (TEio)  mode.  The  parameters  of  these  antennas  are 
collected  for  ezisy  reference  in  Table  1.  Observe  that  Antennas  1  and  3  are  circular 
(see  Fig.  5a),  while  Antennas  2  and  4  are  rectangular  (see  Fig.  5b).  a\1so.  .Antennas  1 
and  2  do  not  have  a  superstrate,  whereas  Antennas  3  and  4  are  covered  by  a  dielec¬ 
tric  slab  made  of  the  same  material  as  the  substrate.  In  all  cases  the  same  X-band 
rectangular  waveguide  (WG-90)  is  used,  with  the  interior  dimensions  a  =  22.86  mm 
and  b=  10.16  mm.  Referring  to  Table  1,  observe  that  the  largest  discrepancy  between 
the  computed  and  measured  values  of  the  resonant  frequency,  where  the  minimum  of 
the  voltage  standing  wave  ratio  (VSWR)  occurs,  is  1%.  Note  also  that  the  resonant 
values  of  VSWR  are  close  to  one  and,  therefore,  the  antennas  are  nearly  matched  at 
their  resonant  frequencies.  Their  bandwidths  (BW  in  Table  1),  however,  are  narrow 
(which  is  characteristic  of  microstrip  antennas)  and  do  not  exceed  5%. 

The  computed  and  measured  VSWR  and  phase  of  the  dominant  mode  reflection 
coefficient  (ZF)  at  the  aperture  for  the  four  antennas  are  shown  in  Figs.  8-11.  In 


24 


Table  1:  Parameters  of  sample  waveguide-fed  microstrip  pah  h  aiilenuas. 


Atiteima 

- — _ L 

‘ _ L 

3  4 

Relative  substrate  t 

2.2 

2.2 

— 

2:1 

2.2  i 

Substrate 

thickness  h  (mm) 

:}.15 

3.15 

3.15 

3.15  1 

Relative  superstrate  t 

N’/A 

.\7A 

2.2 

2.2 

Superstrate 
thickness  t  (mm) 

N/A 

N/A 

1.57 

1.57 

tan  6 

0.001 

0.001 

0.001 

O.OOl 

Patch 

dimensions  (mm) 

Dp  ^  14.0 

Ip  =  14.0 
Wp  =  11.7 

Dp  =  13.7 

Lp  =13.7 
Bp  =  11.4 

Aperture 
dimensions  (mm) 

Wa  =  H.O 

La  =  0.3 

Wa  =  14.0 

La  =  0.3 

tv;  =  13.7 

La  =  0.3 

B;  =  13.7 

La  =  0.3 

Measured  resonant 

frequency  (GHz) 

10.09 

10.07 

9.95 

9.90 

Computed  resonant 
frequency  (GHz) 

10.17 

10.17 

10.0 

10.0 

Error  in  resonant 

frequency  (%) 

0.79 

0.99 

0.5 

1.0 

Measured  gain  (dB) 

8.63 

8.8 

7.9 

8.2 

Measured 

lowest  VSWR 

1.22 

1.137 

1.114 

1.013 

Calculated 

lowest  VSWR 

1.23 

1.151 

1.075 

1.012 

Measured  BW  (%) 

4.0 

4.24 

4.74 

4.92 

Computed  BW  (%) 

3.8 

4.03 

4.4 

4.52 

Figure  9:  Plots  of  (a)  VSWR  and  (b)  phase  of  F  for  a  waveguide-excited  rectangular 
patch  antenna  without  a  superstrate  (Antenna  2). 


(a)  (b) 

Figure  11:  Plots  of  (a)  YSWR  and  (b)  phase  of  F  for  a  waveguide-excited  rectangular 
patch  antenna  with  a  cover  layer  (Antenna  4). 
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Figure  12:  Resonant  currents  of  Antenna  1.  (a)  Patch  electric  current  (imaginary 
part),  (b)  Slot  magnetic  current  (real  part  of  the  longitudinal  component). 


Figs.  12-13  we  show  plots  of  the  computed  patch  electric  current  and  the  slot  magnetic 
current  for  Antennas  1  and  2  at  their  resonant  frequencies.  Since  these  currents  are  in 
general  complex-valued,  only  the  dominant  part  (real  or  imaginary,  as  the  case  may 
be)  is  shown  in  each  case.  These  figures  also  illustrate  the  triangular  mesh  models 
used  in  the  analysis.  The  circular  patch  of  Antenna  I  and  the  rectangular  patch  of 
Antenna  2  were  approximated  by  394  and  288  triangular  elements,  respectively.  In 
both  cases  the  slot  was  modeled  by  44  elements.  For  Antenna  2,  this  resulted  in 
a  451-by-451  global  system  matrix  in  (54),  and  a  computation  time  of  15  minutes 
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Figure  13:  Resonant  currents  for  Antenna  2.  (a)  Patch  electric  current  (imaginary 
part),  (b)  Slot  magnetic  current  (real  part  of  the  longitudinal  component). 
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per  frequency  point  on  a  12-MlPS  computer.  In  Figs.  11  I.').  we  show  the  measured 
and  computed  far  field  patterns  for  the  rectangular  .Vntennas  2  and  1.  1  he  .scallu[) 
observed  in  the  measured  E-plane  patterns  is  almost  certainly  caused  by  the  surla<  <* 
wave  diffraction  at  the  edges  of  the  finite-size  ground  plane  [26),  This  effect  is  more- 
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Figure  14:  Measured  and  computed  far  field  patterns  for  Antenna  2. 
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Figure  15:  Measured  and  computed  far  field  patterns  for  Antenna  4. 

pronounced  for  Antenna  4,  because  the  superstrate  increases  the  intensity  of  the 
excited  surface  wave. 
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Finally.  we  present  sample  results  for  Antenna  1  under  plain'  wavt'  I'xrital  ion. 
In  fig.  16  we  plot  vs.  frequency  the  magnitude  of  llu*  tlominant  nnxh'  aperture 
transmission  coefficient  T,  where  the  latter  is  defined  in  (11).  .\s  e.xpected.  the  peak 
of  transmission  occurs  at  the  resonant  frequency  of  the  waveguide-ilriven  antenna 
(see  Table  1).  In  Fig.  17  we  plot  vs.  frequency  the  monostatic  RCS  refcured  to  1  nr‘ 


Figure  16:  Magnitude  of  the  dominant  mode  aperture  transmission  coefficient  vs. 
frequency  for  Antenna  2. 

(the  units  are  dBsm — 'decibels  above  a  square  meter’  [27,  p.  160])  for  a  polarized 
incident  plane  wave  with  —  I  V/m.  In  addition  to  the  total  RCS  of  the  antenna, 
we  also  plot  the  contributions  from  the  microstrip  patch  and  the  slot.  We  note  that 
the  first  RCS  peak  is  clearly  due  to  the  first  resonant  mode  of  the  patch.  At  the  second 
peak,  which  occurs  near  the  resonant  frequency  of  the  waveguide-driven  antenna,  the 
slot  radiation  is  the  dominant  effect.  We  also  note  that  there  is  an  RCS  minimum 
between  10  and  12  GHz,  caused  by  a  destructive  interference  of  the  patch  and  slot 
radiated  fields. 
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Figure  17:  Monostatic  RCS  of  Antenna  2. 
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6  Summary  and  Conclusions 


We  have  presented  a  rigorous  integral  ecpiation  analysis  ol  a  micrustrip  pat<  h  anlenna 
excited  through  an  aperture  in  the  mouth  of  a  rectangular  waveguide,  ur  illuniinal<‘<l 
by  a  plane  wave.  The  substrate  and  superstrate  may  comprise  any  numix'r  of  isotropic 
or  uniaxial  material  layers.  The  patch  and  the  aperture  may  both  Ix'  of  arbitrary 
shape.  Hence,  waveguide-fed  microstrip  patch  antennas  of  various,  possibly  irregular 
shapes,  residing  in  multilayer,  possibly  uniaxial  media,  may  be  analyzed  within  a  sin 
gle  formulation,  using  the  same  computer  program.  The  analysis  has  been  validated 
against  experimental  and  published  data.  We  have  also  shown  that  a  good  impedance 
match  may  be  achieved  in  this  antenna  configuration,  but  only  in  a  narrow  frerpiency 
band. 


Appendix  A 

Transmission-Line  Analog  of  Layered  Medium 

In  deriving  the  integral  equations  of  Sec.  2.1,  we  have  employed  a  transmission-line 
network  analog  of  the  layered  medium,  in  which  each  layer  is  represented  by  a  trans¬ 
mission  line  section,  as  illustrated  in  Fig.  18.  This  analog  comprises  two  networks. 
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Figure  18;  Layered  dielectric  medium  and  its  transmission-line  network  analogue. 
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which  arise  from  th*'  decomposition  of  the  elect romagiu'tic  Held  into  tw(j  [jartial  held- 
that  are  transverse-magnetic  (  IM)  and  transverse-electric  (  I  K)  to  c.  1  he  quant iiii-s 
corresponding  to  these  networks  are  flistinguished  l)v  the  superscripts  <  and  h.  re¬ 
spectively,  I  he  characteristic  impedance  and  propagation  constant  of  the  nth  -ectuui 
of  the  riv  and  I  M  transmission  lines  are  given  as 


Let  the  network  be  excited  by  a  lA  shunt  current  source  located  at  in  the  nth 
line  section  of  length  t/^.  as  illustrated  in  Fig.  19a.  Lhen.  the  voltage  V[{z  1  :')  and 


•  cLn  * 


(b) 


Figure  19:  Typical  transmission  line  section  with  (a)  1 A  shunt  current  source  and  (b) 
IV  series  voltage  source. 
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current  /,(c| 


at  a  point  r  witliin  this  line  se<-f i<jn  ohcy  the  ('(juatimis  [2ts.  p.  .  h 


^i;{c|c')  =  -jk,n  ld:\:') 

dz 

4-i.{  =  \=')  =  -2A-,„>;v;(r|c')  +  ,'U-c') 

(ir 


( f  iN  ! 


where  6  denotes  the  delta  function.  Here  and  below,  we  omit  the  the  sup<'rs<  ript  p 
for  notational  simplicity.  From  (68)  (69).  V'(;|;')  may  be  expressed  in  the  traveling 
wave  form 


~=r 


2[l-  {\{z,)  [\(r,)] 

•  [l+  T„(~'o)6-^'*-<-’<-->1  11+ 


ro) 


where  2<  =  min(^,c'),  =  max(z,c'),  and  Fnf^o)  and  rn(co)  are  the  reflection 

coefficients  'looking  to  the  left’  and  ‘looking  to  the  right.’  respectively,  at  any  location 
Zo  within  the  line  section.  By  means  of  the  translation  formula 


(71) 


where  the  right  and  left  arrows  correspond  to  the  upper  and  lower  signs,  respectively, 
these  reflection  coefficients  can  always  be  expressed  in  terms  of  the  terminal  reflection 
coefficients,  r„  and  Fn,  which  are  related  to  the  corresponding  terminal  impedances 
Zn  and  n  (  see  Fig.  19)  as 


Zn  +  Zn 


The  voltage  Vv{z  \  z')  and  current  ly{z\z'),  excited  by  a  IV  series  voltage  source  v 
in  the  nth  line  section  (see  Fig.  I9b),  satisfy  equations  dual  to  (68)-(69),  which  are 
obtained  from  the  latter  by  making  the  substitutions:  Vi  — +  /„,  /,  — >  Zn  — >■  Vy,, 

and  Yn  —*  Zn-  Furthermore,  it  can  be  shown  that  the  following  symmetry  and 
reciprocity  relations  hold  (cf.  [28,  p.  194]) 


\\{z\z')^Y{z'\zY  ln{z\z')  =  In{z'\zl  Vn(z  \  z')  -  h{z' \  z) 


78) 
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As  a  result,  |  c')  ran  l)e  ohtaiiie<l  from  (70)  by  replarins  i"  bit’*'!'  ^ b\-  7f, 
(which  also  causes  the  reHectioii  coefficients  to  cliange  signs).  [,{z  [  ;')  follows  from 
(68)  and  (70).  and  ll,(-!-')  ‘Uay  then  be  obtained  from  the  last  relation  in  (73). 

The  voltage  and  current  on  the  nth  transrni.ssion  line  section  that  is  source-fre'e 
satisfy  the  homogeneous  form  of  (68)-(69).  Front  these  equations,  the  voltage  at  any 
point  z'  within  the  same  line  section  may  be  found  as 

S’©}  r  pH; 

V'(.')  =  V(.,)— ^ -  1+  l\iz, 

1+  r„(t:(j) 

where  the  upper  (lower)  sign  corresponds  to  r  >  :r„  (^  <  Zo).  In  the  above,  z.j  —  or 
2o  =  depending  on  whether  the  source  is  located  to  left  or  right,  respectively,  of 
the  nth  line  section,  and  V(Zo)  is  is  the  voltage  across  the  line  terminals  at  r  —  Co-  VVe 
have  omitted  the  subscript  of  V  in  (74),  because  the  latter  applies  irrespective  of  the 
nature  of  the  source,  and  dropped  z'  from  its  argument,  because  (74)  only  implicitly 
depends  on  the  source  location,  which  is  outside  the  line  section.  I'he  current  l{z) 
corresponding  to  (74)  may  be  obtained  by  substituting  the  latter  into  (68). 

Appendix  B 

Fourier  Transform  of  a  Triangle  Shape  Function 

The  Fourier  transform  of  a  shape  function  L,  associated  with  the  nth  triangular 
element  is  given  as 

L,=  I  (75) 

JSn 

where  kp  =  ukp.  To  evaluate  the  integral  in  (75),  it  is  helpful  to  first  convert  it  to 
a  line  integral  around  the  boundary  contour,  5S„,  of  Sn  (cf.  [29],  [30].  [31]).  This  is 
most  easily  accomplished  by  noting  that 
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an<l  by  making  use  of  the  divergence  fln'orem.  As  a  result,  wIkmi  L],  /  0.  we  (obtain 

^-'1  ~  •  +  V  ■  ■  hdf  (77 1 

2Ak^  Jos„ 

where  i)  =  z  x  u  and  h  denotes  a  unit  vector  normal  to  OS^  at  <  in  the  plane  and 
pointing  out  of  Sn-  The  integral  in  (77)  is  easily  evaluated  in  the  lo<-al  coordinati-ds. 
if  one  notes  that  when  (  is  on  edge  i,  h  =  n,,  0  <  <  1.  L,  —  0.  =  1  -  and 

d(  =  ^idLk  (see  Fig.  3).  Also,  in  that  case  p  —  +  g^,  where  g^  is  a  vector  from 
node  j  to  the  point  i  on  edge  i  of  the  element.  The  resulting  close-forn\  expression 
for  Z,  is  given  in  (65). 
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Analysis  of  Multiconductor  Transmission  Lines  of  Arbitrary 
Cross  Section  in  Multilayered  Uniaxial  Media* 

by 

(  'huna,  I  (i.  H.sir,  Roat-r  I' .  llarriiiatoir .  Kr/ys/ioi  A.  NRchai'.kiy  and  DaJiaii  Zhc!!*’ 

Abstract  A  niixed-pott'iitial  {*l<*rtri<-  ticdd  intcarai  aquatiiai  is  and 

applied  in  conjunction  with  the  method  of  moments  to  analyze  a  t ransmissi<ni  inn' 
system  consisting  of  multiple  conducting  strips  ot  arbitrary  cr<;ss  section  erniiedded 
in  a  stratified  medium  with  or  without  top  an<i/or  l)ottom  ground  planes,  Kacii  ia>er 
of  the  medium  is  possibly  uniaxially  anisotroj)ic.  with  its  ojjiica!  axis  irerpendini 
lar  to  the  dielectric  interfaces.  (’ompute<l  liispersion  curves  and  modal  currents  are 
presented  and.  when  possible,  are  compared  with  <iata  available  in  the  literature. 

1  Introduction 

Recent  advances  in  integrated  circuit  technology  have  made  microstrips,  siriplincs, 
coplanar  strips,  and  similar  wave-guiding  structures  attractive  not  only  in  microwave 
and  rnillimeter-wave  applications,  but  also  in  high-speed  digital  computers,  l  l.e  con¬ 
ductors  used  as  interconnects  between  \  LSI  devices  may  be  very  close  to  one  another, 
which  necessitates  treating  them  as  a  single  transmission  line  capable  of  ■'Ui»p<>rtinc 
several  modes,  rather  than  .several  isolated  transmission  lines.  The  interconnects  in 
modern  microwave  and  millimeter-wave  integrated  circuits  tend  to  have  trapezoidal 
cross  sections  due  to  etching  undercuts  or  as  a  result  of  the  epitaxial  growth  process 
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[1],  [2],  and  cannol  always  be  considered  iidinitely  tliin.  Tiiest*  interconnects  arc'  sup¬ 
ported  by  a  dielectric  substrate,  which  often  exhibits  uniaxial  anisotropy,  introducc'd 
in  the  manufacturing  process  [3j. 

Many  numerical  procedures  have  been  successfully  applied  in  the  past  twcj  decades 
to  obtain  frequency-dependent  characteristics  of  rnicrostrips  and  striplinc's.  but  most 
of  them  are  only  applicable  to  (or  optimized  for)  planar  conducting  strifes  cef  zent 
thickness  (cf.  [4],  [o],  [6],  [7j,  [8],  [9|.  [10],  [11],  [12).  [13].  [14],  [15].  [Kij.  [17],  tee  name 
just  a  few).  Relatively  few  papers  have  considered  laterally  open  microst  rij)  struct  ures 
with  conductors  of  other  cross  section  shapes,  such  as  rectangular  [ISj.  trapezoidal 
[19],  [20],  circular  [21],  [22],  or  rectangular  with  semi-circular  edges  [23].  The  concept 
of  equivalent  width  has  often  been  employed  to  approximately  take  into  account  the 
strip  thickness  [24].  However,  it  has  recently  been  demonstrated,  using  a  rigorous 
n'ixed-potential  integral  equation  (.MPIE)  approach  [19],  that  the  dispersion  curve 
for  a  finite-thickness  microstrip  lies  below  that  of  a  rnicroslrip  with  zero  thickness, 
which  is  opposite  to  what  is  observed  when  the  concept  of  equivalent  width  is  used. 

In  this  paper,  we  use  an  MPIE  approach,  which  was  originally  developed  for 
objects  in  isotropic  media  [25],  [19],  [26].  and  recently  extended  u.  .jl)jects  in  uniaxial 
media  [27],  to  analyze  a  transmission-line  system  compo.sed  of  multiple  conductors 
of  finite  thickness  and  arbitrary  cross  section,  embedded  in  a  medium  consisting 
of  an  arbitrary  number  of  planar,  possibly  uniaxially  anisotropic,  dielectric  layers. 
Computed  dispersion  curves  and  modal  currents  for  bound  modes  are  presented  and. 
when  possible,  are  compared  with  data  available  in  the  literatire. 

2  Formulation 

The  cross-sectional  view  of  the  structure  under  consideration  is  shown  in  Fig.  !. 
The  medium  consists  of  N  planar,  homogeneous  dielectric  layers,  with  the  interfaces 
parallel  to  the  xy  plane.  Each  layer,  say  the  nth,  is  characterized  by  permeability 
and  by  transverse  and  longitudinal  permittivities  and  re.'^pectively.  all  relative 
to  free  space.  The  top  layer  of  the  medium  may  extend  to  -f  oc^  along  the  ;  axis,  or  1><’ 
shielded  by  a  ground  plane  made  of  a  perfect  electric  conductor  (PEC).  Similarly,  the 


2 


A^rl) 


d, 

(optional  ground  plane) 


Figure  1:  Cross-sectional  view  of  a  multiconductor  transmission  line  embedded  in  a 
stratified  uniaxial  medium. 
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bottom  layer  may  extend  to  — oc  along  the  c  axis,  or  be  shielded  by  a  Pl-,('  ground 
plane.  There  are  A;,  PEC  strips  embedded  in  the  layered  tnediuiu,  all  uniform  and 
of  infinite  extent  along  the  //  axis,  but  of  arbitrary  cross  section  shap*'.  An  f lime 
dependence  is  assumed  and  suppresse<l  throughout. 

Since  we  are  interested  in  modes  propagating  in  the  //  direction,  we  may  assume 
that  the  phase  factor  is  common  to  all  the  fields  and  currents,  where  J  is  the 

propagation  constant  to  be  determined.  Hence,  we  may  express  the  surface  current 
density  as 

J{r)^  (j) 

where  (:  is  the  arc-length  coordinate  on  the  contours  of  the  conductor  cross  sections. 
By  enforcing  the  condition  that  the  tangential  electric  field  must  vanish  at  the  surface 
of  the  conductors,  we  obtain  an  electric  field  integral  equation  (  EFIE)  of  the  form 

.Vc  . 

«nxj]/  Q^-^{x.z\x\z')-J{e)df'  =  0,  rein.  n  =  1.2.  •■■..V.  (2) 

where  Qf''^{x,  z\x'z'),  with  x  =  x{i)  and  c  =  ~(f),  is  the  electric  field  dyadic  Green's 
function  of  the  layered  medium  [27],  and  where  iin  denotes  an  outward  unit  vector 
normal  to  the  boundary  I„  of  the  nth  conductor.  In  the  above  and  throughout, 
primed  quantities  denote  source  coordinates,  unit  vectors  are  distinguished  by  carets, 
and  dyadics  by  double  underlines. 

The  severe  source-region  singularity  of  the  kernel  of  the  EFIE  (2)  makes  it  un¬ 
suitable  for  a  direct  application  of  the  method  of  moments  [28].  [29].  Hence,  we  first 
transform  it  into  the  MPIE  form, 

itn  X  =  0,  r  G  L„.  n  =  1.2.---..V-  {:{) 

i=i 

where  Vj?  is  the  transverse  (to  y)  part  of  the  operator  nabla.  and  wl 

A,{x,z)=  f  (T) 

Jl,  — 

and 

^.{x,.G=  /  K‘^ix.=\x\z'){V',-yj3)-J{ndf'  (.3) 

J  Li 
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are  the  magnetic  vector  potential  and  the  electric  scalar  pot('ntial.  respc'c  lively,  due 
to  the  surface  current  on  the  dh  conductor.  These  potentials  are  not  unique,  as 
discussed  in  [30],  [26],  [27],  In  the  latter  reference,  two  different  .\1PIE  formulations, 
referred  to  as  the  "traditional''  and  the  "alternative."  are  developed  for  arbitrarily 
shaped  conductors  in  layered  unia.xial  media. 

The  e.xpressions  for  the  dyadic  kernel  and  the  scalar  kernel  l\  -’  comprise 
improper  spectral  integrals  of  the  form 


s^uik,)]  =  I  r  Rk,)  " 

TT  Jo  SI 


COS  kj;{x  —  x') 
sin  kj.{x  —  s') 


A-l'  dk, 


where  kj:  is  the  Fourier  transform  domain  counterpart  of  .r.  and  where  the  subscripts 
c  and  s  are  associated  with  the  cosine  and  sine  functions,  respectively,  and  n  assume.^ 
the  values  0  or  1.  Using  this  notation,  the  nonzero  elements  of  and  for  the 
traditional  MPIE  formulation  [31].  [27]  can  be  expressed  as 


Ay^Cx.-Sx'..-')  =  -  = 

Ay,(i,2|x'.x')  =  [/.yj-'i--')  - 

Ay,(x..-ix',--')  =  -  [i-'.yj--!--')  - 


fS) 

P)) 

(10) 

(11) 
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where  ijo  and  Aq  denote,  respectively,  the  intrinsic  inipedance  and  w'av<>munhtT  of 
free  space,  and  kp  =  (cl  +  iP.  The  subscripts  in  and  n  in  the  above  u  <ii(ate  that  the 
observation  point  (x,  r)  is  in  the  mth  layer,  and  the  source  point  f.r'. in  the  /;th 
layer.  In  deriving  (7)-(l3),  use  has  been  made  of  the  transmission-line  network  ana¬ 
log  of  the  layered  medium  [32.  Ch.  2],  which  is  illustrated  in  Fig.  2  for  a  three-layer 
geometry.  This  network  actually  represents  two  networks  (having  identical  configura¬ 
tions.  but  in  general  different  propagation  constants  and  characteristic  impedance.N) 
that  arise  from  the  decomposition  of  the  electromagnetic  field  into  partial  fields  that 
are  transverse-magnetic  (TM)  and  transverse-electric  (TE)  to  z  [32],  [33].  The  su¬ 
perscript  p  in  Fig.  2  stands  for  e  or  /?,  which  designate,  respectively,  the  rpianlities 
associated  with  the  TM  and  TE  networks.  The  propagation  constants  of  the  nth 
transmission  line  section  are  found  as 


~  \  ^Q^tnP-rn  ~  •  k^^  —  yk^S^nJhn  k 


14) 


where  the  branch  of  the  square  root  is  determined  by  the  condition  that  }  <  0. 

The  corresponding  characteristic  impedances  (and  admittances)  are  given  as 


ye 

"  V®  ko£tn 


yh.  _  ^  _  kQTJoPrr, 


(15) 


In  (7)-(13),  and  Il^^{z\z')  denote,  respectively,  the  voltage  and  current  at 

2  on  the  mth  transmission  line  section,  due  to  a  1  A  current  source  at  z'  on  the  nth 
line  section.  Similarly,  denote,  respectively,  the  voltage 

and  current  at  2  on  the  mth  transmission  line  section,  due  to  a  1  V’  voltage  source  at 
z'  on  the  nth  section.  These  transmission-line  Green’s  functions  are  deriverl  in  the 
appendix  for  a  medium  with  an  arbitrary  number  of  layers. 


3  Numerical  Method 

In  this  section,  the  method  of  moments  [28],  [29]  is  employed  to  solve  the  .MPIE  (3) 
for  the  multiconductor  transmission-line  problem  of  Fig.  1.  .As  the  first  step  of  the 
numerical  procedure,  we  approximate  the  cross  section  contours  of  the  conductors  by 
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Figure  2;  Transmission-line  analogue  of  a  layered  medium. 
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piecewise  linear  segments,  as  illustrated  in  Fig.  3.  The  arc-length  c,K>r(linate  i  will 
now  be  associated  with  the  approximated  contours,  instead  of  tlie  original  ones.  The 
method  of  moments  requires  that  the  unknown  currents  be  expaiKhnl  in  tertns  <of  a 
set  of  known  basis  functions  with  unknown  coefficients,  viz: 


./r(0  =  ECT('').  = 


where 


A^{() 


f  —  (f‘ 

/p  -  ( 

fP  -1  +  ^ _ _  fP  <-  (  c  fp 

j-sf/P  _/-P'  ‘  <  'j+l 

0 .  otherwise 


/P  ^  _Z±]_Ll  \p  —  1 -fp 

ei- 


7  +  1  ^7 


n;(o  = 


i,  (U<t<p, 

0 .  otherwise 


In  the  above  expressions,  the  subscript  j  and  the  superscript  p  signify  quantities 
related  to  the  jth  expansion  function  or  segment  on  thepth  conductor  (the  superscript 
p  is  omitted  in  Fig.  3  for  simplicity). 

The  MPIE  (3)  is  next  tested  with  yllf  and  /If.  In  this  process,  the  transverse 
nabla  operator  in  (3)  is  transferred  to  operate  on  the  testing  function  by  applying 
Green’s  first  identity,  viz: 


J  AUi)-'^oMx.z)de  =  -J[V,-A’:ini^ix,z]di 

=  _  [  [HM  _  ax 

J  Af  ^ 


X,  z)d[ 


To  save  computer  time,  the  following  approximations  are  used  in  the  testing  proce¬ 


dure: 


J n^{i)mde  « 
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where  /(<)  represents  the  scalar  potential  <l>(j-.r).  or  a  component  of  the  vector  pce 
tential  A(.r.r).  and  where  ^  is  the  arc-length  coordinate  of  a  point  specified  by 
the  position  vector  (cf.  [•'ig.  3) 


I 


4- 


T,+r 


14-1 


Furthermore,  when  computing  the  magnetic  vector  potential  due  to  we  approxi 
mate  the  resulting  integral  as 

1 


J  K^(.r.  :\s'.  z')  ■ 


dC 
(24) 

As  a  final  step,  we  substitute  the  expansions  (16)  into  the  tested  form  of  (3)  to 
convert  it  into  a  homogeneous  matrix  equation  for  the  current  expansion  coefficients. 
.Assuming,  for  simplicity,  that  there  are  only  two  conductors,  this  equation  has  the 
form 


;2oi 


2IU. 

^12, n 
,■^0 

yl2,ty 

■[or 

yll.yt 

yU.yy 

7l2,y^ 

'  7^2. yy 

[^] 

[27'1 

[0] 

721, yf 

'  y2'i,yy 
^ij 

^2^,yt 

^22, yy 

. 

“ 

.  (0). 

with  the  matrix  elements  given  as 
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2.7'"  =  f 


-^j+1 


^14-1,44-1 


) 


(tf’  -  fit  ,d 


27'"=  Y  -JD  (wfl’- ^07+, 


(26) 

(27) 

(28) 
(29) 


where  the  index  j,  when  used  as  a  subscript,  should  not  be  confused  with  the  imagi¬ 
nary  unit.  In  (26)-(29),  we  have  introduced  the  notation 
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z'yc 


!:il  j 


where  [t.t])  =  {t.y).  and  i  =  i{()  —  ^  for  (^  <  (  <  In  the  general  case  witii 

J+  2  ^  ^  ^ 

.\c  conductors,  the  matrix  in  (25)  comprises  4.\"  blocks  of  submatrices. 

The  kernel  functions  in  the  above  have  been  defined  by  the  spectral  integrals 
(7)-(13),  which  are  evaluated  by  a  composite  Gauss  quadrature.  lo  accelerate  the 
convergence  of  these  integrals,  we  first  subtract  from  the  integrands  their  large  ar¬ 
gument  asymptotic  forms.  Furthermore,  when  |t  —  .r'i  is  larger  ;  han  —  r'j.  we 
also  employ  the  method  of  averages  [25],  [34].  Finally,  the  closed-form  integrals  of 
the  asymptotic  forms  are  added  back  to  compensate  for  the  subtracted  terms.  The 
former  explicitly  exhibit  their  source-region  logarithmic  singularities,  which  are  inte- 
grable  and  are  easily  taken  care  of  in  evaluating  the  kernel  integrals  (30)  -(31). 

The  equation  (25)  has  nontrivial  solutions  only  for  those  values  of  d.  which  make 
the  matrix  determinant  vanish.  These  values  are  found  by  the  Muller  method  [35]. 
and  the  corresponding  modal  current  coefficients  are  then  determined  from  (25). 


4  Numerical  results 

In  this  section,  we  present  sample  numerical  results  for  the  propagation  constants  and 
modal  current  distributions  for  a  variety  of  transmission-line  configurations.  In  all 
examples  considered,  the  media  are  assumed  lossless  and  nonmagnetic  (i.e..  fJm  —  1 
for  all  layers).  Some  of  the  structures  analyzed  comprise  both  uniaxial  and  isotropic 
dielectric  layers.  If  a  layer  is  isotropic,  its  relative  permittivity  is  denoted  by  cr-  Only 
the  proper,  bound  modes,  which  propagate  unattenuated  with  a  real  propagation 
constant  (3,  are  considered.  The  dispersion  curves  are  given  either  for  T/Aq.  or  for  the 
effective  dielectric  constant  eeff  = 

In  Fig.  4,  we  present  dispersion  curves  for  a  circular-wire  transmission  line  embed¬ 
ded  in  a  grounded  two-layer  isotropic  medium  with  or  without  a  top  ground  olane. 
The  latter  configuration  was  first  analyzed  by  Fache  and  De  Zutter  [21],  using  an 
approach  especially  developed  for  v/ire  conductors,  and  their  results  are  shown  by 
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square  symbols  in  Fig.  4.  In  this  figure,  c^fr  is  plotted  \('rsus  the  elet  trieal  thi<  kne>'- 
of  the  substrate.  d/Xo.  where  Ao  is  the  iVee-space  wavelength.  I  he  wire  is  eonipleiely 
embedded  in  the  dielectric  slab  for  h/d  =  O.T-o  or  ().•’).  and  in  tin*  air  r<‘gioii  when 
h/d  —  1.25  or  1.5.  In  the  analysis,  the  circular  cross  section  contour  of  the  wire  was 
appro.xirnated  by  sixteen  linear  segments  of  equal  length.  We  not(>  that  the  restilfs 
for  B/d  =  oc-  (unshielded  structure)  are  indistinguishable  from  those  for  Bjd  =  10. 
This  is  to  be  expected,  since  the  field  of  the  bound  mode  is  tnainly  trap|)C‘d  in  the 
substrate  and  near  the  conductor.  The  presence  of  the  top  ground  plane  [)resents  a 
noticeable  disturbance  to  the  bound  wave  only  if  it  is  close  to  the  conductor  or  to  the 
dielectric  interface. 

In  Fig.  5,  we  present  dispersion  curves  for  a  three-conductor  microstrip  trans¬ 
mission  line,  which  supports  three  fundamental  modes.  The  dielectric  is  made  of 
ceramic-impregnated  teflon,  known  as  Epsilam  10.  which  is  uniaxial,  with  ti  =  1'5 
and  =  10.2.  .A.s  a  check  for  the  computer  code,  we  have  further  divided  the  di¬ 
electric  slab  into  two  layers.  Two  configurations  have  been  analyzed,  one  without  an 
air  gap  {h/d  =  0),  and  the  other  with  an  air  gap  {h/d  =  0.04).  and  the  correspond¬ 
ing  results  are  shown  in  solid  and  dashed  lines,  respectively.  The.se  configurations 
have  been  previously  analyzed  by  Kitazawa  [15].  the  first  one  by  a  full-wave  method 
(square  symbols),  and  the  second  by  a  quasi-static  approach  (dotted  lines).  It  is  of 
interest  to  note  that  a  small  air  gap  between  the  ground  plane  and  the  dielectric  slab 
results  in  big  changes  in  the  dispersion  curves.  The  longitudinal  (transverse)  current 
distributions  for  modes  1,  2,  and  3  are  found  to  be  even  (odd),  odd  (even),  and  even 
(odd),  respectively. 

In  Fig.  6,  we  present  dispersion  curves  for  a  three-wire  transmission  line  embedded 
in  a  grounded  dielectric  slab,  which  supports  three  fundamental  modes.  The  solid 
lines  represent  the  results  computed  by  Fache  et  al.  [22]  for  the  case  of  an  i.sotropic 
substrate  {Sr  =  4),  whereas  our  results  for  the  same  substrate  are  illustrated  by  the 
dotted  lines.  Note  that  a  logarithmic  scale  is  used  for  the  frequency.  Our  results  for 
a  uniaxial  substrate,  where  €t  =  3.9  and  cj  =  4.1,  are  indicated  by  three  different 
symbols.  As  can  be  seen  from  the  figure,  even  this  slight  anisotropy  has  a  noticeable 
effect  on  the  dispersion  curves.  In  the  isotropic  case,  all  three  modes  remain  in  the 
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Figure  5:  Dispersion  curves  for  a  three-strip  microstrip  transmission  line. 
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bound  regime  in  the  tre({uency  range  (■onsidere<l.  In  tin*  uniaxial  <  a'^e.  mode  2  oniiu-. 
the  leaky  regime  [36j  at  a  frequency  Iretween  d()  to  IGCIllx,  abtne  winch  the  dominant 
slab  mode,  indicated  as  TMq  in  Fig.  (>,  is  excited.  I  h<‘  dispersii>n  curve  of  the  latter 
is  obtained  by  finding  the  zero  ol  (2M)  (see  the  appendix).  1  he  longitudinal  ami 
transverse  modal  current  distributions  at  /  =  UXIIIz  are  .diuwn  in  Fig.  7.  For  modes 
1,  2.  and  2,  the  longitudinal  (transverse)  curnmts  are  «>ven  (odd).  e\en  iodd  i.  and  o<!d 
(even),  respectively.  Because  of  this  symmetry  pro[)erty.  w<-  otdy  plot  the  (urrents 
on  the  left  and  center  conductors.  We  note  that  the  longitudinal  and  trans\'e!se 
currents  are  in  phase  quadrature,  which  is  characteristic  of  bound  modes  on  lossless 
transmission  lines. 

We  next  consider  three  transmission  line  structures,  which  differ  in  the  cross  sec¬ 
tion  shape  of  the  conductors,  as  illustrated  in  Fig.  <S.  I  he  cross  sections  of  the 
conductors  are  {a)  trapezoidal  (which  may  arise  as  a  result  of  an  <‘pitaxial  growth 
process).  (6)  rectangular  (the  ideal  case),  and  (c)  inverted  trapezoidal  (which  may  be 
due  to  etching  undercuts).  The  dispersion  curves  foi  the  three  fundartiental  modes 
that  each  of  the  three  transmission  lines  may  support  are  plotted  in  Fig.  !).  Fhe 
longitudinal  (transverse)  current  distributions  for  modes  1.  2.  and  3  are  found  to  be 
even  (odd),  odd  (even),  and  even  (odd),  respectively.  In  Fig.  9.  we  also  show  the 
ciuasi-static  results  obtained  by  Schroeder  and  Wolff  [2]  for  the  same  transmission 
lines,  but  having  a  finite-width  substrate.  We  observe  that  the  dispersion  curves  for 
configurations  (a)  and  (6)  differ  less  than  those  for  (6)  and  (r).  This  is  expected, 
since  in  the  former  two  geometries  the  conductor  widths  adjacent  to  the  dielectric 
slab  (where  there  is  a  highly  concentrated  field)  are  the  same.  We  also  note  that 
the  quasi-static  results  of  [2]  are  very  close  to  the  low-frequency  limits  of  mode  3  of 
our  results.  However,  it  is  not  clear  from  [2]  to  which  mode  these  quasi-static  values 
correspond. 

In  Fig.  10,  we  present  dispersion  curves  for  a  two-strip  transmission  lint'  in  an 
unshielded  medium  comprising  both  uniaxial  and  isotropic  layers  (see  the  inset ).  Since 
there  are  only  two  conductors  and  there  are  no  ground  planes  in  this  structure,  we 
expect  it  to  have  only  one  non-cutoff  fundamental  mode  (called  mode  1  in  Fig.  10). 
Nevertheless,  an  additional  non-cutoff  mode  (called  mode  2)  has  been  found.  In 
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Figure  7:  (a)  Longitudinal  and  (b)  transverse  current  distributions  at  /  =  lOCiHz 
for  the  three-wire  transmission  line  in  the  configuration  of  Fig.  6.  with  =  3.9  and 
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Figure  8:  Geometry  of  three  transmission  line  configurations  with  conductors  of  (a) 
trapezoidal,  (b)  rectangular,  and  (c)  inverted  trapezoidal  cross  section.  The  dimen¬ 
sions  are  d  =  120  pm,  Wx  —  15  pm,  W2  =  S  —  10  pm,  and  t  =  8  pm. 


18 


151 


Figure  9:  Dispersion  curves  of  modes  I,  2,  and  3  for  the  transmission  line  configura¬ 
tions  (a),  (b),  and  (c)  of  Fig.  8.  The  quasi-static  results  of  Schroeder  and  Wolff  are 
for  a  structure  with  a  substrate  of  a  finite  width  of  130 /zm. 
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Figure  10:  Dispersion  curves  for  a  two-strip  transmission  line  in  an  unshielded  layers 
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addition  to  the  dispersion  curves  of  these*  two  modes,  we  also  plot  in  Fit;.  lU  ihr 
dispersion  curves  of  the  first  two  slab  modes  (indicated  as  I  Eo  and  IMi,  t.  computed 
from  (34)  of  the  appendix.  The  longitudinal  currents  of  mode  1  on  the  two  strij^v 
(not  shown)  are  found  to  be  in  different  directions,  whereas  those*  of  tiiode  2  are  in 
the  same  direction.  Hence,  mode  2  is  similar  to  the  fundamental  nunle  of  a  I'oated 
conducting  cylinder,  where  the  surface  current  on  the  circumference  of  the*  cylinde'r 
flows  in  the  same  direction,  and  whose  suitability  as  a  single-conducteu'  i  ransmissjon 
line  was  studied  by  Coubau  [37j.  .A  salient  feature  of  this  Coubau  mode  is  that  its 
is  very  close  at  low  frequencies  to  that  of  free  space,  and  that  its  field  is  very  loosely 
bound  to  the  dielectric.  In  fact,  when  we  replace  the  two  strips  in  the  configuration 
shown  in  Fig.  10  by  a  single  strip,  we  still  can  find  a  mode,  which  behaves  simi'arl} 
to  the  aforementioned  mode  2. 

5  Conclusion 

.A  mixed-potential  integral  equation  (MPIE)  formulation  has  been  implemented  in 
conjunction  with  the  method  of  moments  to  compute  the  propagation  constants  and 
modal  currents  of  a  multiconductor  transmission  line  embedded  in  a  laterally  open 
multilayered  uniaxial  medium.  The  approach  is  general  and  flexible,  and  can  handle 
both  open  and  shielded  structures.  It  is  applicable  to  conductors  of  arbitrary  cross 
section,  including  trapezoidal,  which  often  arises  in  practice  due  to  underetching  or 
as  a  result  of  the  epitaxial  growth  process.  .Sample  numerical  results  have  been  pre¬ 
sented  for  several  transmission-line  configurations  and.  when  possible,  compared  with 
available  published  data,  obtained  by  specialized  techniques  not  easily  extendable  to 
conductors  of  arbitrary  cross  section. 
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Appendix 

Transmission  Line  Green’s  Functions 


Consider  a  transniission-line  network  analog  of  the  layered  medium  of  Fig.  1.  form(*d 
by  a  tandem  connection  of  transmission  line  sections,  each  corresponding  to  a  di<>!ec 
trie  layer  (cf.  Fig.  2).  Let  the  network  be  e.Kcited  by  a  I  .A  current  sotirre  located  at 
z'  in  the  nth  line  section.  Then,  the  voltage  and  current  at  c  within  any  line  section, 
say  the  mth,  satisfy  the  transmission-line  equations  [32] 


r  dV^u.n 


dF 

»,mn 

dz 


—  —  7P  fP 


(32) 


where  the  propagation  constant  and  characteristic  impedance  (and  admit¬ 
tance  V’^)  have  been  defined  in  (14)  and  (15),  respectively.  The  superscript  p.  which 
stands  for  e  or  k.  as  explained  in  Section  2,  will  henceforth  be  left  out  for  brevity. 

When  m  =  n,  i.e..  the  source  and  observation  points  are  within  ihe  same  nth  line 
section,  the  voltage  V,,nn  is  readily  found  from  (32)  as  (cf.  [32.  p.  213).  [26]) 


Vt,n,.(~|c')  =  Z, 


1  4-  Tl 


(33) 


(34) 


where 

4—  . ♦ 

and  c<  =  min(c,i:'),  =  max(-,r')-  In  the  above,  F^  and  F^  are  the  vc 

reflection  coefficients  “looking  to  the  left”  and  “looking  to  the  right."  respectively,  at 
the  two  interior  ends  of  this  line  section,  as  illustrated  in  Fig.  11.  I'hese  reflection 
coefficients  can  be  found  as 

Z„  -  Zn 
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(35) 
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with  the  terminal  impedances  given  by  the  recursive  relations 


—  ^n±l 


^n±\  '^jZn±l  tan  (  ^■;.n±l  dn±i  ) 
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)  I 


~1“  J  tan  (  A ^ I  dfj ^  j ) 

where  the  upper  and  lower  signs  correspond  to  the  right  and  left  arrows,  respectively. 
For  m  ^  71.  V’,nn(r|r')  is  readily  found  from  (32)  and  (33).  by  (Uiforcing  the- 
continuity  of  the  voltage  and  current  at  the  interfaces.  As  a  result,  we  obtain  ;J(i] 


VU.(c|y}  =  { 
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i=m+i  1  +  p 


It  is  understood  in  the  above  that  the  product  terms  are  equal  to  one  if  the  lower 
limits  exceed  the  upper  limits. 

The  transmission-line  Green's  functions  can  be  efficiently  implemented  into  a  com¬ 
puter  program,  as  explained  below.  First,  we  recognized  that  (33)  may  be  written 
as 

V^„n(^ic')  -  Z„  /i  (n;  c;  Z;  T,;  (40) 

which  serves  to  define  the  function  fi.  The  corresponding  current  can  then  easily  be 
found  from  the  first  of  (32)  as 

>  , 


h,Uz\z')  =  ±f2{r7-z-,z':±%-.±'rny  =^z'  (41) 

which  defines  the  function  In  a  like  manner,  we  may  write  (37)  as 
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which  defines  the  functions  and  i'he  corresponding  current  can  again  Ik*  ionini 
from  the  first  of  (32)  as 


^  t .  m  n  I  1  ; 


:  ^  ("ah  .  f  I  .  ^  —  T'j  j  I  <  ^  <  r;  ) 


rn  >  /) 
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(  J3i 


The  current  /,,,mn  and  voltage  due  to  a  unit-strength  voltage  source  in  the 

nth  line  section,  satisfy  a  set  of  equations  dual  to  (32).  Hence,  vve  may  obtain  these 
voltage  and  current  transmission-line  Green's  functions  from  (10)  (43)  by  tnaking 
the  substitutions  V  — *•  /,  I  — ►  V’,  and  Z  — *•  Y.  Note  that  the  last  substitution 
causes  all  reflection  coefficients  to  change  signs.  We  observe  that  only  four  subrou¬ 
tines,  corresponding  to  the  functions  /i  through  /,.  are  required  to  implement  all  the 
transmission-line  Green's  functions  in  the  computer  code. 
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